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Final  Report:  Grant  No.  DAAL  03-89-G-0039 
XIII.  International  Workshop  on  Condensed  Matter  Theories 

The  grant  enabled  10  scientists  from  the  U.S.A.  to  participate  at  the  XIII. 

International  Workshop  on  Condensed  Matter  Theories  held  at  Campos  do  Jordao,  Brazil  in 
August  1989.  Without  this  grant  the  U.S.  Participation  at  the  Worshop  would  have  been 
very  minimal.  The  grant  allowed  U.S.  scientists  to  interact  effectively  with  their 
international  counterparts. 

The  list  of  the  scientists  from  the  U.S.A.  who  have  been  supported  by  the  grant  is 
attached  in  Appendix  1  and  the  titles  of  their  talks  presented  at  the  Workshop  are  listed  in 
Appendix  2.  Participants  from  the  U.S.A.  receiving  the  grant  have  acknowledged  that  in 
their  respective  articles. 

Copies  of  the  articles  containing  the  subject  matter  presented  by  the  participants 
supposed  by  the  grant  are  attached  in  Appendix  3.  All  presentations  were  invited  talks. 
These  articles  have  been  accepted  for  publication  in  Condensed  Matter  Theories  Vol.  5.  The 
publisher,  Plenum  Publishing  Corporation,  expects  to  bring  out  the  volume  shortly. 

In  keeping  with  the  objective  of  the  grant,  emphasis  was  placed  on  theories  related  to 
high  temperlure  superconductivity,  molecular  dynamics  and  physics  of  strongly  correlated 
systems.  In  addition,  talks  on  the  use  of  maximum  entrophy  principle  to  describe  cooperative 
phenomena  were  presented.  A  large  number  of  scientists  from  the  host  country,  Brazil, 
participated  actively.  Two  topmost  scientists  from  the  U.S.S.R.,  Professor  Bashkin  and 
Pitaevski  of  their  National  Academy  of  Science’s  Institute  for  Physical  Problems,  fook  part. 

The  Workshop  also  helped  to  foster  another  key  objective  of  bringing  together 
scientists  working  in  different  areas  of  physics  in  order  to  fasrilitate  interchange  of  ideas. 


The  Workshop  was  very  effective  in  cross  fertilization  of  ideas  and  development  of  new 
interests.  For  example,  Dr.  Proto,  originally  a  nuclear  physicist,  presented  the  talk  on  the 
role  of  maximum  entrophy  principle.  Dr.  Campbell,  originally  a  solid  state  theorist, 
presented  the  talk  on  correlations  in  atoms.  The  principle  investigator,  originally  an  atomic 
and  nuclear  physicist,  presented  a  talk  on  superconductivity.  Another  objective  of  the 
workshop,  namely  to  serve  as  a  forum  to  initiate  collabortive  research  in  the  area  of 
condensed  matter  theories,  has  also  been  fulfilled. 

Thus,  the  purpose  and  objective  of  the  grant  have  been  successfully  achieved. 
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THOMAS-FERMI  EQUATION  OF  STATE-THE  HOT  CURVE 


George  A.  Baker,  Jr.  and  J.  D.  Johnson 

Theoretical  Division,  Los  Alamos  National  Laboratory 
University  of  California,  Los  Alamos,  N.  M.  87545,  USA 


ABSTRACT 

We  derive  the  high-temperature  limit  of  the  equation  of  state  based  on  the 
Thomas-Fermi  statistical  theory  of  the  atom.  The  resulting  “hot  curve”  is  in  fact 
the  ideal  Fermi  gas.  We  expand  the  thermodynamic  properties  of  this  gas  in  powers 
of  the  fugacity  and  use  this  expansion  to  construct  a  representation  of  the  pressure, 
accurate  to  about  0.1  %.  This  representation  is  compared  with  the  actual  theory  for 
aluminum  and  the  “hot  curve”  is  found  to  represent  it  well  over  a  large  region  of 
interest  in  applications. 


1.  INTRODUCTION  AND  SUMMARY 


The  Thomas-Fermi  (T-F)  statistical  theory  of  the  atom1  as  well-as  the  modifica¬ 
tions  due  to  Dirac2  have  long  been  used  as  a  basic  starting  point  for  the  computation 
of  approximations  to  the  equations  of  state.3,4  In  order  to  maJke  use  of  this  proce¬ 
dure,  computer  programs  have  been  written  to  compute  the  nunv  rical  content  of  the 
theory.  They  consume  a  sufficient  amount  of  computer  time,  even  today,  so  that  it 
is  impractical  to  use  them  to  compute,  ab  initio ,  the  value  of  the  pressure,  internal 
energy,  eic.,  every  time  that  a  new  value  is  required  inside  an  application  computer 
program.  Besides,  as  these  efforts  represent  only  approximate  equations  of  state,  some 
adjustment  is  necessary  to  bring  them  into  accord  with  physical  reality.  Consequen¬ 
tially,  to  date  largely  empirical  fits  have  been  used  to  represent  the  equations  of  state 
for  the  purposes  of  applications. 

In  this  work,  we  are  concerned  with  beginning  an  analysis  of  the  physical  struc¬ 
ture  of  the  equations  of  state  of  real  matter.  As  a  start,  we  will  study  the  Thomas- 
Fermi  model  equation  of  state  which  represents  a  fair  amount  of  the  physics,  at  least 
in  some  regions.  One  method  which  is  normally  fruitful,  is  to  consider  various  limits. 
There  are  currently  two  which  are  known.  The  first  is  the  low-density  limit.  Here 
there  is  complete  ionization  when  the  system  is  in  equilibrium  and  the  pressure  for 
an  element  of  nuclear  charge  Z  is 

PCl/N  =  (Z  +  1  )fcT,  (1.1) 

the  ideal  gas  equation  of  state.  Here  P  is  the  pressure,  Q  is  the  volume  of  the 
system,  N  is  the  number  of  atoms,  k  is  Boltzmann’s  constant  and  T  is  the  absolute 
temperature.  The  second  limit3  is  the  low- temperature  limit,  or  the  “cold  curve.” 
Here  the  pressure  is  of  the  form, 

PSl/N  =  Z*<t>(ZQ/N)t  (1.2) 

where  6(x)  is  a  well  defined  function.  If  we  think  of  the  temperature-density,  quarter- 
plane,  these  results  give  the  limiting  behavior  of  the  T-F  model  along  the  zero- 
temperature  and  the  zero-density  edges.  There  remain  the  high-density  and  the 
high-temperature  regions  to  examine  for  physical  structure. 

One  might  think  that  in  the  high- temperature  limit  it  would  be  appropriate  to 
describe  the  system  in  purely  classical  terms.  Indeed  if  such  were  the  case,  Baker5 
has  proven  that  the  pressure  would  be  of  the  form, 

PCl/N  =  kT  f(ClT3/N ,  Z).  (1.3) 

The  Debye-Huckel  correction6  is  of  just  this  form.  Also  Baker  has  shown  for  this  case 
that  the  internal  energy  has  the  particularly  simple  form, 

u  =  3PQ  -  | (Z  +  1  )NkT.  (1.4) 

The  statistical  mechanics  of  Coulombic  systems  have  been  much  studied.7  It  is  now 
well  known  that  there  does  not  exist  a  classical  ( i.e .  Planck’s  constant  h  =  0)  gas 
because  atoms  with  a  Coulomb  interaction  collapse  to  E  =  — oo.  Thus  if  we  are  to 
ever  introduce  a  Coulomb  attraction  between  the  atomic  nucleus  and  the  electrons, 
we  must  necessarily  include  some  account  of  the  quantum  effects  that  are  needed 
to  stabilize  the  system.  As  is  also  well  known  there  are  two  important  physical 
lengths  to  be  considered.  The  first  is  the  de  Broglie  length  which  is  proportional 

to  h/ymkT ,  where  m  is  the  electron  mass,  and  which  measures  in  a  noninteracting 
gas  the  importance  of  quantum  effects.  The  Coulomb  interaction  does  not  by  itself 
piovide  the  second  length  and  the  difficulty  of  its  long  range  can  not  be  circumvented 
by  studying  dilute  systems  because  it  contains  no  parameter  with  the  dimensions  of 

i 
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a  length.  The  second  length  is  the  Debye  screening  length  which  is  proportional  to 
c2/kT.  This  length  is  however  a  statistical  cflcct.  and  should  follow  from  the  theory, 
but  unfortunately  is  not  there  ah  initio.  Thus  when  we  look  to  the  high-temperature 
and  high-density  regions,  if  we  consider  the  cases  where  Q./N  »  ( c2/kT )3,  then  wc 
can  hope  to  start  with  a  noninteracting  election  gas  (with  a  background  gas  of  atomic 
nuclei)  as  the  basic  system. 

In  the  second  section,  we  derive  the  limit  of  Thomas-Ferim  theory  when  the 
Debye  screening  length  is  negligible  compared  to  the  interparticle  distance,  and  the 
de  Broglie  length  remains  arbitrary.  We  find  that  it  correctly  reduces  to  the  ideal 
Fermi  gas.  We  call  this  limit  the  “hot  curve,”  because  it  is  readied  if  one  either  fixes 
the  density  and  lets  the  temperature  go  to  infinity,  or  mudi  less  restrictively,  it  is  also 
readied  if  one  fixes  the  de  Broglie  length  and  then  lets  the  temperature  go  to  infinity. 
In  the  third  section  we  review  the  theory  of  the  ideal  Fermi  gas  and  describe'  how 
to  calculate  its  properties  in  a  practical  manner.  We  derive  lengthy  fugacity  series 
and  find  that  the  pressure  function  can  be  approximated  to  within,  say  0.1%,  by  a 
low-order,  two-point  Pade  approximant.  In  the  final  section  we  compare  the  ideal  gas 
approximation  to  results  for  aluminum  and  map  out  its  region  of  validity  to  various 
degrees  of  accuracy. 


2.  HIGH  TEMPERATURE  LIMIT  OF  THOMAS-FERMI 
THEORY 

Thomas-Fermi  theory  has  been  applied  to  compute  equations  of  state  at  finite 
temperature  by  Feynman  et  al.z  They  begin  with  an  application  of  the  statistical 
analysis  of  Fermi  and  Dirac  which  leads  to  the  equation 


2  •  Airp2dp/hz 

exp[(p2/2m  —  eV)/kT  +  ?/}  -f-  1  ’ 


(2.1) 


where  —eV  is  the  potential  energy.  We  follow  them  in  defining  for  convenience  the 
auxiliary  functions 


(2.2) 


Then  one  uses  Poisson’s  equation  to  determine  V(r).  It  yields 


i£_ 

r  dr 2 


(rV(r)) 


167T2 


e(2  mkT)^Ii 


(2.3) 


Note  that  in  the  case  of  no  interaction  that  the  right-hand  side  of  (2.3)  vanishes  (e=0) 
and  so  the  equation  implies  that  V  =  a  +  b/r  where  a  and  b  are  constants.  In  order  to 
simplify  the  above  equation,  Feynman  et  al?  introduce  dimensionless  variables.  First 
they  define  a  length  scale, 


hz 


c  = 


K32z2e2m(2mkT)* 


(X  T~4, 


where  s  =  r/c.  Then  since  t?  is  independent  of  r ,  (2.3)  becomes 


ds2 


where 

/?/«  =  (eV(r)/kT)  - 


(2.4) 

(2.5) 

(2.6) 


t 


The  boundary  conditions  of  (2.5)  become,  as  at  the  origin  V(r)  must  behave  as  Ze/?-, 


fi{0)  =  n  =  ZoljkTc  <x  T~*.  '  (2.7) 

The  scheme  employed  is  to  suppose  that  each  atom  is  confined  to  a  sphere  of  volume 
equal  to  the  volume  per  particle.  This  is  clearly  an  approximation.  The  other  bound¬ 
ary  condition  is  to  require  that  the  number  of  electrons  in  the  sphere  is  exactly  equal 
to  the  nuclear  charge.  A  little  manipulation  serves  to  show  that  the  condition, 

=  fi/s  at  s  =  5,  (2.8) 

as 

imposes  this  normalization  in  the  sphere  of  radius  r  =  cb.  Feynman  et  al.3  derive, 
among  other  things,  the  formula  for  the  pressure  as 

PSl/N  =  \(ZkT)  ■  (A)  .  (2.9) 


where  fib  is  the  value  of  fi  on  the  boundary  s  =  b. 

In  a  parallel  way  we  may  set  out  the  corresponding  formulae  for  the  ideal  Fermi 
gas.  In  this  case  the  electron  density  is  simply  given  by  (2.1)  with  e  =  0.  As  rj 
is  independent  of  r,  one  sees  immediately  by  (2.6)  that  the  equation  for  the  density 
(2.5)  is  simply  satisfied.  Since  by  (2.4)  and  (2.7)  both  the  length  and  magnitude  scales 
depend  on  the  electronic  charge  e  =  0,  the  normalization  equation  (2.8),  in  leading 
order,  is  automatically  satisfied,  and  so  does  not  determine  the  number  of  electrons 
in  this  limit.  Returning  to  (2.1),  we  may  impose  the  normalization  condition  by 
integrating  the  density  over  a  sphere  of  radius  r.  It  gives 


(2.10) 


which  implies  r\.  In  this  limit,  the  pressure  equation  (2.9),  becomes, 

pa/N  =  \{ZkT)  QQ  /,(-,).  (2.11) 

a  parametric  expression  for  the  pressure  in  terms  of  the  77  of  (2.10).  Note  is  made 
that  cza  is  independent  of  the  electronic  charge  e  =  0,  so  tills  form  is  valid  in  this 
noninteracting  limit.  Comparison  with  the  results  of  Huang8  for  the  ideal  Fermi  gas, 
reveal  complete  agreement,  when  it  is  remembered  that  for  our  case  the  spin,  s  = 


Now  we  are  ready  to  consider  the  “hot  curve”  limit  of  the  Thomas-Fermi  theory. 
In  the  basic  equations  of  the  theory,  (2.5,  7-8),  we  make  the  following  change  of 
variables, 


a  =  s/ a* ,  j  =  fi/a^ . 

(2-12) 

We  thus  obtain 

(2.13) 

CNj« 

II 

O 

V 

(2.14) 

d'y  7 

-j--  =  — ,  at  the  boundary. 

(2.15) 

lu  the  limit  a  — >  0  (by  (2.7)  this  limit  is  equivalent  to  T  -•*  oo),  we  obtain  the  result 
that  7  =  Aa  solves  (2.13-15).  Again,  as  at  (2.10)  above,  wc  have  an  undetermined 
normalization  constant  to  be  determined  because  m  our  high-temperature  limit  (2.15) 
is  satisfied  automatically.  Again  referring  to  (2.1)  we  obtain  the  normalization  con¬ 
dition. 

2  =  (I)  .  (2.10) 

which  determines  A  and  thus  the  solution  of  (2.13-15).  When  we  note  the  comparison 
,4  =  — 77,  we  find  that  this  limiting  solution  is  the  same  as  the  one  we  obtained  for 
the  ideal  (noninteracting)  Fermi  gas.  This  result  completes  our  demonstration  of  the 
proposition  that  the  “hot  curve”  for  Thomas-Fermi  theory  is  the  ideal  Fermi  gas! 

3.  PROPERTIES  OF  THE  IDEAL  FERMI  GAS 

The  basic  theory  of  the  ideal  Fermi  gas  is  d.  abed  by  Huang.*  To  establish  a 
correspondence  between  the  results  of  the  previous  section  and  more  standard  nota¬ 
tion,  we  note  that  in  ,(2.16)  7 ja  —  A\  therefore  we  introduce  the  notation  z  —  e“  . 
We  can  then  rewrite  (2.16)  and  (2.11)  as 


/  2? rmkT\7  2  f00  zy?e  ydy 

\  h2  )  s/x  Jo  1  +  2' 

27rmkT\3  4  f°°  zy?c~ydy 


h 2  )  s/k  . 
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f  ZJ 
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where  P  is  the  pressure  due  to  the  electrons  only  and  does  not  take  account  of  the 
effect  of  the  motion  of  the  center  of  mass  of  the  atom.  If  we  introduce  the  further 
notation, 

A-(s^sV.  (M) 


27t  mkT 


2  [~zyU-'iy  f-(-l  )'+'z' 

AM -Try.  tt^  =  L— —  • 


— .1V+1*' 


zy*e  ydy  ^  (  — 1  ),+1zl 


r,x  4  f  zy*e~ydy  ^ 
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where  the  series  expansions  are  convergent  for  \z\  <  1.  We  may  now  rewrite  (3.1-2) 


ZN  A3 


-  =  /»w, 


AA.-IM  (37) 

ZNkT  /$(*)'  K  } 

where  £  is  the  de  Broglie  density.  The  procedure  to  calculate  the  pressure  of  the  ideal 
Fermi  gas  is  now,  in  principle,  quite  straightforward.  Eq.  (3.6)  is  solved  for  z  and 
then  that  value  is  substituted  into  (3.7). 

To  evaluate  these  expressions  numerically  we  choose  the  following  method.  First 
we  revert  the  series  expansion  (3.6)  to  give  z(£)  as  a  series  in  £.  Then  we  substitute 
it  into  (3.7)  to  obtain 

PQ. 

ZNkT  =  S®'  ( 3 ’8) 


We  have  calculated  the  leading  3G  terms  of  the  series  expansion.  The  method  used  is 
the  classical  Lagrange  formula  for  the  reversion  of  series.9  The  only  point  of  difficulty 
is  that  a  large  number  of  decimal  places  are  lost  in  the  computation  in  this  case.  We 
ha  'f  therefore  taken  the  precaution  of  using  at  least  58  decimal  places  to  carry  out 
these  computations.  The  results  are  listed  in  Table  1. 

TABLE  1.  ( PQ/ZNkT )  as  a  series  in  the  do  Broglie  density 

0  1.0000000000  0000000000  0000000000  0000000000  00000000000000000E+000 

1  1.7677669529  6636881100  2110905262  1225982120  89844221 18509 147E-001 

2  -3.3000598199  1683655758  8617889323  8790328003  8917U3930S782E-003 

3  1.1128932846  6542504524  9253533917  1305775999  1875768224181 E-004 

4  -3.5405040951  9736538278  3050093233  4626176046  46439677965E-006 

5  8.3863470395  69257296i9  7125843681  6218474298  42743624SE-008 

6  -3.6620617873  4852703663  1688233937  9045807824  8643167E-010 

7  -1.0280607154  3957929799  3273512206  9735581999  5254513E-010 

8  7.0550978435  72S3454626  0275709452  8261969773  091S8E-012 

9  -2.6859639507  9285424406  0526716388  7926863588  4377E-013 

10  4.0571834908.0612166197  1056127182  3091151601  35E-015 

11  2.7970439770  9162019148  3071234746  1358106846  6E-016 

12  -2.8379673439  5952590529  6631787032  9726025304  E-017 

13  1.3992940717  5922219970  7552151122  203412696E  -018 

14  -3.6303052861  0821033013  0082398676  2418074E-0  20 

15  -6.0257400821  7251347692  8112664253  67093E-022 

16  1.2989538153  2549763684  7035089386  73544E-022 

17  -8.1719971340  6344259697  7319803759  795E-024 

18  2.9413082494  4946667164  3606073469  73E-025 

19  -2.0285711098  2088612486  4658243931  E-027 

20  -5.7410636166  7615749309  984730023E  -028 

21  4.8461575378  3763503589  33968480E-  029 

22  -2.2369786852  5871386652  1846940E-0  30 

23  4.7888680538  7474310454  78772E-032 

24  2.0304880286  8391265410  8553E-033 

25  -2.7811009124  7360566430  414E-034 

26  1.6149810555  1163427972  12E-035 

27  -5.2554355032  5730228297  E-037 

28  -1.3309033541  33284697E-  039 

29  1.4721238409  86015824E-  039 

30  -1.1062516681  9956070E-0  40 

31  4.7267873838  86169E-042 

32  -7.6386716803  536E-044 

33  -6.5324794996  62E-045 

34  7.1193401844  5E-046 

35  -3.8268661579  E-047 

36  1.097950074E  -048 


The  above  series  expansion  was  derived  for  \z\  <  1,  but  the  above  series  plainly 
corresponds  to  a  larger  range.  In  the  limit  as  z  — >  co  Huang  shows  that 

AW* j^*0**)1  [I  +  i(i^F  +  '”]  +0(2~')-  <3'9) 

From  the  identity,8  zjzfz(z)  —  fz(z)  one  can  easily  also  derive  the  asymptotic  be¬ 
havior  of  f&(z),  and  thus  from  (3.7)  the  asymptotic  behavior  of  g(Q-  We  obtain, 


I 

i 


(3.10) 


With  this  information  and  the  series  of  Table  1,  we  may  construct  a  two  point 
Pade  approximant10  to  [ <7(01 3  °f  the  form  (Ar  +  2 /N]  which  is  exact  through  order 
£2A’+1  a(  (|1C  origin,  and  is  also  asymptotically  correct  as  £  — *  oo.  We  find  excellent 
cor  •  ,  yen ce  for  this  method  and  that  for  0  <  £  <  oo  we  get  an  accuracy  of  about  0.1 
pei  for  57(C)  from  the  approximation, 


9(0 


1  -1-  0.010948S0C  +  0.12660436£2  +  0.0091177644£31 3 
1  +  0.08061S739C 


Thus  the  total  pressure  would  be  (including  the  center  of  mass  motion) 


Nk-T 

P  =  —{1  +  Zg(0}- 


(3-11) 


(3.12) 


In  the  case  where  the  temperature  is  fixed  and  — >  00,  the  low-density  limit,  not  only 
does  the  Debye  density  go  to  zero,  as  required  to  obtain  the  ideal  Fermi  gas  limit  of 
Thomas-Fermi  theory,  but  also  £  — *  0.  In  this  case,  as  g(0)  =  1,  (3.12)  reduces  to  (1.1) 
and  thereby  supplies  an  alternate  derivation  of  the  low-aensity  limit  of  Thomas-Fermi 
theory. 

As  Huang8  points  out,  the  internal  energy,  U ,  for  this  case  follows  simply  from 
(3.12)  as, 

U  =  |  Pa.  (3.13) 


Epstein11  shows  from  the  thermodynamic  relation  dS  =  (dU  +  PdV)/T ,  the 
above  results,  and  Nernst’s  heat  postulate  that  the  entropy  of  the  ideal  Fermi  gas  is 
simply  given  by 

s.  =  ZNk  (|S(C)  -  logz(o)  .  (3.14) 


where  the  limit  as  T  — *  0  is  the  limit  £  — >  00  by  (3.6)  and  as  Epstein  further  points 
out  Sc  —*  0  in  this  limit.  If  we  add  the  contribution  of  the  motion  of  the  center  of 
mass  to  the  entropy,  we  get 


5  =  Nk 


-(z + 1)  iogc + 1 + z  (|j(o  -  iogM<)/<)' 


-[-constant,  (3.15) 


The  Helmholtz  free  energy  is  now  given  directly  by  A  =  U  —  IS.  The  Gibbs  thermo¬ 
dynamic  potential  is  also  directly  given  and  is  G  =  U  —  TS  +  PSL. 

It  now  remains  to  give  a  representation  of  log  z(0  =  log  C  +  log[.z(C)/C] 
complete  the  representation  of  the  thermodynamic  quantities  for  the  ideal  Fermi  gas. 
Since  logr  x  £3,  the  problem  of  deriving  a  representation  for  log[jr(C) / C]  should 
be  similar-  to  that  of  the  representation  (3.11).  We  give  in  Table  2  the  necessary 
series  coefficients  in  £  for  z(()  to  work  on  this  representation,  but  we  will  leave  it  for 
the  future.  Thermodynamic  consistency  depends  on  the  equation  between  the  two 
representations 

S(0  +  C*'(C)  =  C^1 

d( 


1 


(3.16) 


TABLE  2.  The  fugacity  r  as  a  series  in  the  de  Broglie  density 


1  1.0000000000  0000000000  0000000000  0000000000  00000000000000000E+000 

2  3.5355339059  3273762200  4221810524  2451964241  7968844237018294E-001 

3  5.7549910270  1247451636  1707316601  4181450799  416243291041327E-002 

4  5.7639604009  1025440341  8852781947  0758923518  58214221729707E-003 

5  4.0194941515  2300959555  6172119656  7773364832  0998466829345E-004 

6  2.0981898872  2604799054  4860297423  5099614729  957102872728E-005 

7  8.6021310842  6030566004  3913343164  3181688359  0277772573E-007 

8  2.8647148623  7664872936  8242210245  0573640824  266032220E-008 

9  7.9528314678  5241689019  4817612245  1032872937  5035650E-010 

10  1.0774425910  0567756220  4988130993  7541387605  437996E-011 

11  3.8247968264  1809029592  4653344686  7070280382  2264E-013 

12  6.8432943010  1907998578  8027623030  3596055059  29E-015 

13  1.0762104093  0537917245  5417733813  6774703889  3E-016 

14  1.5124110216  1988369105  9052478125  978137640E  -018 

15  2.0715738792  9770436279  3713783632  7032961E-0  20 

16  1.3846671521  9900108771  8574969994  14568E-022 

17  5.3288541784  7605238410  1301497951  755E-024 

18  3.5079551301  2368023505  6432045696  E-027 

19  -5.9656175104  9257472195  3065263300  E-027 

20  5.2969138512  2627670501  874181389E  -028 

21  -2.5226985875  2718441504  10473445E-  029 

22  6.0209616883  8744484633  512535E-03  1 

23  1.8543035351  4383646428  76522E-032 

24  -3.0176817670  7158240262  6353E-033 

25  1.8757233170  6238133052  809E-034 

26  -6.7714760730  2256395698  9E-036 

27  3.7182598930  255841378E  -038 

28  1.4954203444  742341364E  -038 

29  -1.2728642729  9966' )53E-  039 

30  6.0377265821  589225E-04  1 

31  -1.3644496192  99721E-042 

32  -5.3539191733  757E-044 

33  7.8650740191  78E-045 

34  -4.7690907071  0E-046 

35  1.6535692458  E-047 

36  -2.3890246E-0  50 

37  -4.2646358E-0  50 


An  alternate  procedure  would  be  to  determine  z(()  directly  from  this  equation 
subject  to  the  boundary  condition  lim^_o  z(C)/£  =  •  •  This  equation  is  an  identity  in 
the  exact  theory  and  not  an  extra  condition. 

From  the  theoretical  point  of  view  the  most  satisfactory  proceedure  would  be  to 
construct  a  sufficiently  accurate  representation  of,  say,  the  Helmholtz  free  energy  A 
that  would  provide  adequately  accurate  derivatives  =  —  P,  and  (44)«/  =  —S. 

Using  (3.13),  (3.15)  (ignoring  the  constant),  and  integrating  (3.16)  we  have  for  the 
Helmholtz  free  energy, 

A  =  -  PQ  +  (Z  +  l)AU;r  log  C  +  ZNkTXo&iQ/C,) 

=NkT  ( Z  +  l)(logC  -  1  )  +  Z  l\g(n)  -  1)^1  ,  (3.17) 

Jo  V 

for  which  the  series  expansion  in  (  can  be  easily  derived  from  Table  1.  The  inability 
to  assign  an  absolute  entropy  for  the  ordinary  ideal  gas,  leaves  A  uncertain  by  a 

i 

j 


linear  term  in  T.  It  remains  to  do  seen  which  of  t.ho  pioccdures  outlined  above  are 
computationally  most  efficient. 


4.  COMPARISON  OF  IDEAL  FERMI  GAS  TO 
THOMAS-FERMI  THEORY 

We  now  show  the  extent  of  agreement  for  aluminum  between  the  ideal  Fermi  gas 
and  the  Thomas-Fermi  theory.  We  use  the  computer  program  of  D.  A.  Liberman12 
to  compute  the  T-F  numbers.  We  present  the  results  in  the  figures  as  contours  of 
percentage  differences  (electron  properties  only). 

For  the  pressure,  Figure  1  shows  in  temperature-density  parameter  space  the 
1%.  10%,  and  30%  contours,  as  one  goes  from  the  top  curve  of  the  figure  to  the 
bottom,  respectively.  The  expected  feature  is  that  for  high- temperature  and/or  low 
density  the  ideal  gas  is  accurate.  The  10%  contour,  for  example,  will  serve  as  our  “hot 
envelope,”  that  is  to  say,  the  limit  of  the  validity  of  the  “hot  curve”  approximation.  For 
low-temperature  and  high-density  the  ideal  Fermi  gas  is  again  a  good  representation 
of  the  T-F  theory  because  the  electrons  are  being  forced  to  the  pressure-ionized, 
degenerate,  free  electron  gas.  Since  as  the  density  increases  the  kinetic  energy  per 
atom  is  forced  by  the  Pauli  principle  to  increase  proportional  to  the  density  to  the 
two-thirds  power  (relativistic  corrections  are  ignored  here)  and  the  potential  energy 
is  expected  to  increase  only  as  the  one-third  power  of  density,  the  free-electron-gas 
energy  becomes  dominate.  This  effect  is  begining  to  be  evident  in  the  behavior  of  the 
30%  contour.  The  ranges  of  temperature  and  density  shown  are  those  of  interest  for  a 
great  many  applications.  Thus  the  ideal  Fermi  gas  well  reproduces  the  T-F  pressure 
over  a  substantial  region. 


Figure  1.  Pressure  contours. 
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Figure  2.  Energy  contours. 

Figure  2  shows  the  results  for  the  internal  energy.  Here  we  see  only  the  10% 
and  30%  contours  because  the  ideal  Fermi  gas  does  not  represent  the  T-F  energy  as 
well  as  it  does  the  pressure.  This  result  is  at  least  partly  due  to  what,  in  effect,  is 
an  extra  term  present  in  the  T-F  energy  and  not  in  the  T-F  pressure.  The  bound 
electrons  do  not  contribute  to  the  pressure  but  do  have  a  large  effect  on  the  energy,  for 
the  temperature  and  density  both  small.  Since  the  free  gas  has  no  bound  electrons, 
there  is  more  difficulty  in  matching  the  T-F  energy.  However,  there  is  again  a  “hot 
envelope.” 

We  did  one  other  study  that  was  beyond  our  original  intent.  Our  goal  is  really 
not  to  find  an  analytic  representation  of  the  T-F  theory,  but  to  obtain  a  fit  to  the 
T-F  with  the  zero- temperature  isotherm  subtracted.  Thus  it  is  of  interest  to  compare 
just  such  a  result  to  the  ideal  gas  with  its  zero-temperature  isotherm  subtracted.  We 
expect  an  even  better  correspondence  between  these  pressures,  with  exact  agreement 
both  at  low-density/high-temperature  and  zero  temperature.  Figure  3  shows  again 
the  1%,  10%,  and  30%  contours  for  pressure  and  indeed  there  is  improvement  over 
Figure  1  with  the  “hot  envelope”  now  at  lower  temperatures.  We  do  not  show  the 
contours  that  appear  at  low  temperature  as  they  are  not  of  interest  to  us  in  this 
study.  The  odd  vertical  steps  arise  because  really  the  two  contours  at  that  point 
loop  back  under  themselves  and  come  back  to  the  lower  curves  due  to  the  forced 
agreement  at  zero  temperature.  We  did  not  put  in  these  loops  because  we  felt  that 
was  a  misrepresentation  of  the  high- temperature  behavior. 

The  energy  contours  with  zero-temperature  isotherm  subtracted  are  not  pre¬ 
sented  because  the  results  did  not  turn  out  as  well  as  for  the  pressure.  This  result  is 
again  caused  by  the  absence  of  the  bound  state  energy  in  the  free  Fermi  gas. 

In  general  we  see  the  “hot  envelope”  and  reasonable  agreement  between  the  free 
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Figure  3.  Pressure  contours  for  the  zero  temperature  isotherm  subtracted. 

Fermi  gas  and  T-F  theory  for  a  large  region  of  pressure.  We  understand  the  difference 
between  the  pressure  and  internal  energy. 

This  work  was  performed  under  the  auspices  of  the  U.  S.  Department  of  Energy. 
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from  the  U.S.  Army  Research  Office  and  from  Fundunesp  to  permit  his  attendence  at 
the  XIII  International  Workshop  on  Condensed  Matter  Theories. 
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INTRODUCTION 

Theories  of  electron  correlations  in  condensed  matter  systems  with  a  high  degree  of  density  vari¬ 
ation  and/or  electron  localization  may  be  tested  in  many-electron  atoms  and  ions,  where  nearly  exact 
correlation  energies  are  known.  We  have  adapted  the  inhomogeneous,  non-local  Feenberg-Jastrow  Euler- 
Lagrange  theory  developed  by  Krotscheck  and  his  collaborators1,2  to  apply  to  atomic  systems.  Correla¬ 
tion  energies  for  four  electron  and  ten  electron  closed  shell  atoms  (Be  and  Ne,  respectively)  are  in  good 
agreement  with  the  known  results.  We  also  calculate  the  pair  correlation  functions,  which  are  found  to 
exhibit  an  extreme  sensitivity  to  the  location  of  the  pair  of  electrons  within  the  atom. 

In  the  next  section  we  define  the  theory  and  the  approximations  which  are  necessary  within  the 
theory  A  more  complete  discussion  of  this  theory  in  the  context  of  other  inhomogeneous  quantum  fluids 
is  given  by  Krotscheck  et  al.  elsewhere  in  this  volume. 

The  third  section  contains  the  results  of  this  work  for  Be  and  Ne.  We  conclude  with  a  brief 


discussion. 


FEENBEKG-JASTKOW  EULEH- LAG  RANGE  THEORY  OF  INHOMOGENEOUS  QUANTUM  FLU¬ 
IDS 


The  Hamiltonian  for  an  inhomogeneous  system  is 

11  ~  H  V'?  +  ]C  +  t,2(r'V).  0) 

t  i  »<; 

where  vj(r)  is  the  external  potential,  v2(r,j)  is  the  two-body  potential,  and  N  is  the  number  of  electrons 
in  the  system.  In  the  present  case,  vi(r)  =  —  Z/r  where  r  is  the  distance  to  the  nucleus,  and  t>2(r)  =  1/r, 
where  we  use  atomic  units  (ft  =  c  =  m  a  1),  and  the  nucleus  is  taken  to  have  an  infinite  mass. 

The  Jastrow-Feenberg  wave  function  space  is  defined  by: 


|tf)  =  expji  ^ui(r,)  +  ^ti2(rf,rj)  |  |4>) 


where  |<I>)  is  a  Slater  determinant  of  one  body  orbitals  which  are  mutually  orthonormal,  and  ui(r)  and 
ti2(r;,r;*)  are  the  one-  and  two-body  Jastrow  pseudopotentials.  Each  of  these  functions  is  determined 
by  functional  variation  to  minimize  the  energy  expectation  value,  E: 


6E  __  6  (gjgjg) 

6 tin  ”  6u„  ($1*)  " 


(nal.2) 


(3) 


The  optimal  depends  solely  on  the  choice  of  the  single-  particle  orbitals  xpi  which  which  are  solutions 
of  the  correlated  Hartree-Fock  equation: 


[_Iv2  +  Uext  (r)  +  V//(r)]^,-(r)  +  J  dir,Vg(r,r')p1(r,r,)ipi(r')  =  €;^;(r), 


where  ^//(r)  is  the  generalized  Ilartree  potential,  V^(r,r')  is  the  exchange/correlation  interaction,  and 
Pi( r,r')  is  the  one  body  density  matrix  of  the  wave  function. 

The  two-body  Euier-Lagrange  equation  involves  the  two-body  Jastrow  pseudopotential  ti2(r»,i;) 
and  the  pair  distribution  function  <7(r,r'),  defined  by 


P(rs,r2)^i(ri)/?i(r2)  =  N(N  -  1) 


/d3r3...(/3r^|^(ri<.,.,rA012 

/d3r! 


(5) 


rwhere  pi(r)  is  the  one-body  density  of  the  system,  as  well  as  other  diagrammaticallv  defined  two-bodv 
and  three-body  functions. 

The  chief  remaining  task  in  the  formal  theory  is  to  find  a  tractable  relationship  between  the  ingre¬ 
dients  of  the  wave  function  (4>,  «j( r),  and  ti2(rMr;))  and  the  derived  quantities  such  as  the  density,  pair 
distribution  function,  exchange/correlation  energy,  etc.  For  highly  correlated  systems,  this  is  achieved 
using  the  fermion  version  of  the  hypernetted  chain  resummation  (FIINC),  which  then  requires  an  ap¬ 
proximation  for  manageable  calculations.  In  this  work  we  use  the  minimal  acceptable  approximation, 
FIINC//0,  whose  chief  ingredient  is  the  “direct”  two-point  distribution  function  <7da(r,r'),  which  is  one 
plus  the  sum  of  all  direct  two-point  Born-Mayer  type  diagrams  where  there  is  no  exchange  between 
the  two  external  points  (r,r')  with  any  internal  points.  Thus  0<^(r,r')  is  a  renormalized  version  of 
exp(t/2(r,r')),  and  r^(r,r')  =  <7d<*(r,r')  -  1  can  be  viewed  as  the  dimensionless  correlation  hole  around 


a  particle  at  r'  as  a  function  of  r.  In  the  FHNC//0  approximation,  the  pair  distribution  function  is  given 
by 

<Kr>r')  =  [1  +  rrfrf(r.r,))|^r(r,r')  +  ' Tp  (r)pi(y)  ^ dd ]  (r>r/)  j>  (6) 

where  Sf( r,  r')  and  <fr(r,r')  are  the  structure  factor  and  pair  distribution  function  of  the  uncorrelated 
state  [A  *  jB](r,r')  means  the  convolution  integral  of  the  two  two-body  functions,  and  the  tilda 
above  a  function  is  defined  by  /l(r,r')  =  ,  r') y/pTC^7)- 

RESULTS  FOR  Be  AND  Ne 

The  energy  of  atomic  systems  is  almost  entirely  accounted  for  by  Hartree  Fock.  Correlations  in 
electronic  systems  are  significant  only  for  low  density  regimes,  which  means  the  outer  part  of  the  atomic 
systems.  The  correlation  energy  Ec  is  defined  as  the  difference  between  the  total  energy  Etot  and  the 
Hartree  Fock  energy  Ejjf' 

Ec  =  Etot  -  Ehf  (7) 

This  correlation  energy  consists  of  two  parts,  a  positive  contribution  which  comes  from  the  fact  that 
the  single  particle  orbitals  are  solutions  of  the  correlated  Ilartree-Fock  equation,  and  thus  are  not 
the  best  nncorrclaied  Ilartree-Fock  orbitals,  and  the  remaining  contributions  which  arise  primarily  from 
the  difference  between  the  particle-hole  interaction  and  the  bare  Coulomb  interaccion,  as  well  as  the 
differences  between  the  fully  correlated  and  the  uncorrelated  <?r(r,r').  This  latter  negative 

quantity  exceeds  the  positive  Hartree  Fock  shift,  giving  a  total  negative  value  for  the  correlation  energy. 

This  correlation  energy  is  of  the  order  of  1%  or  less  of  the  total  energy.  Nevertheless,  Clementi  and 
Veillard3  have  obtained  this  energy  from  experiments  on  four  and  ten  electron  atomic  systems.  Their 
results,  shown  in  the  Table,  have  been  corrected  for  center  of  mass  energy  and  relativistic  effects,  and 
thus  can  be  compared  directly  to  our  results  for  non-relativistic  atoms  with  an  infinite  mass  nucleus.  It 
is  seen  from  the  table  that  our  results  are  in  good  agreement  with  these  experiments;  some  differences 
should  be  expected  from  the  FHNC//0  apprximation,  which  produces  approximately  10%  errors  in  the 
case  of  jellium. 

Table:  Correlation  energies  for  Neon  and  Beryllium  atoms  in  atomic  units,  a:  ex¬ 
perimental  data3;  b:  present  work;  c:  LSD  results4;  d:  generalized  gradient 
expansion  of  LM5;  and  e:  generalized  gradient  expansion  of  Perdew6. 


a 

b 

c 

d 

e 

Be 

-0.0944 

-0.096 

-0.224 

-0.099 

-0.094 

Ne 

-0.39 

-0.33 

-0.74 

-0.41 

-0.39 

The  density  profile  of  these  atoms  is  very  close  to  the  Hartree  Fock  densities,  although  the  small 
differences  which  appear  in  the  low  density  tail  of  the  atoms  is  a  significant  effect  in  the  correlation 
energies. 


Fig .  lc.  Same  as  Fig.  la  for  a 
parties  located  at  a  distance  2  a. 
u.  from  the  nucleus. 


Fig .  2a .  Same  as  Fig.  lb  for 
the  Ilartree-Fock  approximation 
9F(  r,r'). 


v 


Fig.  2b.  Same  as  Fig.  lb 
for  the  correlation  contribution 
9dd(*>  r')  =  1  +  r^(r,r/). 

The  pair  distribution  function  is  a  familiar  tool  for  examining  two-body  correlations  in  classical  and 
quantum  fluid  systems,  and  gives  a  similar  but  more  complex  view  of  these  correlations  in  inhomogeneous 
systems.  The  spherical  symmetry  of  closed  shell  atomic  systems  simplifies  the  exhibition  of  this  structure 
somewhat.  Thus  <7(r,r')  depends  only  on  the  distances  r  and  r9  from  the  origin  and  the  angle  between 
the  two  vectors.  In  Figure  1  we  present  this  function  as  a  function  of  r  for  fixed  r',  where  the  origin  (the 
nucleus)  and  r'  define  the  x  axis,  and  ^(r,rx)  is  plotted  above  the  plane  defined  by  r  and  r (For  the 
spherical  atom,  all  such  planes  are  equivalent).  It  can  be  seen  in  Figure  1  that  the  exchange-correlation 
hole  (the  depression  of  the  value  of  <?(r,r')  as  r  approaches  r')  becomes  increasingly  important  as  r 9 
increases  into  the  low  density  tail  of  the  atom.  Of  course  even  the  Hartree  Fock  result  must  contain 
an  exchange  hole  by  the  Pauli  principle,  so  the  question  of  how  much  of  the  exhange-correlation  hole 
can  be  attributed  to  exchange  ^  e.,  statistical  correlations)  and  how  much  to  correlation  (i.e.,  dynamical 
correlations)  requires  further  examination.  Fortunately,  as  Eq.  6  shows,  these  effects  factor  tc  a  very 
good  approximation.  Thus  #(r,r')  is  approximately  the  product  of  <7<w(r,r')  and  <7F(r»r0>  where  the 
former  is  dominated  b\  dynamic  correlations  and  is  1  if  the  t\\o-bod\  interaction  vanishes,  and  the  latter 


is  the  free  particle  pair  distribution  function,  therefore  containing  only  statistical  correlations.  These 
two  quantities  and  g  are  each  shown  in  a  particular  case  in  Figure  2,  where  it  is  seen  that  these  two 
factors  each  make  their  own  distinct  contribution  to  the  exchange-correlation  hole  as  well  as  the  rest  of 
the  two-  body  structure. 

One  version  of  the  local  density  approximation  is  obtained  by  approximating  <7(r,r')  by 
gj(\r  -  r'|,  {/?)r(r*),  where  0j(r,/>)  is  the  pair  distribution  function  for  jellium  at  density  {p)r,r'i  which 
is  some  representative  density  in  the  vicinity  of  the  two  points  r,!*7,  e.g.,  the  geometric  or  arithmetic 
average  of  p(r)  and  p(r').  It  is  clear  from  the  figures  that  these  approximations  are  not  appropriate  for 
<j(r,r');  indeed,  it  also  fails  for  gr(r}rf). 

DISCUSSION 

Small  atoms  provide  one  of  the  most  stringent  tests  of  a  theory  of  electron  correlations,  since  the 
density  falls  ven  rapidly  to  zero  away  from  the  nucleus.  Consequently  our  results  for  neon  and  beryllium 
are  very  encouraging.  Preliminary  results  for  other  closed  shell  four  and  ten  electron  systems  (Ne+6, 
Mp+2,  and  Cn+1°)  are  similarly  encouraging.  In  principle,  the  theory  can  also  be  applied  to  open  shell 
systems;  however  the  lack  of  spherical  symmetry  will  require  additional  efiort. 

It  is  abundantly  clear  in  these  simple  systems  that  the  local  density  approximation  is  inadequate 
for  the  correlation  energies.  This  was  expected  by  those  who  introduced  this  approximation7  and  was 
demonstrated  by  comparisons  between  LDA  theories4  and  experiments  for  the  four  and  ten  electron 
closed  shell  atoms  as  is  seen  in  the  Table.  It  is  also  seen  there  that  adjusted  gradient  expansions5,6  give 
agreement  with  experiment  which  is  comparable  to  our  theory.  However  it  was  necessary  in  each  of 
these  gradient  expansions  to  fit  the  data  either  at  one  atom  (Perdew6  adjusts  his  theory  to  fit  Are)  or 
to  several  systems  (Langreth  and  Mehl5,  who  fit  Be}  Ne+6,  and  Are);  an  unadjusted  gradient  expansion 
actually  gives  worse  agreement  with  experiment  than  the  LDA.  Thus  we  conclude  that  the  Jastrow 
Euler-Lagrange  theory  is  a  very  promising  method  for  including  none-local  density  effects  in  atomic  and 
molecular  systems. 

Since  this  theory  is  also  a  theory  of  the  pair  distribution  function,  it  provides  much  information 
about  the  two-body  structure  of  the  system,  and  permits  a  fairly  clean  separation  between  dynamic 
correlations  and  statistical  correlations. 
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INTRODUCTION 

One  of  the  most  exciting  prospects  on  the  current  nuclear  scene  is  the  promise  that 
precision  high-energy  electron  scattering  experiments  will  reveal  new  (and  perhaps  unfore¬ 
told)  aspects  of  nuclear  structure  and  dynamics.  The  search  is  on  for  distinctive  signatures 
of  subnucleonic  degrees  of  freedom,  and  especially  fo  manifestations  of  the  underlying 
quarkic  substructure  of  nuclei.  However,  to  reach  any  definite  conclusions  regarding  such 
effects,  it  is  necessary  that  we  know,  with  precision,  the  values  which  are  predicted  for  the 
measured  quantities  by  the  conventional  picture  of  nuclei.  In  the  conventional  picture,  a 
nucleus  is  composed  of  nucleons  alone,  moving  nonrelativistically.  The  nucleonic  consti¬ 
tuents  are  considered  to  interact  via  bare  potentials  which  reproduce  the  few-nucleon  data 
while  obeying  certain  constraints  imposed  by  fundamental  symmetries  and  by  meson- 
exchange  theory.  Even  at  this  rather  superficial  level,  one  is  confronted  with  a  very 
difficult  many-body  problem,  essentially  nonperturbative  because  of  the  strong  short-range 
interactions  among  the  nucleons.  It  should  therefore  be  no  surprise  that  mean-field  theory 
(in  old  language,  the  shell  model)  fails  in  experimental  settings  where  large  momentum 
transfers  take  place  and  the  high-momentum  components  of  the  nuclear  wave  function  are 
being  probed.  A  proper  understanding  of  the  quantitative  implications  of  the  conventional 
nuclear  picture  requires  a  careful  and  coherent  treatment  of  the  fine-scale  spatial  correla¬ 
tion  structure  and  the  collective  properties  which  emerge  from  solution  -  or  approximate 
solution  -  of  the  many-nucleon  Schrodinger  equation  of  the  nucleus.  Only  when  this 
refined  microscopic  description  is  achieved,  and  the  conventional  picture  still  found  want¬ 
ing,  can  a  legitimate  claim  be  made  that  the  relevant  experiments  have  penetrated  v  a 
more  fundamental  level  of  nuclear  hadronic  dynamics. 

One  focal  point  of  the  ongoing  refinement  of  conventional  nuclear  theory  is  the  longi¬ 
tudinal  response  function  Ri(q, CO)  of  heavy  nuclei  in  the  quasielastic  regime.  In  treating 
this  property,  mean-field  theory  has  proven  inadequate  and  it  is  necessary  (but  perhaps  not 
sufficient)  to  incorporate  the  dynamical  correlations  among  nucleons  in  a  consistent  and 
accurate  manner.  A  microscopic  understanding  of  the  dynamic  structure  function  S(q, oo) 
of  the  hypothetical  problem  of  infinite,  symmetrical  nuclear  matter  should  yield  valuable 


insights  into  the  effects  of  correlations  in  the  observed  longitudinal  response.  In  particular, 
Fantoni  and  Pandharipande1  have  argued  for  an  approximate  proportionality  of  Ri(q,(£>) 
and  5(<7,co)  at  large  q,  the  proportionality  factor  being  given  by  the  absolute  square  of  the 
proton  form  factor  (see  also  Ref.  2). 

The  work  of  Fantoni  and  Pandharipande1  and  Fabrocini  and  Fantoni,2  based  on  the 
method  of  correlated  basis  functions  (CBF),3,4  represents  a  substantial  advance  over  all 
previous  efforts  on  this  problem,  providing  state-of-the-art  evaluation  of  the  longitudinal 
response  function  of  nuclear  matter  (and  more  directly,  its  dynamic  structure  function 
S(q,  (£>)),  with  a  realistic  two-nucleon  interaction  and  simulated  three-body  potential  as 
input.  Within  the  CBF  scheme,  the  calculation  was  performed,  roughly  speaking,  at  the 
Tamm-Dancoff  level.  The  correlated  random-phase  approximation5,6,4  (CRPA),  which 
performs  the  ring  sums  within  CBF  perturbation  theory,  opens  the  way  for  considerable 
improvement  of  some  aspects  of  this  calculation,  particularly  at  the  lower  momentum 
transfers  where  the  effects  of  the  particle-hole  force  are  most  evident.  In  this  paper  we 
present  an  initial  application  of  CRPA  theory  to  infinite  nuclear  matter.  The  goals  of  the 
calculation  are  modest:  Preliminary  to  a  full  calculation  with  a  realistic  nucleon-nucleon 
interaction  with  elaborate  state-dependence,  we  shall  explore  the  predictions  of  CRPA  for 
a  simplified  model  of  the  nucleon-nucleon  interaction  which  has  frequently  been  used  as  a 
test  case,  namely  the  V2  potential.7  Moreover,  we  shall  apply  a  simplified,  local  version  of 
CRPA  which  has  proven  successful  in  applications  to  spin-polarized  liquid  3 He  and  the 
electron  gas.  It  is  found  that  this  approach  already  leads  to  results  of  qualitative  or  semi- 
quantitative  value,  which  may  be  systematically  improved  as  the  techniques  of  CRPA,  and 
more  generally  the  CBF  theory  of  dynamical  response,  are  further  developed. 

In  a  broader  context,  microscopic  evaluation  of  the  density-density  response  function 
n(<7,co)  of  nuclear  matter  (whose  imaginary  part  gives  S(q,  co)),  together  with  consistent 
evaluation  of  the  self-energy  £(/:,£),  yield  fundamental  information  about  the  elementary 
excitations  of  this  hypothetical  hadronic  system.  The  properties  of  collective  modes, 
typified  by  the  zero-sound  dispersion  relation,  may  be  extracted  from  r%,  co),  while  the 
nature  of  single-particle  excitations  is  revealed  by  £(£,£),  from  which  one  may  derive  an 
energy-dependent  effective  mass.  These  properties  have  obvious  importance  for  a  deeper 
understanding  of  nuclei.  They  are  likewise  basic  to  a  description  of  the  structure,  dynam¬ 
ics,  and  thermal  history  of  neutron  stars  -  being  essential  to  the  evaluation  of  such  quanti¬ 
ties  as  the  specific  heat,  viscosity,  superfluid  gap,  etc.  Since  empirical  constraints  on  the 
properties  of  neutron-star  material  are  limited  in  the  extreme,  such  astrophysical  applica¬ 
tions  make  it  doubly  important  to  hone  and  test  our  many-body  calculational  methods. 

LOCAL  CORRELATED  RANDOM-PHASE  APPR  'XIMATION 

The  correlated  random-phase  approximation  (CRPA)  extends  ordinary  RPA  to 
strongly-interacting  systems  like  liquid  3  He,  nuclear  matter,  and  nuclei.  Since  this 
approach  has  been  developed  in  detail  in  other  places,5,6,4,8,9  we  shall  only  describe  its 
main  features  and  display  the  working  formulas  of  a  local  approximation  proposed  by 
Krotscheck.6  Here,  “local”  implies  that  the  panicle-hole  force  is  taken  to  depend  only  on 
the  momentum  transfer  Tiq  =  |  p  -  h  |  in  the  direct  panicle-hole  channel  (apart  from 
momentum-conserving  delta  functions).  We  confine  our  attention  to  the  uniform,  infinite 
Fermi  medium,  for  which  the  appropriate  model  states  are  Slater  determinants  of  plane 
waves,  and  complete  the  specification  of  the  underlying  correlated  basis  by  the  adoption  of 
a  Jastrow  correlating  factor  F  =  IT/  <;/  (r,y). 

In  terms  of  an  irreducible  particle-hole  (ph)  interaction  U  which  CBF  theory  gen¬ 
erates  for  a  given  bare  (and  possibly  singular)  two-body  potential,  correlated  RPA  for  our 
problem  looks  exactly  like  the  familiar  RPA  for  a  weak  potential,  aside  from  a  minimal 
energy  dependence  of  U  which  can  usually  be  ignored.  Implementation  of  the  full  CRPA 


approach  including  exchange  nonlocalities  is  in  principle  straightforward,  but  in  practice 
requires  considerable  numerical  effort,  entailing  (for  example)  discretization  of  the  matrix 
eigenvalue  equations  on  a  grid  in  momentum  space.8  To  provide  a  simpler  alternative, 
Krotscheck6  has  constructed  a  local,  energy-independent  approximation  to  U  which  is 
designed  to  preserve  certain  fundamental  relations  of  the  ingredients  of  U  to  the  static 
structure  function  S(q)  of  the  Jastrow  ground-state  trial  function  and  to  its  graphical 
derivative  S'(q).  One  of  these  relations  is  the  optimization  condition  on  the  Jastrow  pair 
correlation  function  /  (r),  which  reads 


5<H  >0 
5ln  f2(r) 


=  A(r)  =  0  . 


(1) 


In  this  expression,  <H  >0  is  the  energy  expectation  value  in  the  Jastrow  trial  ground  state 
II  f(rij)<$0,  where  d>0  represents  the  noninteracting  Fermi  sea.  The  proposed  local  approx¬ 
imation  to  the  particle-hole  interaction  is  simply 

U(q)  =  mS~2(q)+  >  (2) 


where  Sp (q)  is  the  static  structure  function  of  the  noninteracting  Fermi  gas  and 
A (q)  =  pj  d3r  exp(/q-r)  A(r)  =  (q)~  1]  +S'(q)  . 


The  Jastrow  S(q)  entering  (2)-(3)  may  be  evaluated  with  good  accuracy  by  solving  the 
(nonlinear)  FHNC/C  equations,  while  its  graphical  derivative  S'(q),  appearing  in  (3),  may 
be  obtained  by  solving  the  (linear)  FHNC/C'  equations.10  The  vanishing  of  A  is  equivalent 
to  the  optimization  condition.  Hence  for  optimal  Jastrow  correlations  the  particle-hole 
force  U {q)  depends  only  on  S  ( q )  and  properties  of  the  noninteracting  system.  For  optimal 
correlations,  the  choice  (2)  for  U  (q)  is  just  what  is  needed  to  regain  S (q)  from  the 
density-density  response  function  through  the  fluctuation-dissipation  theorem,  if  the  col¬ 
lective  approximation  n£  =  (fi2q2/m)[7i2(£>2-(1i2q2/2rnSp(q))2]~l  is  used  for  the 
particle-hole  propagator.6,11 


Having  adopted  a  local  particle-hole  force,  one  has  quite  standard  algebraic  RPA  for¬ 
mulas  which  lead  to  the  familiar  RPA  expressions  for  the  density-density  response  func¬ 
tion  and  self-energy12,13: 


n(4,  co) 


no(?,co) 

l-U(q)n0(q,a>) 


(4) 


T(k,Q)  =  U(0)  +  Upxk(k)  +  -^r-jd2qdClGo(k-q,Cl-^)U2(q)U(q,Q)  .  (5) 

(27t)4p 

In  (4),  the  response  function  Tl0(q,  co)  is  the  particle-hole  propagator  of  the  free  Fermi  gas, 
i.e.,  the  Lindhard  function.  In  (5),  Ga  is  the  free  single-panicle  Green’s  function,  U0  is  a 
constant  (related  to  the  chemical  potential),  and  Up ock(^)  is  lhe  term  of  the  panicle- 
hole  force  U(q). 

The  dynamic  structure  function  and  the  properties  of  zero  sound  (if  present)  are 
derived  from  the  relation  (4)  in  the  usual  manner.  Trivially, 

S(q,  co)  =  Im  nfo,  co)  .  (6) 

7Z 

The  zero-sound  dispersion  relation  co  =  co zs(q)  is  determined  by  the  roots  of  the  denomina¬ 
tor  of  (4),  i.e.  the  roots  of 


l-U(q)  Re  Tl0(q,  co)  =  0 


(7) 


in  the  region  where  Im  Tl0{q,  co)  =  0,  The  strength  Z1S  of  the  zero-sound  mode  is  given  by 


Z;!(q)  =  U\q) 


d 

da> 


Reno(<7,co) 


(8) 


Note  that  we  are  considering  only  the  density  channel;  more  generally  a  spin-isospin 
decomposition  can  be  made,  and  collective  modes  other  than  zero  sound  can  be  studied. 

A  version  of  the  single-particle  energy  spectrum  can  be  obtained  from  the  self-energy 
(5)  using  the  on-shell  prescription 

e*  =  tk  +  E(fc,  Cl) ,  Q.  =  tk  =  T?k2 Hm  .  (9) 

The  on-shell  effective  mass  m*  is  then  given  by 

(m*Tl  =k~l  dek/dk  .  (10) 


A  qualitative  shortcoming  of  the  local  correlated  RPA  (LCRPA),  seen  in  the  form  (4) 
of  the  polarization  propagator,  is  that  the  effects  of  dynamical  correlations  on  the  continu¬ 
ous  portion  of  S(q,(Ai),  in  regions  of  the  07,co)  plane  where  no(q,  co)  vanishes,  are  not 
accessible.  It  is  in  fact  just  such  effects  which  were  examined  many  years  ago  by  Czyz 
and  Gottfried.14  Thus  any  useful  comparison  with  their  work  is  precluded.  Although  the 
(<7,G))  domain  corresponding  to  individual  \plh  excitations  is  the  same  in  LCRPA  as  in 
the  free  system,  the  RPA  denominator  in  expression  (4)  introduces  nontrivial  correlation 
effects  in  that  region;  moreover,  outside  that  region,  for  lmno((?,0))  =  0,  zero  sound  may 
emerge  as  a  distinct  collective  mode,  corresponding  to  vanishing  of  the  denominator. 

Local  CRPA  will  also  suffer,  at  a  quantitative  level,  from  the  static  nature  of  the 
effective  interaction  U  ( q )  appearing  in  (4).  By  contrast,  dynamic  screening  is  known  to 
be  important  in  the  electron  gas  at  metallic  densities.15  Moreover,  one  does  not  expect  the 
momentum  dependence  of  the  self-energy  to  be  faithfully  predicted  within  the  LCRPA 
scheme,  especially  in  the  very  delicate  example  of  unpolarized  liquid  3  He  (Ref.  11).  In 
spite  of  its  shortcomings,  LCRPA  offers  a  simple  and  straightforward  microscopic  touch¬ 
stone  for  phenomenological  theories  of  comparable  structure,  such  as  the  polarization- 
potential  model. 

In  the  polarization-potential  approach  of  Aldrich  and  Pines,16  which  has  been  adapted 
to  nuclear  problems  by  Pines,  Quader,  and  Wambach,17  the  density-density  response  func¬ 
tion  n(<7,co)  is  expressed  in  a  form  similar  to  (4).  However,  the  Lindhard  function  appear¬ 
ing  in  (4)  is  replaced  by  a  more  complicated  propagator  accounting  both  for  a  single-pair 
effective  mass  different  from  the  bare  value  and  for  the  presence  of  multipair  excitations. 
Secondly,  the  static,  local  particle-hole  interaction  fls{q)  of  this  approach  is  supplemented 
by  a  wave-vector  and  frequency-dependent  contribution  (c o2/q2)f.l(q)  corresponding  to 
backflow.  Beyond  these  structural  differences  there  is  the  important  conceptual  distinction 
that  our  panicle-hole  interaction  is  determined  microscopically,  whereas  the  polarization 
potentials  of  the  Aldrich-Pines  theory  are  determined  by  a  combination  of  sum  niles  and 
phenomenology. 


RESULTS  FOR  v?  HOMEWORK  MODEL 

Based  on  the  local  correlated  RPA  scheme,  we  have  have  carried  out  a  numerical 
study  of  the  dynamical  response  of  a  model  of  nuclear  matter  which  is  extremely  simple, 
yet  may  capture  important  aspects  of  its  correlation  structure.  The  bare  interaction 
between  the  nucleons  is  taken  as  the  V2  “homework-model”  potential.7  This  two-nucleon 
potential  consists  of  the  central  part  of  the  3Si  -3Di  component  of  the  Reid  soft-core 
interaction,  assumed  to  act  in  all  partial  waves.  It  has  seen  wide  use  in  tests  of  many-body 
methods  (see,  for  example,  Refs.  18,19).  Specifically, 


v2(r)  =  [9924.3exp(-4.2r)  -  3187.8exp(-2.8r) 

+  105.468exp(-1.4r)-  10.463exp(-0.7r)]/(0.7r)  .  (11) 

Our  calculations  are  based,  primarily,  on  the  parameterized  Jastrow  correlation  factor 

f  (r)  =  cxp[-Ac~Br(l-e~r/D)/r]  ,  (12) 

with  parameters  which  minimize  the  Jastrow  ground-state  energy  expectation  value  Ej  at 
single-particle  level  degeneracy  v  =  4  and  given  density,  Ej  being  computed  by  the 
Metropolis  Monte  Carlo  algorithm.20  Additionally,  we  have  examined  the  particle-hole 
force  generated  by  two  versions  of  the  correlation  function 

/  (r)  =  (1  -e_r2/fc2)/*  +  grmt~r1'^  (13) 

studied  by  Benhar  et  al 21  In  the  simpler  version  (Bl),  the  parameter  g,  which  measures 
the  overshoot  of  /  (r)  above  unity,  is  set  zero  and  the  remaining  parameters  b  and  n  fixed 
by  minimization  of  the  energy  expectation  value  truncated  at  three-body  cluster  order.  In 
the  other  version  (B2),  the  three-body  cluster  approximation  to  the  energy  is  minimized 
(approximately)  with  respect  to  the  five  parameters  b,  g,  y,  m  and  n,  subject  to  the  sequen¬ 
tial  (or  normalization)  condition22  as  constraint. 

The  three  correlation  functions  are  labeled  (in  order)  C,  Bl,  and  B2.  None  of  these 
choices  is  optimal  in  the  sense  of  being  a  solution  of  the  Euler  equation  (1).  The  sequential 
condition  imposed  in  the  case  of  B2  is  probably  too  stringent,  while  B 1  is  presumably  too 
crude  a  choice,  although  both  alternatives  have  the  virtue  of  promoting  good  convergence 
of  the  cluster  expansion  of  the  Jastrow  energy  expectation  value.  The  function  C,  being 
determined  through  a  Monte  Carlo  evaluation,  may  provide  a  decent  representation  of  the 
true  optimal  function,  except  in  the  large- r  region  which  contributes  little  to  the  ground- 
state  energy.  The  deviation  at  large  r  will,  however,  be  reflected  in  the  behavior  of  U(q)  at 
small  q.  In  particular,  it  will  produce  a  significant  departure  of  A (q)  from  zero  at  low  q 
values.  In  the  present  work,  as  in  Ref.  9,  we  have  chosen  to  set  A(^)  equal  to  zero,  since 
the  behavior  of  the  U  ( q )  at  small  q  is  already  suspect  due  to  the  local  approximation  itself. 
This  conforms  with  a  strategy  suggested  by  Krctscheck,6  who,  in  considering  the  Bijl- 
Feynman  dispersion  relation,  found  that  it  is  a  better  approximation  to  assume  that  the 
correlation  function  has  been  optimized  -  even  if  the  actual  optimization  has  not  been  car¬ 
ried  through  -  than  to  include  the  A (q)  correction. 

Numerical  results  with  the  C  correlations  have  been  obtained  for  symmetrical  nuclear 
matter,  i.e.,  level  degeneracy  v  =  4,  at  kp  =  1.39  fm-1,  corresponding  to  a  density 
p  =  0.182  fm-3  near  nuclear  saturation;  and  for  pure  neutron  matter,  i.e.,  v  =  2,  at 
kp  =  1.75,  2.25,  2.90  fm”1,  corresponding  to  three  densities  p  =  0.182,  0.386,  0.822  fm-3 
of  relevance  to  the  study  of  neutron  stars.  For  both  level  degeneracies  we  have  examined  a 
range  of  wave-number  transfers  q  from  0  to  about  4  fm-1 .  Results  with  Bl  and  B2  corre¬ 
lations  are  available  only  for  symmetrical  nuclear  matter,  at  Fermi  wave  number 
kp  =  1.4  fm-1. 

For  reference,  we  have  also  generated  nuclear-matter  (v  =  4)  results  for  the  case  that 
the  two  nucleon  interaction  is  replaced  by  a  pure  hard-core  potential  of  core  radius  c  -  0.5 
fm.  As  in  an  earlier  calculation,9  where  c  -  0.4  fm  was  studied,  the  simple  form 

f(r)  =  0,  r<0  , 

=  1 -exp[-|i(r-c)]  ,  r>c  (14) 

is  assumed  for  the  Jastrow  pair  correlation  function.  An  optimal  determination  of  the 
parameter  (I  has  been  carried  out  by  Flynn23  within  the  FHNC/C  approximation,  for  both 
core  sizes  ( c  =  0.4  and  0.5  fm). 


We  confine  our  discussion  to  a  selection  of  results  for  symmetrical  nuclear  matter. 
(A  more  detailed  report  will  be  published  elsewhere.)  The  CBF  particle-hole  interaction 
appropriate  to  the  v2  potential  at  kp  =  1.39  fm**1,  for  the  C  choice  of  /  (r),  is  shown  in 
Fig.  1,  along  with  the  c  =  0.5  fm  hard-sphere  result  at  the  same  density.  Due  to  the  differ¬ 
ence  between  the  respective  correlation  functions  (the  v2  correlations  being  slightly 
stronger  beyond  1  fm),  the  panicle-hole  force  for  v2  is  actually  somewhat  more  repulsive 
at  low-to-medium  q's  than  that  for  the  hard-core  potential.  Both  potentials  support  a  col¬ 
lective  mode  corresponding  to  zero  sound.  In  both  cases,  this  mode  is  found  to  emerge 
from  the  particle-hole  continuum  around  0.3-0.4  fm-1  and  sink  back  into  it  at  about  1.53 
fm-1.  The  zero-sound  dispersion  relation  oo7S(q)  and  the  strength  Z^iq)  of  the  zero-sound 
mode  are  plotted  in  Fig.  2,  for  the  v2  potential.  The  corresponding  curves  for  the  hard¬ 
core  potential  are  nearly  coincident  with  those  for  v2. 


Fig.  1.  Wave-number  dependence  of  local  particle-hole  interaction  U(q )  for  v2  model  of  nuclear 
matter,  based  on  correlation  functions  C,  Bl,  and  B2  (as  labeled);  and  for  hard-sphere  nucleons 
with  c  =0.4  and  0.5  fm  (as  labeled),  based  on  optimized  correlation  function  (14).  Density  near 
0.18  nucleons/fm3. 


Fig.  2.  Zero  sound 
dispersion  relation  (£>7Ji(q) 
(solid  curve)  and  wave 
number  dependence  of 
zero  sound  strength 
ZK{q)  (dashed),  for  v2 
model  of  nuclear  matter 
at  kp  =  1.39  fm-1 ,  based 
on  C  correlation  func¬ 
tion. 


q  (fm"1) 


It  may  seem  surprising  that  the  effective  interaction  derived  from  the  soft-core  v2 
potential  is  so  close  to  that  for  the  hard-core  potential.  The  concurrence  implies  that  the 
two  potentials  are  similarly  effective  in  lifting  particles  from  below  to  above  the  Fermi 
surface,  at  least  when  the  momentum  transferred  to  the  particle-hole  pair  is  not  excessive 
and  the  density  is  near  the  saturation  value  for  nuclear  matter.  In  turn  this  observation  sug¬ 
gests  that  a  useful  measure  of  the  strength  of  the  particle-hole  force  may  be  furnished  by 
the  wound  parameter  of  Brueckner  theory,  which  may  be  interpreted  as  the  average  deple¬ 
tion  of  the  Fermi  sea  produced  by  the  given  bare  interaction.  Indeed,  the  wv.und  parame¬ 
ters  for  the  two  potentials  (v2  and  c  =0.5  fm  hard  core)  are  quite  similar.  As  is  well 
known,18  the  Jastrow  variational  analog  of  the  Brueckner  wound  parameter  is 
k  =  pjd3r  | f  (r)~  1]2[1  -  v~ll2(kpr)),  where  l(x)  =  3.x-3(sin;t-;ccos;t).  For  the  C  choice 
of  the  correlation  function,  we  have  k  =  0.232,  which  is  to  be  compared  with  the  value 
k  =  0.263  for  correlation  function  (14)  with  p.  =  =  2.5  fm-1 . 

We  may  cany  these  considerations  somewhat  further  by  examining  the  curves  in  Fig. 
1  labeled  B1  and  B2.  The  particle-hole  force  U  ( q )  for  B1  or  B2  is  seei  to  differ  substan¬ 
tially  from  that  corresponding  to  conrelation  function  C,  being  much  weaker  at  smaller  q 
values  (below  2  fm-1 )  and  somewhat  stronger  at  large  q.  The  large  differences  cannot  be 
due  to  the  slight  discrepancy  in  densities  (0.185  fm-1  (B1  and  B2)  compared  to  0.182  fm-1 
(C)).  Instead,  it  must  be  a  reflection  of  the  reduction  in  the  flexibility  of  /  (r)  entailed  by  a 
procedure  employing  a  truncated  cluster  expansion,  the  range  of  variation  being  restricted 
to  functions  for  which  the  cluster  expansion  of  the  energy  is  rapidly  convergent.  The  sub¬ 
stantial  difference  in  predictions  for  U{q)  is  (as  expected)  accompanied  by  substantial 
differences  in  the  wound  parameters  for  B 1  and  B2,  from  those  for  the  correlation  function 
C.  The  specific  values  are  k  =  0.1 46  (Bl)  and  k  =  0.143  (B2).  Our  thesis  that,  the  wound 
parameter  is  a  major  determinant  of  the  overall  behavior  of  U (q),  gains  additional  support 
from  the  observation  that  correlation  functions  Bl  and  B2,  with  approximately  the  same 
wound  parameters,  yield  U  (q)  results  which  are  quantitatively  as  well  as  qualitatively 
similar.  To  complete  the  comparison,  we  have  included  in  Fig.  1  a  U  ( q )  curve  for  a 
c  =  0.4  fm  hard-sphere  interaction  at  kp  =  1.39  fm-1 ,  described  by  means  of  the  correlation 
function  (14)  with  p  =  PoPt  =  3  fm-1.  The  corresponding  wound  parameter  is  k  =  0.139. 
Once  again  it  is  seen  that  correlation  functions  with  similar  wound  parameters  yield  simi¬ 
lar  lesults  for  the  panicle-hole  force  U (q). 

The  discussion  now  focuses  on  the  results  obtained  for  the  C  choice  of  correlations. 
The  dynamic  structure  function  S(q,  co)  of  the  v2  model  of  nuclear  matter,  derived  from 
the  LCRPA  treatment,  is  plotted  as  a  function  of  energy  transfer  7ko  in  Figs.  3-6,  at  fixed  q 
values  of  1.47,  2.03,  2.76,  and  3.87  fm-1,  respectively.  These  choices  for  q  are  close  to 
values  at  which  experimental  data  on  the  longitudinal  response  function  Ri(q,  co)  of 
medium-  and  large-A  nuclei  are  available24-26  and/or  are  (reasonably)  close  to  q  values 
examined  in  the  microscopic  calculation  of  Fantoni  and  Pandharipande.1  In  the  figures  we 
trace,  for  comparison,  the  Fermi-gas  structure  function  Sp(q,  co)  as  well  as  the  LCRPA 
result  for  the  c  =  0.5  hard-core.  It  has  been  argued  by  Fantoni  and  Pandharipande  that  the 
longitudinal  dynamic  structure  function  Si(q,oS)  of  nuclear  matter  should  not  be  very  dif¬ 
ferent  from  S{q,(£>)  over  the  range  q  -  1. 5-2.8  fm-1,  and,  accordingly,  that  the  theoretical 
nuclear-matter  S(q,(£>)  multiplied  by  the  square  of  a  suitably  chosen  proton  form  factor, 
may  be  sensibly  compared  with  the  experimental  response  function  R^(q,(a)  of  heavy 
nuclei  (cf.  also  Refs.  2,27).  At  a  deeper  level,  such  comparisons  exploit  the  generally 
accepted  view  that  meson-exchange  currents,  pion  production,  delta  excitation,  etc.,  should 
not  be  very  important  in  the  longitudinal  response  function,  for  the  q  values  considered 
here.28  In  this  view  one  must  look  elsewhere  for  a  dramatic  breakdown  of  the  conven¬ 
tional  picture  of  the  nucleus  as  a  nonrelativistic  many-body  problem  with  only  nucleonic 
degrees  of  freedom. 


Fig.  3.  Dynamic  structure  function  S(q, (0)  versus 
energy  transfer  Hca  at  fixed  wave  number  transfcr 
q  =  1 .47  fm-1  and  kp  =  1.39  fm-1.  Solid  curve: 


for  v2  model  of  nuclear  matter  based  on  C  correla¬ 
tion  function.  Dot-dashed:  for  hard  sphere 
nucleons  with  c  =  0.5  fm.  Dotted:  for  free 


Fig.  4.  Same  as  Fig.  3,  except 
q  =  2.03  fm_l  and  there  is  no  zero  sound. 


nucleons.  Zero  sound  contribution  in  v2  case  is 


Itg)  (MeV)  l,w  iMeV) 


Fig.  5.  Same  as  Fig.  3,  except  q  =  2.76  fm-1  and  there  Fig-  6.  Same  as  Fig.  3,  except 
is  no  zero  sound.  <?  =  3.87  fm-1  and  there  is  no  zero 

sound.  Vertical  scale  as  in  Fig.  5. 


As  seen  in  Figs.  3-6,  the  LCRPA  calculation  shows  the  well-known  quenching  of  the 
response  at  low  energies,  compared  to  the  independent-particle-model  result.28-30  This 
effect  is  perhaps  too  severe  in  our  results  at  the  lower  q  values.  The  predictions  for  5  (<7,0) 
are  very  similar  for  V2  and  hard-core  potentials,  suggesting  that  the  response  in  the 
relevant  q  regime  is  insensitive  to  the  choice  of  bare,  state-independent  central  potential, 
provided  the  wound  parameter  is  kept  roughly  the  same.  For  a  given  <7,  the  strength  is 
shifted  to  values  of  Tico  higher  than  the  experimental  peak  energy.  This  effect,  as  well  as 
the  excessive  quenching  at  low  co,  is  presumably  due  to  the  overly  repulsive  character  of 
the  v2  and  hard-core  potentials,  which  act  equally  in  all  partial  waves,  in  contrast  to  the 
strong  partial-wave  dependence  of  realistic  nucleon-nucleon  interactions.  Indeed,  in  the 
cases  q  =  1.47  firf 1  and  q  =  2.03  fm-1  (corresponding  respectively  to  290  MeV/c  and  400 
MeV/c),  there  is  a  very  pronounced  peak  at  the  upper  end  of  the  allowed  range  of  energy 
transfer,  which  is  rather  unphysical.  As  mentioned  earlier,  aside  from  zero  sound,  the 
LCRPA  gives  no  response  outside  the  ( q ,  (D)  region  relevant  for  the  noninteracting  system, 
so  the  strength  which  would  otherwise  be  distributed  to  higher  energies  seems  to  be  “pil¬ 
ing  up”  at  the  high-co  boundary.  The  similar  trends  seen  in  the  results  of  Fantoni  and  Pan- 
dharipande  and  of  Pines,  Quader  and  Wambach  are  much  milder.  As  we  go  to  high  <7,  the 
departures  of  the  LCRPA  result  for  S  (<7,0))  from  the  Fermi  gas  curve,  and  from  the  micros¬ 
copic  prediction  of  Fantoni  and  Pandharipande,  become  less  noticeable.  This  is  an  obvi¬ 
ous  consequence  of  the  decreasing  importance  of  LJ(q)  at  larger  momentum  transfers. 
Since,  quite  generally,  we  may  expect  the  particle-hole  force  to  become  negligible  in  the 
large-<7  regime,  predictions  for  S  (q,(£>)  should  be  relatively  insensitive  to  the  choice  of  the 
bare  interaction,  or  to  the  treatment  within  conventional  nuclear  many-body  theory,  at 
large  enough  q.  Of  course,  excitation  of  subnucleonic  degrees  of  freedom  may  become  a 
significant  factor  in  this  domain. 

For  symmetrical  nuclear  matter  at  kp  =  1.39  fm-1,  a  preliminary  evaluation31  of  the 
LCRPA  self-energy  (5)  using  the  techniques  of  Ref.  32  yields  m*  =  1.1m  for  the  on-shell 
effective  mass  at  the  Fermi  momentum  and  shows  no  significant  enhancement  of  m*(k)  in 
the  vicinity  of  the  Fermi  surface  (cf.  Ref.  33).  If  the  v2  interaction  is  allowed  to  act  only 
in  S  waves  (more  appropriate  in  the  nuclear  context),  a  value  of  m*(kp)  near  0.8m  should 
be  obtained.34,19 

CONCLUSIONS 

The  density-density  response  function  of  symmetrical  nuclear  matter  and  pure  neu¬ 
tron  matter  has  been  studied  under  simplifying  conditions,  leading  to  a  microscopic  theory 
-  the  local  correlated  random-phase  approximation  -  requiring  as  input  only  the  static 
structure  function  corresponding  to  a  Jastrow  description  of  the  ground  state.  This  theory 
has  the  virtues  of  easy  application  and  straightforward  interpretation  of  its  predictions. 
Relative  to  an  ordinary  local  RPA  treatment,  it  takes  account,  approximately,  of  important 
dispersive,  polarization,  and  geometrical  effects  arising  from  the  strong  interactions.  In 
spite  of  its  limitations,  the  theory  has  considerable  value  in  establishing  qualitative  trends. 
For  example,  it  is  found  that  predictions  for  the  dynamic  structure  factor  do  not  depend 
very  strongly  on  the  choice  of  central,  state-independent  potential,  so  long  as  the  wound 
parameter  remains  essentially  unchanged. 

One  serious  limitation  of  the  theory  as  implemented  here  lies  in  the  choice  of  the  v2 
model  interaction  as  the  bare  force.  It  is  well  known  that  state  dependence  of  the  nuclear 
interaction  and  especially  the  presence  .of  a  tensor  component  have  important  conse¬ 
quences  for  most  nuclear  properties.  The  dynamic  response  is  no  exception,  particularly  at 
low  momentum  transfer.35  Nevertheless,  before  proceeding  to  a  realistic  interaction  like 
Urbana  or  Argonne  V14,  and  to  the  introduction  of  three-nucleon  interactions,  the  intrinsic 
limitations  of  the  local  correlated  RPA  approach  should  be  quantified.  As  implemented 


here  for  parameterized  correlation  functions  which  do  not  strictly  obey  the  Euler  optimal¬ 
ity  condition,  one  such  limitation  is  the  neglect  of  the  A (q)  term  in  the  formula  (2)  for  the 
panicle-hole  force.  We  have  argued  that  this  is  not  likely  to  produce  significant  additional 
errors  in  the  description  of  the  panicle-hole  interaction,  since  the  correlation  function 
assumed  in  our  primary  calculations  (of  form  (12))  is  believed  to  be  nearly  optimal.  Devi¬ 
ations  from  optimality  will  mainly  affect  U  ( q )  at  small  q  -  but  in  this  region  the  local 
approximation  must  itself  be  considered  questionable.  These  arguments  notwithstanding, 
the  importance  of  the  A (q)  term,  for  the  correlations  used  here,  is  currently  under  numeri¬ 
cal  investigation,  and  optimal  correlations  will  be  implemented  in  future  calculations. 

The  next  step  is  to  assess  the  effects  of  nonlocalities  in  the  particle-hole  force. 
Within  the  CBF  description  of  elementary  excitations,  this  will  involve  numerical  solution 
of  the  full  CRPA  equations,  with  explicit  inclusion  of  exchange.8  Such  a  treatment  is 
needed  to  remove  the  unphysical  feature  of  local  CRPA  that  (apart  from  possible  collec¬ 
tive  modes)  the  response  is  confined  to  the  ( q ,  co)  domain  in  which  the  dynamic  structure 
function  of  the  noninteracting  system  is  nonzero. 

A  truly  quantitative  description  of  the  response  and  elementary  excitations  in  nuclear 
systems  will  require  a  theory'  which  goes  beyond  correlated  RPA  to  the  explicit  inclusion 
of  correlated  multipair  effects.  Work  in  this  direction  has  begun. 
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INTRODUCTION 

The  two-body  density  matrix  P2(ri,r2,r,1,i*2)  of  the  ground  state  of  a  quantum 
fluid  is  a  rich  repository  of  information  about  its  correlation  structure.  For  exam¬ 
ple,  the  restricted  version  Pi(rj, r2, r^)  =  /^(I’lj^ri)^)  plays  a  crucial  role  in  the 
description  of  final-state  interactions  in  a  novel  theory  of  the  deep-inelastic  neutron 
scattering  from  the  helium  liquids,1  capturing  the  correlation  effects  which  give  rise  to 
deviations  from  the  impulse  approximation  out  to  very  high  momentum  and  energy 
transfers.  This  contribution  will  focus  on  the  microscopic  evaluation  of  p2(ri  i  r2:  ri ) 
for  uniform  Fermi  and  Bose  fluids.  For  simplicity,  we  assume  that  the  ground  state  of 
either  system  is  adequately  described  by  a  Jastrow  trial  function  of  the  appropriate 
symmetry. 

The  formal  structural  analysis  of  the  two-body  density-matrix  elements  is  most 
efficiently  carried  out  in  the  configuration-space  representation,  and  will  be  pursued 
using  well-established  cluster-expansion  procedures,2,3  followed  by  hypernetted-chain 
resummation  of  cluster  diagrams.4,5  On  the  other  hand,  the  physical  interpretation  of 
the  results  is  more  vividly  expressed  in  the  momentum  representation,  i.e..  in  terms 
of  the  generalized  momentum  distribution 

«(P>q)  =  53  <  -laI+q4-qaPa*l*  >  ’ 

jfc 

which  is  related  to  /?2(ri ,  r2,  rj)  by  Fourier  transformation.  (Here,  k  denotes  the 
single-particle  orbital  with  quantum  numbers  k,cr,  where  o  is  the  spin  projection, 
while  k  —  q  =  (k  —  q,  a).)  Using  the  formal  results  obtained  for  the  structure  of 
P2(ri>  r2,  rj)  -  involving  its  closed  expression  in  certain  sums  of  irreducible  cluster 
diagrams  -  we  have  been  able  to  achieve  a  clean  decomposition  of  n(p,q)  into  con¬ 
tributions  from  various  scattering  processes  occurring  in  the  medium.  In  comparing 
Bose  and  Fermi  cases,  certain  features  arising  from  exchange  (Pauli  kinematic  effects) 


will  also  be  evident.  Quantitatively,  the  contributions  from  these  assorted  physical 
effects  are  determined  by  a  set  of  form  factors,  which  are  susceptible  to  evaluation  by 
hypernetted-chain  techniques.  We  shall  present  numerical  results  for  the  form  factors 
depending  on  a  single  momentum  variable,  and  test  the  quality  of  an  estimate  for 
rc(p,q)  used  by  Silver1  as  input  for  his  theory  of  final-state  effects  in  deep-inelastic 
neutron  scattering. 

The  microscopic  treatment  of  n(p,  q)  will  be  developed  for  a  Fermi  system  of 
arbitrary  single-particle  level  degeneracy  v.  The  corresponding  results  for  the  Bose 
case  may  then  be  obtained  by  taking  the  limits  1/-+00  and  fcp  — ►  0+  while  keeping 
the  density  p  =  uk^/Gir2  constant. 

The  role  of  n(p,q)  in  deep-inelastic  neutron  scattering  becomes  more  tangible 
when  we  write  this  quantity  as 

rc(p,q)  =<  $|/v4-qapl'J'  >  -n(p)  ,  (2) 

where  pq  is  the  density  fluctuation  operator  (with  q  ^  0)  and  n(p)  is  the  single¬ 
particle  momentum  distribution  function.  The  expectation  value  (2)  may  be  in¬ 
terpreted  as  a  transition  matrix  element  for  scattering  of  a  particle  out  of  orbital 
P  =  (Pi a>)  into  another  orbital  p  —  q  =  (p  —  q,  a'),  the  process  being  mediated  by  a 
density  fluctuation  of  wave  vector  q.  Thus,  an  evaluation  of  n(p  —  q)  amounts  to  a 
calculation  of  the  rightmost  single-atom  vertex  in  the  Fig.  1. 


Fig.  1.  Deep-inelastic  scattering 
of  neutrons  from  a  (normal)  he¬ 
lium  liquid,  involving  final-state 
interactions  mediated  by  a  den¬ 
sity  fluctuation  of  the  target  sys¬ 
tem  (a  phonon  of  wave  vector  q). 


STRUCTURAL  RESULTS  FOR  p2(ri,r2,r') 

The  ground-state  wave  function  of  the  iV-particle  Fermi  system,  supposed  to 
constitute  a  uniform  fluid,  is  taken  in  the  Jastrow-S later  form  'D  =  C'-1n,-<j/(r,-y)$, 
where  $  is  the  usual  Slater  determinant  of  plane-wave  orbitals,  filling  a  Fermi  sea 
characterized  by  Fermi  wave  number  kp  and  level  degeneracy  is,  while  /(r,-y)  is  a 
Jastrow  pair  correlation  function  to  be  determined  by  minimization  of  the  energy 
expectation  value.  The  constant  C  is  introduced  to  normalize  \f'  to  unity.  The 
results  to  be  described  below  can  be  generalized  in  a  straightforward  manner  to  the 
case  that  the  Jastrow  correlating  factor  is  replaced  by  a  Feenberg  function.6 

The  first  step  of  the  formal  development  is  to  recognize  that  the  generalized 
momentum  distribution  function  (1)  may  be  expressed  as  follows 

n(p,q)  =  Sq0(N  -  1  )n{p)  +  (l  -  Sq0)  <  ^>\N(p,  q)|tf  >  (3) 

in  terms  of  the  one-body  momentum  distribution  n(p)  and  the  expectation  value  of 
a  symmetric  sum  of  two-body  operators, 

N 

N(p,<l)  —  ^T^[eiq  r,op-qip(j)  +  e  q  J<?p_q(p(z)] 

«<i 


(4) 


The  action  of  the  one-body  operator  Op_qip(i)  is  specified  by 

>=  ^pplP  ~  q  >  i  <  ^l°r-q,p  =  fy,p-q  <  P\  ‘  (®) 

Standard  procedures  are  available  for  cluster-expanding  the  expectation  value  of  a 
sum  of  two-body  operators,  the  most  familiar  case  being  the  potential  energy  corre¬ 
sponding  to  a  sum  of  pair  potentials.7  These  procedures  have  been  applied  to  yield 
explicit  results  for  the  cluster  contributions  to  n(p,q)  through  three-body  order,  plus 
selected  four-body  contributions.  The  generalized  Ursell-Mayer  diagrammatic  repre¬ 
sentation  introduced  in  Ref.  2  facilitates  the  bookkeeping.  At  two-body  order  there 
are  7  cluster  diagrams,  at  three-body  there  are  138,  etc.,  the  number  growing  very 
rapidly.  The  analysis  becomes  more  transparent  when  we  go  over  to  the  position 
representation,  inverting  the  Fourier  relation 

n(p,q)  =  J P2(rir2r,1)e-,p^ri~r^e-lq^r,~r^dridr2dr'1  (6) 

to  obtain  the  corresponding  cluster  expansion  for  the  two-body  density-matrix  ele¬ 
ments  /?2(rir2r'i).  There  is  a  very  simple  recipe  for  accomplishing  this  digrammati- 
cally.  One  simply  removes  the  arrows  representing  specific  plane- wave  orbitals  and 
changes  the  field  points  which  they  originally  intersected  into  root  points.  It  then 
becomes  evident  that  most  of  the  diagrams  are  reducible  (factorizable),  i.e.,  they 
consist  of  products  of  simpler  graphs. 

In  the  context  of  extensive  technical  experience  with  cluster  expansions,  the 
diagrams  which  have  been  explicitly  generated  are  sufficient  to  reveal  the  structure 
of  out  to  infinite  cluster  order.  An  exact  representation  of  this  quantity 

is  provided  by  the  expression 

/>2(rir2ri)  =  p2D(rir2r,1)[I(r1ri) +  L(rir2r,1)]  .  (7) 

The  first  factor  collects  the  direct-direct  portions  of  the  full  set  of  diagrams  contribut¬ 
ing  to  P2(i‘ir2r'1).  By  definition,7  direct-direct  diagrams  do  not  have  exchange  lines 
attached  to  any  of  the  root  (or  reference)  points  n,  r2,  and  r[.  The  complementary 
set  of  graphs  contains  only  diagrams  with  exchange  lines  beginning  and/or  ending 
at  two  or  three  reference  points.  Of  these  graphs,  the  ones  with  exchange  lines  at 
two  reference  points  combine  to  form  the  two-point  exchange  factor  jD(rjrj),  while 
the  graphs  with  exchange  lines  at  all  three  reference  points  compose  the  three-point 
function  £(rj  r2r'] ).  Either  of  these  exchange  functions  vanishes  if  any  one  of  the 
coordinates  in  its  argument  recedes  to  an  infinite  distance  from  the  others.  In  the 
Bose  limit  (u  — *  co,  kp  — »  0+,  p  constant),  the  function  £(riri)  approaches  unity 
and  the  three-point  function  )  goes  to  zero.  As  expected,  only  the  direct- 

direct  contribution  p2£>  survives.  Returning  to  (7)  at  finite  u,  this  component  may 
be  compared,  at  a  diagrammatic  level,  with  the  structural  result  which  was  derived 
for  the  Bose-fluid  P2(rir2r,1)  in  Ref.  8.  We  arrive  thereby  at  the  representation 

/>2D(r,r2r[)  =  p/»iD(rir',)/(|r,-r2|)/(|r'l-r2|)exp|-P(r,r2)-/3(r'1r2)-/’(r,r2r'1)]  . 

(8) 

The  generating  functions  P(rir2)  =  -P(|ri  -r2|)  s  P(r)  and  P(rir2r,1)  are  irreducible 
quantities  -  sums  of  irreducible  diagrams  -  just  as  in  the  Bose  case.  To  give  at  least 
a  hint  of  their  character,  Fig.  2  shows  the  leading  cluster  contributions  to  these 
functions.  We  may  remind  the  reader  of  the  relevant  diagrammatic  conventions.2,8 
Open  circles  represent  root  points  (lower  left  circle:  r5,  lower  right  circle:  r\,  upper 


circle:  r2).  Solid  circles  denote  field  points  and  imply  an  integration  p  f.  A  wavy 
[dashed]  line  stands  for  the  correlation  bond  f(rai>)  —  1  [respectively,  f2(rab )  —  1], 
where  r„&  =  |ra  —  r^),  and  r„  and  rj  are  the  coordinate  points  (solid  or  open  circles) 
conned  —1  by  the  bond.  The  diagrams  shown  explicitly  in  Fig.  2  are  in  fact  the 
same  as  the  leading  diagrams  which  occur  for  the  Bose  fluid.  However,  as  in  the 
familiar  examples  of  the  momentum  distribution  function  n(p )  and  the  one-body 
density  matrix  p\  (rj  r't ),  additional,  non-Bose  diagrams  will  arise  from  the  systematic 
introduction  of  exchange  insertions  at  the  field  points  of  the  Bose  diagrams  of  higher 
orders.7’2  These  insertions  are  characterized  by  the  presence  of  exchange  lines  at  one 
or  both  end  points.  An  exchange  line,  drawn  solid  with  an  arrow  attached,  represents 
the  Slater  statistical  bond  l(kpra *),  l(x)  =  3x-3  sin  a;  —  x  cosx. 


-Q(r)  = 

+--- 

-P(r)  = 

cx 

o>  +-" 

Fig.  2.  Leading  cluster  diagrams 
of  the  irreducible  direct  functions 

C?(r),  P(r),  and  P(rir2ri)  and  of 

9 

« 

the  exchange  functions  Li(r)  and 

•P(iW'i)  = 

+" 

L2(r). 

L!(r)  = 

I'D 

-O*- 

L2(r)  =  -  ^  +••• 

The  function  piz?(rirj)  appearing  as  a  factor  in  (8)  is  just  the  direct-direct 
component  of  the  full  Fermi  one-body  density  matrix  pi(riri).  The  structure  of  p\o 
is  well  known  from  previous  work.2-4  This  quantity  is  generated  by  the  irreducible 
phase-phase  correlation  function  <5(ri  ri )  (see  Fig-  2,  and  Figs.  8  and  10  of  Ref.  2), 
according  to 

Pi£>(rir'i)  =  pn0exp[— Q(rirj)]  ,  (9) 

where  n0  =•  expQ(rj  r) )  is  an  overall  strength  factor. 

At  this  point  we  have  identified  three  sums  of  irreducible  diagrams,  namely  Q(r), 
P(r),  and  P(rir2rj)  (where,  and  throughout,  r  =  |rj  —  r2|).  As  in  the  Bose  case 
analyzed  earlier8  (and  in  the  more  familiar  example  of  the  generating  function  U{r ) 
of  the  hypernetted-chain  representation  g(r)  =  /2(r)e-t,^r^  of  the  Bose  radial  dis¬ 
tribution  function7),  these  functions  may  be  decomposed  into  nodal  and  elementary 
components.  Fig.  2  shows  in  fact  the  leading  diagrams  of  nodal  type;  elementary 
diagrams  first  appear  in  the  next  cluster  order.  In  practice,  the  diagram  sums  (Q(r), 
P(r ),  and  P( rir2r'j)  may  be  evaluated  approximately,  to  all  orders,  by  solving  sets 
of  Fermi-hypernetted-chain  (FHNC)  equations,7’4,5  incorporating,  stepwise,  larger 
and  larger  classes  of  elementary  diagrams.  The  simplest  approximation,  denoted 
FHNC/0,  involves  the  neglect  of  all  elementary  components. 

Fermi  exchange  effects  arise  implicitly  from  insertions  at  the  field  points  of  the 


Bose  diagrams  contributing  to  Q(r),  P(r),  and  P^i^r'i).  Exchange  manifests  itself 
more  explicitly  in  the  second  factor  of  the  structural  relation  (7).  We  now  examine 
the  functions  P(rjri)  and  I(rir2r'1)  in  more  detail. 

Consider  that,  due  to  the  disappearance  of  correlations  at  large  distances,  the 
quantity  p2(r l^r'j)  must  reduce  simply  to  ppi(rir'j)  when  r2  — *  oo.  On  the  other 
hand,  consider  that  the  one-body  density  matrix  Pi(i‘i,  r',)  is  known  to  have  the 
structure* 

Pi(riri)  =  /JiD(riri)[Pi(riri)  +  I2(riri)]  (10) 

in  terms  of  the  direct-direct  component  p\£>  of  (9)  and  the  exchange  functions  L\ 
and  L2  (denoted  N\  and  iV2  in  Refs.  2,4).  The  exchange  functions  L\  and  L2  are  in 
turn  known,  both  formally  and  numerically,  within  FHNC  theory.4,5  The  diagrams 
contributing  to  these  functions  through  three-body  cluster  order  are  indicated  in 
Fig.  2.  The  foregoing  considerations,  together  with  the  vanishing  of  the  function 
Zf(rir2r'j)  as  r2  — >  oo,  lead  to  the  identification 

L(riri)  =  Z-Kiqri)  +  L2(rlr[)  .  (11) 

The  structure  of  the  three-point  exchange  function  L{ rji^r'i)  cannot  be  deter¬ 
mined  from  asymptotic  properties;  it  is  instead  inferred  from  the  raw  results  for 
the  cluster  expansion  of  p2(riv2r\),  along  with  the  resolution  (7).  Many  diagrams 
are  seen  to  factorize,  consistent  with  the  following  expression  in  terms  of  irreducible 
diagram  sums  Pap,  PQpy'- 

P(r1r2r,]i)  =  -v~x  /(r)l(r') 

+  f(r)[PCc(r')  +  PdccWr.r,)]  +  l(r')[Pcc(r)  +  P^r^r',)] 

+  /(^r'^fPde^)  +  Pde(r')  +  Pded(rir2r,1)] 

+  Pcdc(**ir2ri)  +  Pcdc(rir2r,1)  +  PCeC(rir2r'1) 

-  ^[Pcc(f)  +  Pdcc(rir2ri)][Pcc(r')  +  P(icc(r,1r2ri)]  .  (12) 

The  two-  and  three-point  irreducible  exchange  functions  Pap  and  Pa/?7  are  clas¬ 
sified  according  to  the  presence  or  absence  of  exchange  lines  at  the  root  points. 
The  category  to  which  a  given  function  belongs  is  indicated  by  its  subscripts  otfi  or 
a/?7,  according  to  the  conventional  scheme':  cc  (circular),  de  (direct-exchange),  dec 
(direct-circular),  ded  (direct-exchange-direct),  ede  (circular-direct-circular),  and  cec 
(circular-exchange-circular).  The  direct-direct  functions  P(r)  and  P(rir2r'j)  fit  in 
the  same  scheme,  with  dd  subscripts  (omitted  for  simplicity).  All  of  these  functions 
may  be  separated  into  their  nodal  and  elementary  parts,  and  the  techniques  of  Fermi- 
hypernetted-chain  theory  invoked  for  their  numerical  evaluation.  It  is  to  be  noted 
that  the  two-point  functions  Pap  arise  already  in  the  theory  of  the  Fermi  one-body 
density  matrix  and  have  been  thoroughly  studied  within  the  FHNC  framework.4,5 

The  requirement  that  the  fully-diagonal  portion  p2( rjr2ri)  coincide  with  p2</(r i2), 
where  g(ri2)  is  the  radial  distribution  function  for  the  assumed  wave  function,  im¬ 
poses  certain  relations  on  the  quantities  P(r),  P(rir2r'j),  P<*/j,  and  PQ/?7,  namely 

—2 P(r)  -  P(rir2rj)  =  Ndd(r)  +  Edd(r )  , 

Pcc(r)  -f  Pdc c(rjr2ri)  =  Ncc(r)  +  Ecc(r)  , 

2  Pde(r)  +  Pded(r\r2r1)  =  Ndc(r)  +  Ede(r)  , 

Pc</c(rir2ri)  =  Ndc(r)  +  Edc(r)  , 

Pce  c(rir2rj)  =  Nce(r )  +  Ecc(r) 


(13) 


The  nodal  ( N )  and  elementary  (E)  diagram  sums  on  the  right-hand  side  are  just 
those  (with  the  corresponding  subscripts)  which  arise  in  the  FHNC  analysis  of  g{r) 
(cf.  Ref.  7). 

FORM  FACTORS 

We  now  determine  the  structure  of  the  generalized  momentum  distribution  func¬ 
tion  n(p,q)  of  (1)  by  exploiting  the  structural  results  (7)  (8),  (9),  (11),  and  (12)  for 
the  two-body  density-matrix  elements  p2(rir2 rj)  appearing  in  Fourier  integral  (6). 
To  arrive  at  a  decomposition  of  n(p,  q)  which  accomplishes  a  clean  separation  of 
contributions  from  differing  physical  processes,  the  function  P2(i‘ir2r'1)  is  first  de¬ 
composed  into  a  part  containing  all  terms  generated  purely  by  two-point  functions, 
and  a  remainder  in  which  the  terms  also  depend  on  the  irreducible  three  point  func¬ 
tions, 

Mri^r'j)  =  p?)(rir2ri)  +  p53)(r1r2r'1)  ,  (14) 

the  notation  being  transparent.  The  last  term  vanishes  if  the  various  three-point 
functions  P(rir2ri)  and  Pap-y(rir2r'j)  are  set  equal  to  zero. 

Next,  we  appeal  to  the  Fermi-hypemet  equations  resulting  from  the  FHNC  anal¬ 
ysis  of  the  one-body  density  matrix.4,5  These  equations  relate  the  bare  correlation 
function  f(r)  to  the  spatial  distribution  functions  defined  by  the  direct  and  nodal 
diagram  sums  XQap(r)  and  NQap(r )  introduced  in  Refs.  4,5,  NQdd(r ),  in  particular, 
being  identified  with  the  nodal  part  of  P(r): 

9Qdd{r )  =  1  +  Fdd{r )  =  1  +  XQdd(r)  +  NQdd(r)  , 

9Qcc(r)  =  1  +  Fcc(r)  -  1  +  XQcc(r)  +  NQcc(r)  , 

9Qde{r)  =  1  +  Fde(r)  —  1  +  XQdc(r)  +  NQde(r)  .  (15) 

The  hypernet  equations  read 


f(r)e  P(r)  =  gQdd(r )  , 

f(r)e-p^Pcc(r)  =  v-ll{r)Fdd(r)  +  Fcc(r)  , 

f(r)e-p^Pde(r)  =  Fde(r)  .  (16) 

Eqs.  (16)  permit  us  to  eliminate  the  bare  correlation  factor  /(r)  from  each  term  of 
the  explicit  expression  for  the  decomposition  (14),  with  the  results 

P22)(rjr2r'1)  =  pp  i  ( r  i  r  j  )gQdd(r)gQdd(r' ) 

+  PPiD(rir\)l(vir[)[gQdd(r)Fde(r')  +  gQdd{r')Fde(r)} 

~  */PPiZ?(rl ri )(t/-1Kr)  -  Fcc(r)][ -  Fcc(r')}  (17) 

and 


P23)(rir2ri)  =  P22)(ri^2r,1){exp[-P(rir2r,1)]  -  1} 

+  PPiD(rir'1)pQ<f(<(r)pQrf(i(r,)exp[-P(r1r2r,1)] 

x  (/(r)Pdcc(r'jr2ri)  +  l{r')Pdcc(r ir2rj) 

/(rj  r\  )Pded(  ri  r2r', )  +  Pcdc(  ri  )  +  Pc2rfc(r,  r2r', ) 

Pocc(r ir2r'j)  -  uPdcc^\^\)Pdcc{^\^\) 

-  i/Pcc(r)Prfcc(rir2ri)  -  t/Pcc(r,)P<fcc(r1r2r,1))  .  (18) 


(Here  r  =  |rj  -  r2|  (as  usual)  and  r'  =  |r'j  -  r2|.) 

Inserting  (14)  into  the  integral  (6),  and  using  the  results  (17)  and  (18),  we  obtain 
the  decomposition 

n(p,q)  =  Mqon(p) 

+  Fdd(q)[n(p)  +  n(\p  -  q\)] 

+  Fde{q)[nDl(p)  +  nDt( |p  -  q| )] 
no[0(kF  -p)-  Fcc(p))[0(kF  -  |p  -  q|)  -  Fcc(|p  -  qj)] 

+  n(2)'(p,q) +n(3)'(p,q)  .  (19) 

Thus,  the  component  n(p,  q)  -  Ar5q0n(p)  of  the  generalized  momentum  distribution 
depending  nontrivially  on  two  momenta  is  written  as  a  sum  of  (i)  separable  contri¬ 
butions  involving  form  factors 

Fa M  =  P I  Fap(r)ei<l'rdr  (20) 

and  either  the  one-body  momentum  distribution 

n(72)  =  v~l  J pi(r)e,p-rdr  ,  (21) 

a  modified  momentum  distribution 

nDl(p)  =  v~l  J  PiD(r)l(r)eip'rdr  ,  (22) 

or  the  the  strength  factor  n0  [second,  third,  and  fourth  terms  of  (19),  respectively],  (ii) 
a  non-separable  integral  n ^  (p,  q)  involving  only  two-point  quantities,  and  (iii)  an¬ 
other  three-point  integral  (p,  q)  generated  from  the  component  (18)  of  p2(r1r2r'1). 
Explicitly,  the  fifth  term  is 

n(2)'(p,q)  =  —  J  I<( r1r2ri)e-,>(ri-ri)e-,'q-(r!-r2)dr1dr2dr'1  ,  (23) 

where 

A'(rir2ri)  =  ppi{ryY\)FQdd{r)FQdd{r') 

+  PPi  D(^ir[)l(rir[)[FQdd(r)Fde(r')  +  FQdd(r')Fde(r)) 

~  vp[piD(rir\)  -  pn0][u~l l{r)  -  Fcc(r)][i/_1/(r')  -  Fcc(r')j  .(24) 

The  contribution  (23)  may  be  reduced  to  a  three-dimensional  integral  in  momentum 
variables. 

In  momentum  space,  the  sequential  relation 

J  p2(rir2r'j)dr2  =  (N  -  l)pi(riri)  (25) 


becomes 


n(P>  q  =  0)  =  (iV  —  l)n(p) 


(26) 


Specializing  (19)  to  q  =  0,  and  substituting  into  (26),  there  results  a  condition 
on  the  form  factors.  If  we  employ  this  condition  in  (19)  itself,  we  may  recast  the 
decomposition  of  n(p,q)  in  the  form 

rt(p,q)  =  (iV-l)<5qon(p) 

+  (1  -  <5qo)i?d<i(a)[n(p)  +  n(|p  -  q|)j 

+  (1  -  8qo)Fde{q)[nDi(p)  +  rc£>/(|p  -  q|)] 

-  n0(l  -  $qO)[0(&F  -  p)  -  -Pcc(p)][^(^F  -  |p  -  q|)  -  ^cc(|p  -  q|)) 

+  (1 -<5q0)n(2),(p,q)  +  (l-<5q0)n(3)'(p,q)  .  (27) 

This  expression  achieves  the  desired  separation  of  contributions  from  the  various 
scattering  processes  underlying  the  generalized  momentum  distribution  function  (cf. 
Ref.  8).  The  first  term  reproduces  the  trivial  result  for  dynamically  and  statistically 
uncorrelated  particles  (except  that  the  momentum  distribution  function  n(p)  appear¬ 
ing  in  (27)  is  that  for  the  fully  correlated  Fermi  system).  The  correlations  prevailing 
in  the  interacting  fluid  permit  the  scattering  of  a  fermion  from  orbital  p  to  another 
orbital  p  —  q,  with  the  intervention  of  a  phonon  to  conserve  momentum.  The  effect 
of  this  process  and  the  corresponding  time-reversed  mechanism  are  described  by  die 
second  term  in  (27).  The  associated  exchange  scattering  effects  are  embodied  in  the 
third  term,  which  is  proportional  to  the  exchange  form  factor  Fde(q)-  To  interpret  the 
fourth  line  of  (27),  we  note  that  if  the  particles  are  noninteracting,  but  the  statistical 
correlations  are  turned  on,  the  totally  uncorrelated  result  (N  —  l)£qon(p)  must  be 
corrected  by  a  Pauli  kinematic  term  —(1  —  6q0)6(kp'  —  p)0(kp  —  |p  —  q|).  For  the 
present  case  of  interacting  fermions,  the  dynamical  correlations,  manifested  in  virtual 
single-particle  excitations  out  of  the  Fermi  sea,  lead  to  tails  on  the  step  distributions 
(the  Fcc  terms).  The  dynamical  correlations  also  produce  an  overall  quenching  of  the 
effect,  through  the  strength  factor  n0  (0  <  n0  <  1).  The  last  line  of  (27)  contains 
terms  of  “higher  order”  which  act  to  correct  the  various  processes  just  considered. 

Upon  taking  the  Bose  limit  of  (27),  we  recover  the  corresponding  decomposition 
of  the  generalized  momentum  distribution  of  a  Bose  fluid  as  derived  in  Ref.  8,  namely 

n(p,q)  =  <5q0(iV  -  l)n(p) 

+  (1  —  <$qo)(6po  +  t>pq)Nn0Fi(q) 

+  (1  ~  £qo)(l  -  £Po)(l  -  <5pq)n'(p,  q)  ,  (28) 


the  non-condensate  portion  being  given  by 
rc'(p,q)  =  Fx(q)[n'(p)  +  n'(\p  -  q\)) 

+  jj  J p[(rir\ )^l (r, r2r, )e-tp^ri “r>  *e"lq'*ri ~ri^dridr2dr\  .  (29) 

Only  the  direct  contributions  to  (27)  survive.  Thus  the  third  and  fourth  lines  are 
to  be  omitted,  and  the  form  factor  Fdd(q)  is  to  be  identified  with  the  function  Fx(q) 
studied  in  the  earlier  paper.  The  function  ^(ri^rj)  entering  (29)  is  specifically 

*i(r, rarj)  =  Fx(r)Fx(r')  +  [1  +  fi(r)J[l  +  F’1(r')]{exp(-P(r1r2r'1)]  -  1).  (30) 

The  three  addends  of  (28)  have  the  following  interpretations  as  contributions  from 
distinct  scattering  processes  taking  place  in  the  many-body  medium  (see  Fig.  3): 
The  trivial  first  term  corresponds  to  the  null  scattering  process  3(a).  The  third  term 


accounts  for  3(b),  in  which  a  boson  is  scattered  from  orbital  p  outside  the  condensate 
to  another  non-condensate  orbital  p  —  q  by  transferring  momentum  tiq  to  a  density 
fluctuation.  The  other  two  processes  shown,  3(c)  and  3(d),  are  described  by  the 
second  line  of  (28)  and  involve  the  zero- momentum  condensate.  They  correspond, 
respectively,  to  creation  of  a  particle  of  momentum  —hq  by  the  condensate  and  to 
absorption  of  a  particle  of  momentum  Tip  into  it,  a  phonon  being  created  in  each  case 
to  conserve  momentum. 


(c)  '  (d) 


Fig  3  Scattering  processes  contributing  to  the  generalized  momentum  distribution  function  ri(p,q) 
of  a  Bose  fluid  with  a  macroscopic  zero-momentum  condensate  (cf.  Eq.  (28)  and  following  text). 


Finally,  we  should  address  the  issue  of  the  practical  validity  of  the  sequential 
relation  (25)  or  (26).  In  the  Bose  limit,  we  may  split  this  relation  into  two  conditions, 
one  applying  at  p  =  0  (the  condensate  condition)  and  the  other  applying  at  finite  p 
(and  involving  terms  smaller  by  a  factor  l/N).  These  conditions  read  l+2F\(q  =  0)  = 
0  and  n ^  (p,0)+n^  (p,  0)  =  0,  respectively.  As  discussed  in  Ref.  8,  they  constrain 
the  choice  of  correlation  factor  in  the  Bose  case.  In  particular,  it  may  be  shown  that 
the  condensate  condition  is  satisfied  for  optimal  Jastrow  correlations,  even  at  the 
HNC/0  level  of  approximation  (i.e.,  neglecting  elementary  contributions),  but  fails 
for  a  generic  /(r).  In  the  Fermi  case,  the  situation  is  rather  different:  the  sequential 
relation  (26)  is  fulfilled  identically  for  any  Jastrow  trial  function,  and  indeed,  even 
if  this  choice  is  extended  to  include  multi-body  correlations  of  Feenberg  type.  The 
universal  satisfaction  of  (26)  is  a  consequence  of  the  presence  of  Pauli  exchange 
correlations,  implied  by  the  Slater  determinant  <£.  Technically,  this  property  may  be 
attributed  to  a  Fermi  cancellation  phenomenon  of  the  type  encountered  in  the  FHNC 
analysis  of  the  static  structure  function  S(q )  (cf.  Ref.  7).  While  formally  exact,  the 
required  cancellation,  and  hence  the  sequential  relation,  is  violated  to  some  degree 
in  any  finite-order  approximation  within  the  HNC/n  sequence. 

NUMERICAL  RESULTS  FOR  SINGLE-MOMENTUM  FORM  FACTORS 

In  correcting  the  impulse  approximation  to  deep-inelastic  scattering  from  the 
helium  liquids  for  final-state  interactions,  Silvei1  has  used  the  simple  approximation 

Pi( fir2ri)  ^  pMor'jMIrj  -  r2|)  (31) 

for  the  required  two-body  density-matrix  elements.  This  approximation  may  be 
tested  microscopically  within  variational  theory,  by  hypernetted-chain  evaluation  of 
the  functions  p2(rir2r'j),  /^(rir'j),  and  g(r). 


In  the  framework  of  the  Jastrow  variational  theory  developed  here  for  Fermi 
systems  of  arbitrary  level  degeneracy,  Silver’s  approximation  (31)  corresponds  to 
certain  rough  first  estimates  of  various  quantities  appearing  in  the  contributions  (17) 
and  (18)  to  (14),  en  route  to  the  construction  of  n(p,  q)  via  (27).  In  more  detail,  it 
is  equivalent  to  the  replacements 


ff,l(rir,j)  —  />i  £?(r  i  r#,  )/(r,  r, )  ,  (32) 

FdJ(r)  +  Fd,(r)  ct  9(r)  -  1  ,  Fdd(r')~  0  ,  Fde(r')~  0  (33) 

v-'l(r)-F«(r)x  0  ,  (34) 

and  to  neglect  of  the  “pure-three-point”  contribution  )  of  Eq.  (18).  In 

momentum  space,  relation  (32),  the  first  member  of  (33),  and  relation  (34)  become, 
respectively, 

n(p)~nDl(p)  ,  (35) 

Fdd(q)  +  Fde(q)  ~  S(q)  -  1  ,  (36) 

v~x0(kp  —  p)  —  Fcc(p)  —  0  ,  (37) 

where  S(q)  is  the  static  structure  function  corresponding  to  the  radial  distribution 
function  g(r).  As  may  be  seen  by  taking  the  Bose  limit,  the  corresponding  replace¬ 
ments  implied  in  the  Bose  case  are 

F,(r)~5(r)-1  ,  F,(r')~0  ,  (33') 

F,(?)~%)-1  ,  (36') 

and  neglect  of  the  three-point  function  F^ri^rj). 

The  most  striking  feature  of  Silver’s  approximation  is,  of  course,  its  violation  of 
time-reversal  invariance,  evident  in  the  asymmetric  treatment  of  at  least  one  of  the 
pairs  Fdd(r),  Fdd(r')  and  F<f«(r),  Fde{r')  by  (33),  and  of  the  pair  Fi(r),  Fj(r')  by 
(33’). 

To  make  a  quantitative  judgment  of  the  efficacy  of  approximations  (32)-(37), 
(33’),  and  (36’),  we  have  calculated  the  distribution  functions  and  form  factors  en¬ 
tering  these  relations,  for  liquid  4  He  and  liquid  3  He  at  their  respective  equilibrium 
densities  (p  =  0.0218  A-3  and  p  =  0.01658  A-3).  For  the  Bose  liquid  we  have  used 
a  Jastrow  correlation  factor  /(r)  optimized  by  a  paired-phonon  analysis9;  a  Schiff- 
Verlet  form10  /(r)  =  exp[— (b/r)5 /2]  was  chosen  in  the  Fermi  case,  with  b  =  2.9547  A 
(cf.  Ref.  11).  The  numerical  evaluations  were  carried  out  in  the  Bose  or  Fermi 
hypernetted-chain  approximation  (HNC/0  or  FHNC/0)  in  which  the  elementary- 
diagram  contributions  to  the  various  quantities  are  set  zero.  This  approximation 
should  be  adequate  for  the  immediate  task  of  testing  the  ansatz  (31). 

Our  results  are  summarized  in  Figs.  4-6. 

Figure  4  compares  the  form  factor  Fj  (5)  (solid  curve)  with  the  overshoot  S(q )  —  1 
of  the  static  structure  function  (dashed  curve),  testing  (36’)  (or  the  .first  member  of 
(33’)]  in  liquid  4He.  The  two  functions  have  roughly  similar  shapes,  but  depart 
substantially  at  small  q.  At  q  —  0,  the  numerical  result  for  F\{q)  goes  to  the  correct 
asymptotic  result  —1/2  (as  it  should,  since  the  ‘condensate’  condition  derived  from 
the  sequential  relation  must  be  met  in  HNC/0  when  /(r)  is  optimized).  On  the  other 
hand,  the  estimate  (36’)  deviates  from  the  correct  limit  by  a  wide  margin,  since  5(0) 
must  vanish  for  the  optimized  f(r). 


Figures  5  and  6  address  the  situation  in  liquid  3  He.  The  shortcomings  of  the 
estimate  (36)  (or,  equivalently,  the  first  member  of  (33)]  are  revealed  in  Fig.  5,  which 
shows  S(q )  —  1  and  the  sum  of  form  factors  Fj^q)  +  Fje(q).  Within  the  FHNC/0 
approximation,  the  departure  from  assumption  (36)  is  clearly  exposed  and  again  is 
particularly  apparent  at  small  momenta.  The  exclusion  principle  alters  the  <7  =  0 
limits  seen  in  the  Bose  case  (Fig.  4),  by  virtue  of  the  Fermi  cancellation  effect  men¬ 
tioned  in  the  preceding  section.  For  a  Fermi  system  described  by  a  Siater-Jastrow 


Fig.  4.  Test  of  (36’).  Form  factor  F\(q)  for  creation  of  a  particle  out  of  the  condensate,  in  HNC/0  (solid 
curve)  and  Silver’s  (dashed  curve)  approximations,  for  liquid  4He  at  equilibrium  density,  described 
by  an  optimized  J  as  trow  wave  function.  Static  structure  function  of  (36’)  was  evaluated  in  HNC/0 
approximation. 


Fig  5.  Test  of  (36)  Sum  of  form  factors  Fdd(q)  -f  Fdc(q)  (solid  curve)  compared  with  overshoot  of 
static  structure  function  ?(q)  (dashed  curve),  for  liquid  3ile  at  equilibrium  density,  described  within 
a  Jastrow-FHNC/O  approximation. 


wave  function,  tin*  cancellation  phenomenon  guarantees,  at  q  =  0,  the  properties 
Xdt{q)  =  -V -ct(q)  =  -1  and,  consequently,  S(q)  =  0  and  F<u(q )  +  Fdc(q)  =  0.  The 
latter  properties  are  (approximately)  reflected  in  our  numerical  results  for  S(q)  and 
Fdd(q)  +  Fjc(q).  However,  one  does  see,  in  Fig.  5,  slight  deviations  from  the  cor¬ 
rect  limiting  value  of  zero,  which  result  from  use  of  the  FHNC/0  approximation. 
The  standard  FHNC  approximants  (/0,  /4,  etc.)  are  known  to  disobey  the  Fermi 
cancellation  rules  as  a  result  of  the  neglect  or  inconsistent  treatment  of  elementary 
diagrams.' 

Fig.  6  shows  the  momentum  distributions  n(p)  and  noi(p)  and  tests  the  replace¬ 
ment  (35)  (equivalent  to  02)].  The  strength  factor  associated  with  both  of  distribu¬ 
tions  is  n0  =  0.2212.  The  two  functions  are  seen  to  have  very  similar  behavior,  but 
their  magnitudes  differ  typically  by  10-15%. 


Fig.  6.  Test  of  (35)  and  (37).  Momentum  distribution  functions  n(p)  (solid  curve)  and  noi(p)  (dot- 
dashed  curve)  of  liquid  3He  at  equilibrium  density,  described  within  a  Jastrow-FHNC/O  approxima¬ 
tion.  Dashed  curve  shows  circular-exchange  function  Fcdp). 


According  to  (37),  the  Pauli  exclusion  corrections  to  n(p,  q)  of  the  circular  type, 
involving  u~lQ{kp  —  p)  —  Fcc(p),  are  entirely  ignored  in  Silver’s  treatment  -  even  the 
trivial  kinematic  statistical  effect  of  the  first  term  is  absent.  The  form  factor  Fcc(p), 
evaluated  in  FHNC/0  approximation,  is  shown  in  Fig.  6  as  the  dashed  line.  This 
function  vanishes  inside  the  Fermi  sea,  jumps  to  a  height  of  about  0.5  at  the  Fermi 
surface,  and  decreases  slowly  in  magnitude  with  further  increase  of  the  wave  number 
p.  In  general  one  may  therefore  expect  such  statistical  effects  to  be  important.  On 
the  other  hand,  we  note  that  their  net  contribution  to  u(p,  q)  of  (27)  is  proportional 
to  the  strength  factor  n0,  which  is  only  about  0.2  in  liquid  3 He. 

In  brief,  an  initial  application  of  the  microscopic  theory  of  the  density-matrix 
elements  /?2(iT  ri )  has  documented  significant  quantitative  deficiencies  of  the  simple 
estimate  proposed  by  Silver.  We  should  point  out,  however,  that  the  final-state 
corrections  evaluated  in  Silver’s  theory  of  deep-inelastic  scattering  at  high  momentum 
and  energy  transfers  may  be  insensitive  to  the  errors  we  have  noted,  and,  in  particular, 
to  the  behavior  of  the  sum  of  form  factors  Fdd(q)  +  Fje(q)  at  small  q.  This  possibility 
is  currently  under  investigation.12 


FURTHER  WORK 

We  have  not  reported  numerical  data  on  the  non-separable  terms  n(2^(p,  q)  and 
(p,q)  of  (27),  nor  on  the  integral  term  in  (29)  involving  the  three-point  func¬ 
tion  Fi(r]r2r,1).  While  these  more  complicated  objects  may  all  be  calculated  from 
quantities  generated  in  the  FHNC/0  or  HNC/0  treatment,  their  detailed  evaluation 
will  be  deferred  until  a  scaling  or  interpolation  procedure13,14  has  been  implemented 
for  the  incorporation  of  elementary-diagram  corrections.  At  the  same  time,  the  Jas- 
trow  ansatz  will  be  supplemented  by  triplet  correlations.  Work  in  these  directions  is 
currently  in  progress.15 

It  is  worth  mentioning  that  some  HNC/0  results  for  the  quantity  F\  (r  i  rorj )  were 
presented  in  Ref.  8.,  results  which  demonstrate  the  failure  of  the  HNC/0  treatment 
to  fulfill  the  condition  P2(rir2ri)  —  P2g{r 12)  and  thereby  provide  another  reminder 
of  the  necessity  of  including  effects  of  elementary  diagrams. 

While  our  derivation  of  the  decomposition  (27)  for  the  Fermi  generalized  momen¬ 
tum  distribution  was  predicated  on  the  Slater- Jastrow  choice  for  the  wave  function 
the  corresponding  decomposition  (28)  for  the  Bose  system  (derived  here  as  a  lim¬ 
iting  case  of  the  Fermi  result)  was  originally  obtained  for  very  general  'F,  including 
the  exact  ground-state  wave  function.8  Indeed,  Ref.  8  contains  a  general  asymptotic 
analysis  of  the  full  Bose  two-body  density  matrix  in  configuration  space,  as  well  as 
the  restricted  version  p2(r l^r'j)  considered  here. 

A  more  detailed  presentation  of  some  of  the  results  of  this  contribution  may  be 
found  in  a  longer  article.16  Further  studies  along  the  same  lines  will  be  concerned 
with  the  full  Fermi  two-body  density  matrix. 
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The  classic  Cooper  electron-pair  problem  is  generalized  via  a  Fermi 
sea  having  two  concentrj.c  surfaces,  rather  than  the  familiar  sphere  of  a 
perfect  (interactionless)  fermion  gas.  Substantially  tighter-bound  pairs 
are  obtained  for  fixed  phonon  coupling  with  the  BCS  model  interaction,  no 
matter  how  weak.  This  will  admit  increased  transition  temperatures  for 
superconductivity  in  the  BCS  theory  and  beyond,  and  suggests  that  designing 
materials  having  multiply-connected  Fermi  surfaces  with  maximal  interior 
area  will  yield  larger  T  values. 

C 


INTRODUCTION 

The  Cooper  electron  pair  problem  is  vitally  central  to  the 
Bardeen-Cooper-Schrieffer  (BCS)  microscopic  theory  of  superconductivity  . 
It  consists  in  a  Schrodinger  equation  for  two  particles  of  opposite  momenta 
and  spins  being  scattered  by  a  particular  attractive  interaction  only  into 
one-electron  plane  wave  (PW)  states  just  above  the  spherical  Fermi  sea 
occupied  by  the  PW  orbitals  of  the  other  N  -  2  background  electrons  in  the 
N-electron  system.  PW  states  are  used  for  simplicity,  instead  of  Bloch 
states.  The  well-known  resulting  eigenvalue  equation  for  the  pair  energy 
E  is  then 
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where  E(k)  =  h  k  /2m  are  the  unperturbed  single-particle  energies  and  V  >  0 
is  the  strength  of  the  effective  attractive  electron-electron  interaction 
induced  by  the  electron-phonon  coupling.  This  is  nonzero  only  within  a 
very  thin  shell  of  thickness  hw  above  the  Fermi  surface  of  energy  E  = 

D  F 

hKV2m,  k,  being  the  Fermi  sphere  radius;  it  is  the  so-called  "BCS  model 
interaction".  Here,  hw  is  the  maximum  energy  possible  for  an  ionic 

-3  -2 

lattice  phonon,  and  is  typically  10  to  10  times  smaller  than  Ep.  The 

prime  over  the  summation  sign  means  restriction  to  those  PW  states  such 

that  E  <  E(k)  <  E  +  hw  ,  i.e. ,  to  states  that  are  unoccupied.  The 

1  /2 

integral  involves  the  free-electron  density  of  states  g(§)  «  &  ,  which 

can  in  turn  be  factored  out  from  the  integral  as  a  constant  g(S^)  due  to 
the  smallness  of  hw^/E^.  This  leaves  an  elementary  integral  to  be 
performed  that  gives  a  logarithm.  Solving  for  the  eigenvalue  E  then  yields 


2hw 
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Putting  e  =  E  /2E  and  v  a  hw  /2E 

°  o  o  F  D  F 

to  read 


the  first  equation  can  be  rearranged 
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The  lhs  (a  constant  for  fixed  coupling  and  density  of  states)  and  rhs  of 

_3 

this  equation  are  displayed,  for  the  special  but  typical  case  2v  =  10  ,  in 
Figure  1  as  function  of  e  s  1  -  A  /2E  ,  where  A  is  the  pair  binding 

o  F  o 

energy.  We  see  that  6  differs  very  little  from  unity  (and  hence  A  /2E 
from  zero)  for  all  but  the  largest  values  of  the  coupling  parameter  X  = 
g(EF)V/2.  Since  (3)  results  from  solving  a  two-body  Schrodinger  equation 
(in  momentum  representation),  its  validity  is  not  limited  to  weak  coupling, 
unless  regarded  as  a  Bethe-Goldstone  equation  (ladder  approximation) 
treatment  of  the  many-fermion  problem,  which  is  exact  only  for  small  k^X. 

2 

Also  restricted  to  weak  coupling  is  the  BCS  many-electron  theory  .  In 
this  formalism  a  temperature-dependent  energy  gap  A(T)  emerges,  which  for  T 
=  0  is  found  to  be 
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"3 

Fig.  1.  Graphical  solution  of  Eq.  (3)  for  2v  -  10  and  several 
values  of  the  coupling  parameter  A  as  indicated. 
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i.e.,  identical  to  the  weak-coupling  limit  of  the  Cooper  pair  binding 
energy  A  in  (2).  The  "normal  metal"  to  "superconductor"  transition 

O 

temperature  T  is  then  determined  by  the  vanishing  of  A(T),  or  A(T  )  =  0, 

C  C 

whose  solution  is  the  celebrated  relation 
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where  %  ^  0.577  is  the  Euler  constant.  Combining  (4)  and  (5)  we  have 


T  =  1.  13  0  e 

c  D 
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(6) 
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where  8d  a  hu./k  Js  the  Debye  temperature.  For  elemental  superconductors 

measured  T  values  range  from  a  very  low  3.25  X  10-4  K  for  Rh  to  9.26  K  for 

Nb.  Since  0d  is  proportional  to  an  appropriately  averaged  sound  speed  c 

—  -1/2 

for  the  lattice,  and  since  c  «  M  ,  with  M  the  mass  of  the  lattice 

Ion  ion 

-1/2 

ions,  we  see  that  T  «  M  .  This  is  known  as  the  isotope  effect, 

cion 

observed  for  many  elemental  superconductors  but  not  for  the  new 

2 

copper-oxide  high-T  materials.  Since  0  ~  10  K,  (6)  severely  limits  T 

c  D  c 

3 

to  a  few  degrees  K  with  acceptable  values  of  X.  Other  T  formulas  , 

C 
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beginning  with  the  MacMillan  formula  based  on  strong-coupling 
Migdal-Eliashberg  theory5,  give  values  of  T  as  high  as  around  40  K. 

C 

Indeed,  a  recent  realistic  tight-binding  calculation  with  the  Eliashberg 
equations  gave  Weber6  T  values  between  30  to  40  K  for  the  copper-oxide 
superconductors  La  (Ba,Sr)  CuO  ,  which  have  empirical  T  values  in  the 

2-x  X  4  c 

range  30  to  36  K.  But  for  YBa  Cu  0  with  observed  T  =  95  K  Weber  and 

2  3  7  c 

Mattheiss,  using  similar  techniques,  were  not  able  to  extract  a  larger 
than  about  30  K.  Furthermore,  the  observed  value  of  the  exponent  a  in  the 

formula  T  «  M~a  is  about  0.18  for  the  lanthanum-based  copper  oxide  just 

c  ion 

cited,  and  0  for  the  yttrium-based  copper  oxide. 

We  show  how  the  restriction  to  low  T  's,  exemplified  by  (6)  and 
i/A  c 

traceable  to  the  e-1  factor  in  the  Cooper  pair  problem  (2),  can  be 
surmounted  through  a  generalized  Fermi  sea.  And,  moreover,  a  negligibly 
weak  isotope  effect  is  then  possible,  even  with  a  phonon-mediated 
interaction. 

The  fact  that  (quasi)  electron-pair  binding  occurs  for  arbitrarily 
weak  attraction — as  well  as  the  non-perturbative  (essential  singularity  in 
V)  property  of  the  energy  A  needed  to  break  up  a  Cooper  pair - occurs  even 

O 

8 

for  one-body  potential  well  problems  in  two-dimensions  (2D).  Recalling 
that  g (g)  =  const  in  2D,  it  is  then  clear  why  the  Cooper  pair  problem 
behaves  like  a  2D  quantum  system,  a  fact  which  may  be  pivotal  in 
understanding  the  superconductivity  of  the  layered  copper-oxide  ceramics. 


ABNORMAL  OCCUPATION 

The  all-important  Fermi  sea  assumed  in  the  Cooper  problem  is 
appropriate  to  the  perfect  Fermi  gas  Slater  determinant  ground  state  wave 
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function  for  N  particles  enclosed  in  a  volume  Q, 


,  n-i/2  ,  ,  \r.-i/2  lk  •  r 

$>  =  (N!  )  det  o  Q  e  1  j 

k  •* 

with  i,j  =  1,2,  ...N,  and  n°  =  0(kF  -  k)  the  unit  step  function,  0(x)  =  i[l  + 
sgn(x)].  For  an  interacting  system  the  most  general  occupation  scheme 
consistent  with  the  Pauli  principle,  however,  is  merely 

nJc  =  0  or  1,  ^  n^,  =  N.  (8) 

k 

9 

We  have  raised  the  general  question  of  what  the  optimum  scheme  might  be 
for  a  non-ideal,  fully-interacting  many-fermion  assembly  in  any  dimension. 
After  all,  Fermi  surfaces  for  many  materials  are  known  which  contrast 
drastically  with  the  familiar  Fermi  sphere.  Overhauser10,  for  example,  has 
considered  multiply-connected  Fermi  seas  associated  with  both  charge  and 
spin  density  wave  states.  These  states  break  the  translation  symmetry  and 
will  not  concern  us  here.  Indeed,  as  early  as  1950  Frohlich11  had  already 
contemplated  a  departure  from  the  simple  Fermi  sphere,  with  ordinary  PW 
orbitals.  He  took  a  spherical  shell  concentric  upon,  but  disjointed  from, 
an  interior  sphere,  of  occupied  electronic  states,  and  used  this  with  his 
new  electron-phonon  hamiltonian  in  (second-order)  perturbation  theory  to 
find  a  lower  energy  state,  if  the  coupling  exceeded  a  certain  critical 
value.  This  behavior,  however,  sharply  differs  from  the  empirical  fact 
that  superconducting  critical  temperatures  can  be  arbitrarily  small,  a 
drawback  eventually  circumvented  by  the  (infinite-order)  treatment*  of 
Cooper  with  the  standard  Fermi  sea.  More  recently,  in  an  infinite 
meson-nucleon  system  a  transition  from  a  Fermi  sphere  to  a  Fermi  "shell" 
[cf.  below,  Eq.  (10)]  distribution  has  been  interpreted  as  a  possible  phase 
transition  from  "nuclear"  to  "quark"  matter  within  a  relativistic 
Hartree-Fock  (HF)  treatment  .  The  (first-order)  transition  is  signalled 
by  a  large  jump  at  high  densities  in  the  resulting  low-temperature  specific 
heat.  Finally,  abnormal  occupation  in  a  finite  nucleus  has  been  discovered 
in  24Mg,  in  constrained  HF  calculations13  (with  good  total  angular  momentum 
and  its  projection),  using  a  realistic  effective  two-nucleon  interaction 
based  on  the  Reid  soft-core  potential. 

Within  the  (nonrelativistic)  plane-wave  Hartree-Fock  (PWHF) 


5 


approximation  with  many-fermion  hamiltonian  H  we  originally  showed9  that 
for  a  simple  one-dimensional  system  under  a  sufficiently  attractive  (but 
non-collapsing  in  the  thermodynamic  limit),  short-ranged,  two-body 
interaction  v  ,  lower  HF  total  energies 


^HFtnJ  5  <  > 


=  V  E(k)n  +  =■  V  (<k  k  |v  |k  k  >  -  <k  k  |v  |k  k  >}  n  n 
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emerge  for  some  particle  densities  N/fl  =  kj./3n  ,  with  a  generalized  Fermi 
sea  defined  by 


n  =  e(k  -  0k  )e(yk  -  k) 

K  r  F 

r3  -  03  =  1,  7  >  0  st  0,  7*1. 


(10) 


This  corresponds  to  a  spherical  shell  in  k-space  of  inner  and  outer  radii 

0k  and  jk  ,  respectively.  In  other  words,  it  was  established  that 

8  [n  ]  <  8  [n  ]  for  a  range  of  densities,  and  in  a  manner 

reminiscent  of  a  (first-order)  gas-liquid  phase  transition.  Any  lowering 

of  the  total  HF  energy  at  fixed  density  comes  entirely  from  the  last 

(exchange)  term  in  (9),  since  the  first  (kinetic  energy)  term  can  only 

increase  for  any  n^  other  than  a  spherical  Fermi  sea,  while  the  second 

(direct)  term  is  unchanged  if  the  interaction  is  local,  as  then  that  matrix 

element  is  independent  of  the  summation  indices  for  PW  states.  The  search 

for  lower-energy,  abnormally  occupied  Slater  PW  determinants  for  a  wide 

14 

variety  of  pair-interaction  cases  was  subsequently  extended  to  three 
dimensions  (3D),  and  to  a  much  larger  class  of  abnormal  occupation  schemes 
n^.  Moreover,  starting  from  a  very-many-shell  structure  modeling  the 
general  distribution  (8),  random-search  and  random-walk  numerical 
techniques14  established  the  "single-shell"  distribution  (10)  as  the 
optimum  one  in  numerous  cases.  In  still  further  work15  it  was  proved,  for 
example,  that  any  hard-core  plus  square-well  two-body  potential  prefers 
abnormal  occupation,  no  matter  how  weak  the  attraction.  This  potential  has 
been  employed  to  model  liquid-  He  semi-realistically.  Finally,  several 

17 

mar.y-boson  systems  were  also  found  which  prefer  abnormal  occupation. 


Summarizing,  even  at  the  PWHF  level  of  approximation,  Fermi  seas  more 
general  than  the  usual  spherical  sea  are  favored  for  sufficiently  strong 
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interparticle  coupling.  A  possible  physical  rationale18  for  this  is  the 
appearance,  as  coupling  is  increased,  of  "particle  clusterings"  of  some 
kind,  since  abnormal  occupation  precludes  small- k  states,  meaning 
suppression  of  particle  orbits  with  larg<.  “>stial  extensions. 


TIGHTER-BOUND  COOPER  PAIRS 

Robustly  tighter  bound  Cooper  pairs  are  possible,  for  any  coupling 

strength,  with  abnormal  occupation.  We  suggest  that  when  suitably 

19  20  21 

incorporated  into  the  BCS-Bogoliubov  -Anderson  -Gor’kov  -Migdal 

-Eliashberg5  formalisms,  this  may  yet  provide  a  comprehensive  understanding 
of  both  low-  and  high-T  superconductivity,  perhaps  in  terms  of  the  phonon 

C 

mechanism  alone.  The  term  "normal"  is  here  used  in  the  sense  of  Landau’s 

Fermi  liquid  theory,  since  a  key  assumption  for  a  first-principles 

22 

derivation  by  Klein  of  Landau’s  theory  is  that  the  lowest-energy 

single-particle  states  be  occupied,  i.e. ,  that  the  familiar 
finite-temperature  Fermi  distribution  holds  for  the  single-particle 
spectra,  instead  of  some  other,  more  general,  distribution. 


Guided  by  our  PWHF  studies9’14’17  we  employ  the  model  (10)  merely  as 

an  illustration.  This  particular  occupation  scheme  is  a  definite  step 

beyond  the  perfect  Fermi  gas  picture,  and  suffices  to  uncover  an 

instability  in  the  Fermi-sphere-induced  Cooper  pair.  The  two  surfaces  in 

(10)  are  situated  at  energies  E  =  p  =  (y  -  1)  E  £  0  and  E  =  y  E  £ 

"IF  F  2  F 

E  .  The  parameter  y  is  presumably  characteristic  of  the  fully-interacting 
many- body  system,  and  can  later  be  fixed,  within  a  BCS-like  theory,  e.g. , 
variationally,  with  respect  to  the  total  superconducting  ground-state 
energy.  The  integral  in  (1)  now  becomes 


1  -  V 


|g(E) 
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28  -  E 

E  -hw 


g<y 


E  +hw 
,  2  D 

d& 

28  -  E| 


D 


(11) 


provided  that  E^^  is  significantly  larger  than  ha>D  (as  was  verified  post 
hoc).  Scattering  is  now  allowed  by  the  BCS  model  interaction  in  the 
vicinity  of  both  surfaces.  Performing  the  integrals  makes  solving  for 
^-2/g(Ep)v  easier  than  for  the  new  eigenvalue  E.  Putting  E/2E^  =  e  and 
ha>D/2EF  s  i>,  Eq.  (11)  leads  to  a  transcendental  equation  for  e  which,  in 
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the  3D  application  presented  here,  is 


-l/A  _  -2/g(E  )V 

e  =  e  F 

Both  quantities  in  the  square  brackets  must  be  non-negative  to  ensure  that 
the  lhs  of  (12)  be  real. 

We  have  carried  out  a  numerical  search  for  values  of  €  and  y 
satisfying  (12),  such  that  y  >  1  and  >  €  =  1  -  A/2Ef.,  with  A  positive, 
for  a  wide  range  of  (X,  w)  values.  However,  a  better  feeling  for  the 
solutions  of  our  model  is  obtained  from  Fig.  2,  where  both  sides  of  Eq. (12) 
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J 

- 1 

CNJ 
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Fig.  2.  Typical  case  illustrating  graphical  solution  of  Eq.  (12). 

“3 

are  plotted  for  the  case  specified  by  X  —  0.5,  2v  =  10  ,  and  y  -  1.2.  The 

lhs  then  equals  *  0.135.  The  rhs  of  (12)  gives  rise  to  the  full  thick 
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curves  shown;  the  dashed  curves  correspond  to  negative  quantities  inside 
both  square  brackets  of  (12)  and  so  are  not  relevant.  The  two  zeros 
(poles)  are  designated  by  dots  (crosses)  on  the  €-axis.  They  are  labeled 
by  their  ~al  values,  for  any  y  £  1  and  Zv,  which  for  the  specific  case 
just  cited  are  respectively  equal  to  0.808259,  0.809259,  1.440  and  1.441. 
The  dot-to-cross  distances  have  been  exaggerated  somewhat  in  the  figure, 
for  clarity.  The  lowest  energy  eigenvalue  occurs  from  the  projection  of 
the  intersection  marked  "NEW"  and  indicated  by  an  open  circle:  it  is  a 
(quasi)  bound  state  since  e  <  1  and  moreover,  by  inspection,  e  <  e  (=  1  - 

O 

10-3  =  0.999).  The  intersection  marked  "OLD"  and  indicated  by  an  open 
triangle  is  now  an  unbound  level  in  the  pair  continuum,  and  becomes  the 
ordinary  Cooper  bound  state  (2)  or  (3)  as  y  -»  1,  when  the  open  circle 
intersection  is  pushed  leftwards  outside  the  6  >  0  region.  The  reason  for 
tighter-bound  Cooper  pairs  with  abnormal  occupation  is  almost  trivial:  an 
interior  Fermi  surface  allows  for  lower  kinetic  energy  electron  pairs  with 
roughly  the  same  potential  energy*  As  in  the  usual  Cooper  pair  case,  the 
new,  tighter-bound  Cooper  pair  solution  will  survive,  no  matter  how  weak 
the  coupling  V,  since  smaller  coupling  merely  lowers  the  horizontal  line 
marked  "Ihs".  More  interesting,  however,  is  that  these  tighter-bound 
Cooper  pairs  can  lead  to  transition  temperatures  T  with  a  robust 

C 

contribution  of  order  T  (  ~  10  to  10  K),  in  addition  to  a  much  weaker 

F  2 

" isotooe-effect"  contribution  of  order  0  (  ~  10  K),  as  we  now  discuss. 

4  n 


HIGH-TEMPERATURE  SUPERCONDUCTIVITY 

In  the  weak-coupling  limit  (A  =  g(EF)V/2  «  1)  the  intersection  marked 
"new"  in  Fig.  2  will  correspond  to 

1  -  A/2Ef  =  e  =  (?3  -  1  i273  -  Zv  -  v,  (13) 

with  0  <  v  «  1.  Inserting  the  last  member  of  this  into  (12)  and  expanding 

about  7)  =  0  gives  7)  ^  Zve~1/^.  Recalling  (5),  (13)  leads  to 

T  «  1.13  (1  -  £2)T  +  1.13  (1  -  e_1/^X)  0n  +  ....  (14) 

c  F  D 

3  1/3 

where  (3  =  (y  -  1)  is  the  radius  (in  units  of  kp)  of  the  inner  Fermi 
surface  of  the  abnormal  distribution  (10).  We  found  that  typically  (3 

ranged  from  about  0.7  to  about  0.9,  so  that  1  -  /32  ranges  from  0.2  to  0.5. 
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Thus  values  well  in  excess  of  room  temperatures  (~300  K)  are  possible 

without  invoking  stronger  electron-phonon  coupling ,  nor  unconventional 

interaction  mechanisms.  The  latter,  however,  are  not  ruled  out  in 

“1/2 

principle.  We  also  note  that  in  (14)  only  9  depends  on  M  so  that  the 

D  Ion 

2  3 

isotope  effect,  if  present,  will  be  10  to  10  times  weaker,  since  0  /T 
10-3  to  10"2. 

CONCLUSIONS 

Tighter-bound  Cooper  pairs  will  arise  by  generalizing  the  assumed 
Fermi  sea  for  the  occupied  background  electrons,  without  invoking  either 
stronger  electron-phonon  coupling  in  the  BCS  interaction  nor  unconventional 
interaction  mechanisms.  The  calculations  performed  in  3D  can 

straightforwardly  be  done  in  any  dimensionality  as  only  the  form  of  g(S) 
(which  is  absorbed  into  A  anyway),  and  the  £  -  r  relation  (10),  will 
change. 

Multiple-surfaced  Fermi  topologies  can  lead  to  T  ’  s  tied  not  to  the 

C 

2 

Debye  temperature  0  (  ~  10  K)  as  in  normal-occupation  BCS  theory,  but 

4  5 

rather  to  the  Fermi  temperature  Tf  (  ~  10  to  10  K).  They  also  lead  to 

the  absence  (or  extreme  weakness)  of  an  isotope  effect,  but,  as  before, 
this  neither  proves  nor  disproves  the  presence  of  an  electron-phonon 
interaction  mechanism. 

If  more  were  known  about  multiply-connected  Fermi  surfaces  and  how  to 
manipulate  them,  perhaps  one  could  design  materials  (most  probably 
compounds)  having  interior  Fermi  surfaces  which,  allowing  for  Cooper  pairs 
with  smaller  kinetic  energies,  would  result  in  higher  T  values. 
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Abstract:  The  nuclear  shell  model  exists  now  for  forty  years.  The  notion  of 
independent  particle  motion  which  is  at  the  basis  of  this  simple  picture  of  nuclei  has 
recently  been  challenged  by  new  experimental  data.  Nevertheless,  a  considerable 
amount  of  experimental  data  has  successfully  been  interpreted  over  the  years  in  terms  of 
this  picture.  On  a  microscopic  level,  the  success  of  the  shell  model  has  been  linked  to 
the  strong  short-range  repulsion  in  the  nuclear  interaction  by  introducing  the  concept  of 
the  healing  of  the  relative  wave  function  of  two  nucleons.  Recent  calculations  of  the 
effective  interaction  in  nuclear  matter  are  discussed  which  demonstrate  that  this  concept 
of  healing  is  not  valid  but  instead  particles  scatter  at  all  energies  in  the  medium.  It  will 
be  argued  that  the  success  of  the  shell  model  picture  should  instead  be  linked  to  the 
validity  of  the  quasi-particle  picture  as  used  by  Landau  for  liquid  3He.  Calculations  of 
the  spectral  function  for  160  are  discussed  which  demonstrate  the  essential  difference 
between  a  finite  nucleus  and  nuclear  matter.  The  finiteness  of  a  nucleus  ensures  that  a 
relevant  single-particle  basis  has  shell  structure,  implying  that  there  is  an  energy  window 
in  nuclei  related  to  the  difference  between  particle  and  hole  energies  in  which  this 
quasi-particle  picture  can  be  applied. 


1.  INTRODUCTION 

The  simplest  version  of  the  shell  model  considers  nucleons  bound  in  some 
attractive  well  moving  independently  from  each  other.  This  well  can  be  thought  of  as 
being  generated  by  the  average  interaction  a  nucleon  experiences  with  the  other 

’This  research  was  supported  by  the  Condensed  Matter  Theory  Program  of  the  Division  of  Materials 
Research  of  the  U.S.  National  Science  Foundation  under  Grant  No.  DMR-85 19077  (at  Washington 
University)  which  also  provided  computer  time  for  the  calculations  which  were  partly  perfonned  aljhc 
JM^urgh^upcrcomputer  Center. ^AftraVcl;  grantTroim  the  U.S.  Army  Rese^-’^^.is^gbt|ftiily|sbk^‘ 
i  nowdedgedr  Additional  support  wds" provided  in  par’t'by  NATO  under  Grant  No.  RG  85/0684  and  in  part' 
Dy*CAYCfr  Gram  No.  PB85-0072-C02-00  (Spain). 


nucleons  in  the  system.  Since  this  average  interaction  is  obviously  attractive,  nucleons 
can  be  assigned  single-particle  (sp)  quantum  numbers  of  a  spherical  well.  When  an 
average  spin-orbit  interaction  is  included,  this  nuclear  shell  model  can  account  for  all 
observed  magic  numbers  as  well  as  a  number  of  other  experimental  data.  This  is 
essentially  the  picture  as  it  has  been  published  40  years  ago.1,2  It  has  been  very  difficult 
to  link  this  mean  field  picture  of  the  nucleus  and  the  strong  repulsion  present  in  the 
nucleon-nucleon  (NN)  interaction  on  a  microscopic  level.  Historically,  Brueckner3 
suggested  to  convert  the  bare  interaction  into  a  renormalized  interaction  (G  matrix) 
which  includes  any  number  of  interactions  (ladder  diagrams)  between  the  panicles  in  the 
medium,  similar  to  the  corresponding  T  matrix  in  free  space.  This  G  matrix  is  obtained 
by  solving  the  Bethe-Goldstone  equation.4  This  interaction  has  played  an  important  role 
in  the  efforts  to  reconcile  the  strong  interaction  at  short  distance  with  the  tranquil  picture 
of  nucleons  moving  in  simple  sp  orbits.  In  ref.  5  the  concept  of  healing  of  the  relative 
wave  function  is  introduced  as  an  explanation  for  the  validity  of  the  shell  model.  This 
healing  property  results  from  the  fact  that  the  conventional  solution  of  the  Bethe- 
Goldstone  equation  does  not  allow  scattering  for  energies  below  2eF  as  a  consequence 
of  the  Pauli  principle  and  is  therefore  purely  real. 

One  of  the  purposes  of  this  contribution  is  to  show  that  this  analysis  must  be 
improved  and  extended.  This  new  analysis  will  rely  on  explicit  calculations  which  have 
been  performed  recently  for  nuclear  matter.6  In  this  work  the  short-range  correlations  of 
the  Reid  soft  core  interaction  were  carefully  treated  and  the  resulting  sp  properties  were 
studied  in  detail.  This  study  is  also  relevant  for  finite  nuclei  since  short-range 
correlations  should  be  similar  in  finite  and  infinite  systems  for  the  same  bare  interaction. 
It  will  be  shown  that  the  concept  of  healing  as  discussed  in  the  literature  has  no  strict 
validity  since  the  particles  in  the  medium  do  scatter  at  all  energies.  In  sec.  2  of  this 
contribution  the  results  of  these  calculations  which  employ  the  Green’s  function  method 
will  be  reviewed  and  its  implications  for  the  nuclear  shell  model  indicated.  Special 
attention  will  be  paid  to  the  nucleon  spectral  function  and  the  relevance  of  the  quasi¬ 
particle  concept.  In  sec.  3  the  modifications  of  these  results  due  to  the  finiteness  of  the 
system  will  be  discussed.  Results  for  calculations  of  the  spectral  function  for  '60  will 
illustrate  the  usefulness  of  the  quasi-particle  idea  for  finite  nuclei  and  its  relation  to  the 
success  of  the  simple  shell  model. 


2.  SHORT-RANGE  CORRELATIONS  AND  SINGLE-PARTICLE  PROPERTIES 

In  ref.  6  the  original  idea  of  Brueckner  how  to  treat  the  influence  of  short-range 
correlations,  is  embedded  in  a  Self-Consistent  Green’s  Function  (SCGF)  approach.  The 
ladder  equation  is  solved  with  inclusion  of  both  particle-particle  (pp)  as  well  as  hole- 
hole  (hh)  propagation.  The  connection  between  this  interaction  and  the  resulting  sp 
properties  is  made  by  calculating  the  resulting  self-energy.  This  self-energy  can  then  be 
used  to  calculate  the  sp  propagator  or  the  particle  and  hole  spectral  functions  in  nuclear 
matter.  In  principle,  this  requires  a  self-consistent  solution  in  which  the  full  dressing  of 
the  sp  propagator  is  included  in  the  determination  of  the  ladder  correlations.  At  present 
this  self-consistency  has  been  established  in  ref.  6  for  the  quasi-particle  energy  which 
corresponds  to  the  inclusion  of  only  the  real  on-shell  part  of  the  self-energy  in  the  self- 
consistency  procedure. 

k2 

e(k)  =  +  Re  I(k,e(k))  .  (1) 

2m 

When  this  average  self-consistency  is  achieved  the  complete  energy  and  momentum 
dependence  of  the  self-energy  can  be  studied.  This  requires  knowledge  of  the  imaginary 


pari  of  the  self-energy  above  but  also  below  the  Fermi  energy.  It  is  this  last 
contribution  which  has  never  been  considered  in  the  conventional  treatment  using  the 
Bethe-Goldstone  equation  which  only  considers  pp  propagation  in  the  ladder  equation. 
Including  hh  propagation  one  naturally  treats  the  coupling  of  sp  degrees  of  freedom  to 
2hlp  and  more  complicated  states  and  therefore  one  obtains  an  imaginary  contribution 
to  the  self-energy  below  eF.  A  Green’s  function  procedure  automatically  treats  hh 
propagation  on  the  same  footing  as  pp  propagation.  Recently,  this  hh  propagation  (to  all 
orders)  has  also  been  shown7  to  be  the  crucial  link  in  connecting  the  coupled  ladder 
self-energy  problem  with  the  description  of  pairing,  going  even  beyond  the  conventional 
BCS  description.8 

Including  hh  propagation  in  the  ladder  equation  implies  that  the  resulting 
interaction  is  complex  not  only  for  energies  above  2eF  but  also  below.  This  conflicts 
with  the  notion  of  healing  which  has  been  related  to  the  G  matrix  being  purely  real  for 
energies  below  2ep  in  the  case  of  the  Bethe-Goldstone  equation.5  Obviously,  collisions 
do  take  place  in  the  correlated  system  at  any  energy.  This  statement  must  be  regarded 
as  general  for  a  normal  Fermi  liquid.  The  validity  of  shell  model  concepts  in  finite 
nuclei  can  therefore  not  be  linked  to  this  healing  property  of  the  relative  wave  function. 

Calculations  have  been  performed  for  the  central  part  of  the  3Si  -  3Dj  channel  of 
the  Reid  soft  core  potential.9  This  avoids  the  strong  pairing  instability10  in  the.deuteron 
channel  and  allows  a  careful  study  of  the  influence  of  short-range  correlations  on  the  sp 
self-energy.6  The  self-energy  of  a  nucleon  contains  all  the  information  necessary  to 
obtain  occupation  probabilities,  quasi-panicle  strength  and  broadening  features  which 
alternatively  can  be  visualized  in  terms  of  hole  and  particle  spectral  functions.  Results 
for  the  hole  spectral  function  are  given  in  Fig.  1  for  k  =  0.828,  1.258,  2.220,  and 
4.376  fm-1  at  a  density  corresponding  to  kp  =  .1.4  fm-1.  Results  for  the  particle  spectral 
function  are  given  in  Fig.  2  for  the  same  momenta.  Clearly  one  observes  the  increasing 
validity  of  the  quasi-particle  concept  when  k  approaches  kp  since  the  peak  of  the 
spectral  function  becomes  sharper  the  closer  the  momentum  is  to  kp.  This  holds  for 
momenta  below  kF  (Fig.  1)  as  well  as  for  momenta  above  kp  (Fig.  2).  One  observes 
therefore  in  this  actual  calculation  which  uses  a  strongly  repulsive  NN  interaction  at 
short  distances,  that  the  notion  of  quasi-particles  is  valid  in  exactly  the  same  sense  as  it 
has  been  used  by  Landau  in  his  description  of  Fermi  liquids.11  In  an  infinite  system  this 
concept  is  strictly  valid  only  for  the  Fermi  momentum.  For  other  momenta  one 
observes  therefore  that  the  further  the  momentum  value  is  removed  from  kp  the  more 
the  quasi-particle  peak  is  broadened. 

For  momenta  close  to  kp  most  of  the  hole  strength  is  under  the  quasi-particle  peak 
when  k  is  smaller  than  kp.  For  k  =  1.258  fm-1  83%  of  the  hole  strength  is  under  the 
peak  whereas  3%  is  found  at  lower  energies.  The  occupation  of  this  momentum  in  the 
correlated  ground  state  is  obtained  by  integrating  the  hole  spectral  function  up  to  £p  and 
therefore  corresponds  to  86%.  The  remaining  strength  is  moved  beyond  the  Fermi 
energy  and  belongs  to  the  particle  spectral  function  shown  in  Fig.  2.  The  remaining 
14%  of  the  strength  is  smeared  out  over  a  very  large  energy  domain  which  is  intimately 
connected  with  the  character  of  the  interaction  which  couples  sp  degrees  of  freedom  to 
very  high-lying  2plh  etc.  states.  This  removal  of  sp  strength  to  high  energy  is  observed 
for  all  momenta  by  the  same  amount  (see  Fig.  2)  emphasizing  the  short-range  nature  of 
the  effect.  The  jump  in  the  momentum  distribution  at  kp  is  equal  to  the  strength  of  the 
quasi-particle  pole  at  kp  and  a  value  of  0.83  is  obtained  here.  The  inclusion  of  tensor 
correlations  is  expected  to  decrease  this  value  further  by  a  few  percent.12 
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Fig.  1 .  Hole  spectral  functions  at  kF  =  1.4  fm-1  for  k  =  0.828  1 .258,  2.220,  and  4.376  fm"1. 


Fig.  2.  Panicle  spectral  functions  at  kF  =  1.4  fm-1  for  k  =  0.828  (solid), 

1.258  (short  dashed),  2.220  (long  dashed),  and  4.376  fm-1  (dot  dashed). 


These  and  other  results  arc  more  fully  discussed  in  ief.  6  together  with  the 
calculational  details.  One  can  conclude  from  these  results  that  even  though  there  is 
strong  scattering  in  the  medium,  the  concept,  of  a  quasi-particle  picture  is  still  extremely 
useful  in  exactly  the  same  way  as  it  has  been  used  in  Landau’s  Fermi  liquid  'theory11  for 
3He.  For  a  normal  Fermi  liquid  this  validity  of  the  quasi-particle  concept  is  almost 
exclusively  a  phase  space  effect  related  to  the  presence  of  2plh  and  2hlp  states  around 
the  Fermi  surface  which  can  mix  with  sp  states.  The  strength  of  the  quasi-particle  pole 
is  then  a  good  measure  of  the  correlations,  the  closer  to  1  the  more  the  system  can  be 
described  in  terms  of  a  mean  field  picture. 

The  role  of  short-range  correlations  should  be  similar  in  finite  nuclei  and  infinite 
nuclear  matter  at  the  corresponding  density.  These  results  have  therefore  direct  bearing 
on  the  situation  in  finite  nuclei  and  depletion  effects  due  to  short-range  correlations  are 
inevitable.13  From  the  present  results  one  can  also  infer  that  the  missing  sp  strength  is 
not  just  “around  the  comer”;  it  is  removed  to  very  high  energies  as  a  consequence  of 
the  short-range  correlations.  This  effect  of  the  short-range  correlations  and  its 
consequences  have  not  been  properly  considered  in  Brueckner  type  calculations. 


3.  FINITE  NUCLEI  AND  IMPLICATIONS  FOR  THE  SHELL  MODEL 

Presently,  the  breakdown  of  the  mean  field  picture  is  an  experimental  reality. 
Combined  analysis  of  (e,e'p)  reactions14  together  with  absolute  charge  density 
distribution  measurements15  indicates  that  the  occupation  of  the  3si/2  proton  shell  in 
208Pb  is  80  ±  10%.  This  should  be  compared  to  the  theoretical  results  of  refs.  13,16  and 
17.  Direct  analysis  of  (e.e'p)  cross  sections  with  shell  model  momentum  distributions 
suggests  that  the  discrete  low  energy  transitions  carry  about  50%  of  the  hole 
strength.18,19  This  would  imply  that  30%  of  the  hole  strength  is  at  even  higher  energies 
in  the  A-l  nucleus  and  cannot  be  detected.  The  analysis  of  ref.  16  for  the  neutron  states 
in  208Pb  indicates  that  about  73%  of  the  hole  strength  is  found  at  low  energy  for  such 
states  with  occupation  numbers  corresponding  to  0.85.  Assuming  similar  results  for 
protons  this  indicates  a  20%  discrepancy  between  this  theoretical  result  and  the 
experimental  analysis 

To  address  this  question  it  is  profitable  to  investigate  the  influence  of  low-lying 
excitations  on  the  redistribution  of  sp  strength  directly  in  finite  nuclei.  One  should  keep 
in  mind  that  the  nuclear  matter  results  discussed  above  indicate  that  somewhat  more 
than  10%  of  the  sp  strength  is  at  high  energy.  This  result  is  due  to  a  basic  asymmetry 
between  the  influence  of  pp  and  hh  contributions  and  should  be  taken  into  account  in 
analyses  like  the  one  in  ref.  16.  The  spectral  functions  in  finite  nuclei  can  be  obtained 
directly  by  solving  the  Dyson  equation  in  a  suitable  basis  for  the  relevant  nucleus. 
Results  are  reported  here  for  160.  First  a  Hartree-Fock  (HF)  basis  is  generated  by 
solving  the  HF  equations  using  the  G  matrix  interaction  G/Jm  which  is  a  good  local  and 
static  approximation  of  a  160  G  matrix 20  As  in  any  HF  calculation  a  mean  field  is 
generated  in  which  one  can  define  fully  occupied  and  completely  empty  sp  states  for 
this  nucleus.  The  G  matrix  that  is  used  here  does  not  include  any  shell  structure 
information  since  it  is  calculated  originally  for  nuclear  matter  with  a  standard  gap  in  the 
sp  spectrum.21 

To  calculate  the  influence  of  low  energy  excitations  on  the  distribution  of  sp 
strength  the  second  order  self-energy  contribution  calculated  with  this  interaction  is  used 
to  solve  the  Dyson  equation. 
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The  sp  quantum  numbers  a,p,  etc.  refer  to  the  HF  basis.  The  sums  in  Eq.  (2)  are 
restricted  to  HF  states  within  50  MeV  of  the  Fermi  energy  of  the  HF  calculation.  This 
ensures  that  indeed  the  influena  of  low-lying  2plh  and  2hlp  states  on  the  redistribution 
of  sp  strength  in  a  finite  system  is  studied.  Contributions  from  higher  shells  will 
eventually  also  lead  to  some  form  of  double  counting  since  these  states  have  already 
been  considered  in  the  construction  of  the  residual  interaction.  The  presence  of  an 
explicit  energy  dependence  in  this  second  order  self-energy  contribution  is  responsible 
for  the  fragmentation  of  sp  strength  both  above  the  Fermi  energy  by  coupling  to  2plh 
states  as  well  as  below  the  Fermi  energy  by  coupling  to  2hlp  states.  This  can  be 
clearly  seen  by  inspection  of  the  explicit  second  order  self-energy  contribution 
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The  solution  to  Eq.  (2)  has  the  form 
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The  poles  in  this  equation  give  the  energies  of  the  states  in  the  A  ±  1  systems  which 
can  couple  to  the  ground  state  of  the  A-body  system  through  the  creation  and 
annihilation  operators  a  and  aT,  respectively.  The  residues  at  the  poles  for  the  diagonal 
sp  propagator  correspond  to  the  spectroscopic  factors  for  the  corresponding  transition 
from  the  ground  state  in  the  A-body  system  to  the  state  n(m)  in  the  A+l(A-l)-body 
system.  The  energies  and  residues  in  Eq.  (4)  are  obtained  by  solving  an  energy 
dependent  eigenvalue  problem  which  can  be  derived  from  Eq.  (2) 22 

In  Fig.  3  results  for  the  sp  strength  distribution  are  shown  for  the  relevant  si/2 
states  in  160.  The  results  are  smeared  by  distributing  the  strength  of  each  peak  by 
0.5  MeV.  In  the  HF  picture,  the  lowest  sl/2  state  has  full  occupation,  i.e.,  a 
spectroscopic  factor  of  1.0  to  one  state  in  the  A-l  nucleus.  This  result  is  also  typical 
for  Brueckner-Hartree-Fock  calculations  in  which  the  energy  dependence  of  the  G 
matrix  is  treated  explicitly.  When  no  hh  contributions  in  the  ladder  equation  are 
included  the  resulting  self-energy  is  real  below  the  Fermi  energy  and  therefore  leads  to 
HF  like  hole  states.  The  next  si/2  state  is  completely  empty  and  has  a  spectroscopic 
factor  of  1.0  to  one  state  in  the  A+l  system.  Including  dynamic  correlations  in  this 
second  order  self-energy  calculation  still  leads  to  a  concentration  of  the  2st/2  strength  in 
one  peak  as  shown  in  Fig.  3.  In  contrast  the  lsi/2  strength  (represented  by  all  the  other 
peaks  in  the  figure)  is  completely  fragmented  with  most  of  its  strength  below  the  Fermi 
energy  leading  in  this  case  to  an  occupation  of  96%.  The  remaining  4%  is  located  at 
higher  energy  shown  by  the  smaller  peaks  to  the  right  of  the  2si/2  peak. 
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Fig.  3.  Distribution  of  sl/2  strength  in  160  as  a  function  of  energy.  The  2sl/2  strength 
is  located  mainly  in  one  peak  reduced  here  by  a  factor  of  1/100.  The  rest  of 
the  2sl/2  strength  is  not  shown.  All  other  peaks  correspond  to  strongly  fragmented 
lsl/2  strength.  Note  that  a  small  portion  of  lsl/2  strength  is  also  found  above 
die  Fermi  energy. 

This  calculation  therefore  has  the  same  typical  features  that  one  encounters  in 
nuclear  matter.  For  sp  states  which  are  close  to  the  Fermi  energy  the  dynamic  coupling 
in  the  self-energy  results  in  one  strong  peak  which  carries  less  than  but  of  the  order  of 
100%  of  the  strength.  In  a  finite  nucleus  like  160  this  is  particularly  obvious  since  there 
are  no  2plh  or  2hlp  states  in  the  immediate  "vicinity  of  the  2si/2  energy  and  the  strength 
remains  mainly  concentrated  in  one  state  since  more  complicated  states  are  too  far  away 
in  energy.  This  feature  emphasizes  the  important  role  of  the  finiteness  of  the  system.  It 
leads  automatically  to  a  finite  ph  energy  gap  which  results  in  a  larger  window  of 
validity  of  the  quasi-particle  approximation  in  finite  nuclei.  For  the  lsl/2  state  one 
observes  a  strong  fragmentation  since  its  energy  (in  HF)  is  located  in  an  energy  interval 
where  (in  this  case)  many  2hlp  are  located  and  it  becomes  relatively  easy  to  mix  with 
such  more  complicated  states.  The  same  result  is  observed  in  infinite  systems  for 
momenta  far  from  the  Fermi  momentum.  Experimental  results  on  finite  nuclei  by  means 
of  the  (e,e'p)  reaction  show  a  smooth  lsi/2  strength  distribution  for  proton  knock-out  of 
160.23  Recent  (p,2p)  and  (p,pn)  experiments  on  160  show  the  same  feature.24  Clearly 
the  knock-out  of  an  lsi/2  particle  is  possible  for  a  broad  range  of  energies. 

The  calculation  gives  similar  results  as  discussed  above  for  other  orbitals.  Panicles 
in  shells  closer  to  the  Fermi  energy  like  the  lp3/2  and  the  lpi/2  are  found  in  large 
fractions  in  the  low  energy  stares  r.f  the  A  =  15  system  but  do  not  carry  100%  of  the 
strength.  On  the  other  side  of  the  Fermi  energy  most  of  the  lds/2,  2si/2  (see  Fig.  3)  and 
ld3/2  strength  is  found  in  the  low  energy  states  of  the  A  =  17  system.  Unoccupied 
states  which  in  HF  are  far  removed  from  the  Fermi  energy  show  a  similar  fragmentation 
as  the  Is  1/2  state  on  the  other  side  of  the  Fermi  energy. 

These  results  are  in  sharp  contrast  with  the  mean  field  picture  of  nuclei.  In  an 
atom  e.g.,  where  the  Hartree-Fock  approximation  works  very  well,  all  hole  states  have  a 
well  defined  energy  and  carry  100%  of  the  strength.  The  situation  in  a  nucleus 


thcrefoie  resembles  much  more  lhat  in  liquid  'He  and  in  nuclear  matter  as  discussed  in 
sec.  2.  Nuclear  single-particle  excitations  as  observed  near  the  Fermi  energy  must 
therefore  be  reinterpreted  as  quasi-particle  or  quasi-hole  excitations  which  carry  a 
considerable  amount  of  sp  strength  Present  experimental  results  put  the  nucleus 
somewhere  in  between  atoms  and  liquid  ?He  as  far  as  the  strength  of  correlations  are 
concerned.  States  which  are  far  away  from  the  Fermi  energy  in  mean  field  are  located 
at  energies  at  which  a  considerable  number  of  2p  1  h  or  2hlp  states  are  found.  The 
strength  of  the  nuclear  interaction  is  such  that  considerable  mixing  between  these 
degrees  of  freedom  occurs  leading  to  the  appearance  of  sp  strength  in  a  broad  energy 
range. 

The  experimentally  observed  smooth  distribution  of  hole  strength18, 19,23,24  further 
emphasizes  the  non-linear  aspects  of  the  problem  of  calculating  the  sp  propagator.  The 
present  calculation  uses  HF  sp  propagators  in  the  calculation  of  the  second  order  self- 
ene  gy  contribution.  However,  also  these  intermediate  2hlp  and  2plh  states  are  not 
good  eigenstates  (see  Eq.  (3))  and  are  broadened  themselves.  This  effect  can  be 
incorporated  by  using  self-consistent  sp  propagators  in  the  calculation  of  the  self-energy 
in  the  same  non-linear  fashion  as  the  HF  problem  is  solved.  This  self-consistent 
formulation  is  necessary  to  explain  the  observed  experimental  results  of  smooth  hole 
strength  distributions.  This  observation,  suggests  that  it  is  not  very  realistic  to  obtain 
information  on  the  widths  of  such  strength  distributions  from  nuclear  matter 
calculations.  The  role  of  the  ph  gap  in  a  finite  system  in  the  fragmentation  of  sp 
strength  should  be  explicitly  treated. 

The  results  obtained  here  for  160  do  r.ot  yet  include  the  effect  of  short-range  and 
tensor  correlations  on  the  sp  propagator  but  instead  emphasize  the  role  of  the  finiteness 
of  the  system  with  its  shell  structure  and  surface.  These  calculations  show  that  up  to 
5%  of  strength  (using  second  order  self-energy)  can  be  found  in  tiny  fragments  below 
the  Fermi  energy  as  a  result  of  the  coupling  of  sp  states  to  low-lying  2h Ip  states. 
Taking  the  global  effect  of  short-range  correlations  (as  discussed  in  sec.  2)  into  account 
and  adding  another  5%  depletion  effect  due  to  tensor  correlations,12,17  one  expects  about 
20%  of  the  sp  strength  to  be  located  above  eF  for  normally  occupied  states.  Still,  most 
of  the  remaining  strength  for  these  valence  hole  states  should  be  found  in  discrete  peaks 
and  not  much  more  than  10%  at  lower  energies  (higher  excitation  energy  in  the  A-l 
system)  according  to  the  present  theoretical  results.  This  is  in  reasonable  agreement 
with  die  results  of  ref.  16,  but  disagrees  with  the  present  analysis  of  the  (e,e'p) 
experiment18,19  which  assigns  about  50%  of  the  hole  strength  to  discrete  states  at  low 
energy  and  from  a  combined  analysis  together  with  charge  density  distributions15  leads 
to  80%  occupation  for  the  3si/2  state  in  208Pb. 

From  the  results  discussed  here  and  in  sec.  2  one  can  conclude  that  the  success  of 
the  conventional  shell  model  is  related  to  the  success  of  the  quasi-particle  picture  for 
nuclear  systems  and  the  finiteness  of  the  nucleus  which  extends  the  window  of  the 
validity  of  this  picture.  In  a  Fermi  liquid  this  picture  pertains  strictly  to  the  Fermi 
momentum  and  energy.  Since  particle-hole  excitations  start  already  at  zero  excitation 
energy,  there  is  immediately  this  broadening  effect  when  one  moves  away  from  the 
Fermi  surface.  In  contrast  to  this,  a  finite  particle-hole  excitation  energy  exists  in  nuclei 
since  these  systems  are  confined  in  space  which  leads  to  discrete  single-panicle 
energies.  For  this  reason  there  is  a  broader  validity  of  the  quasi-particle  picture  in 
nuclei  and  the  shell  model  makes  sense  when  it  is  not  pushed  beyond  its  window  of 
■validity.  Finally,  one  should  also  remember  the  difficulty  to  explain  the  experimental 
charge  density  distribution  of  nuclei25  when  a  mean  field  picture  is  employed.  From 
partially  occupied  si/2  states  in  a  nucleus  one  obtains  an  automatic  reduction  of  the 
central  charge  density  which  can  qualitatively  explain  the  experimental  data. 


The  conclusion  one  can  draw  from  the  above  discussion  is  that  a  marriage  ot  the 
idea  of  Brueckner  to  treat  short-range  correlations  and  the  notion  of  quasi-particles 
employed  by  Landau  both  developed  already  in  the  fifties,  leads  to  a  fascinating  picture 
of  the  nucleus  as  a  correlated  many-body  system  and  also  provides  a  basis  for  a 
theoretical  foundation  of  the  nuclear  shell  model  in  the  presence  of  strong  interactions. 
Present  experimental  and  theoretical  evidence  suggests  that  nuclei  are  considerably  more 
correlated  than  atoms  and  the  mean  field  picture  for  nuclei  is  not  valid. 
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1.  INTRODUCTION 

In  recent  years,  mixed  quantum-classical  systems  consisting  of  excess 
electrons  interacting  with  classical  many-body  systems  at  finite  temperatures 
have  been  studied  extensively  with  computer-simulation  techniques1.  The 
simulation  methods  for  these  systems  include  the  path  integral  Monte  Carlo1  or 
molecular  dynamics2  and  dynamical  simulated  annealing3.  The  latter  can  only 
provide  the  ground-state  static  properties  of  the  quantum  particles.  The  path 
integral  approach  has  been  used  successfully  to  calculate  the  equilibrium 
properties,  but  the  study  of  time  correlation  functions4  is  not  reliable  at  long 
times.  However,  the  recently  developed  quantum  molecular  dynamics  method, 
which  deals  directly  with  the  time-dependent  Schrodinger  equation,  contains  all 
the  dynamical  information  for  quantum  particles. 

2.  QUANTUM  MOLECULAR  DYNAMICS  METHOD 

Quantum  molecular  dynar'cs  method5  provides  the  real-time  dynamics  of 
quantum  and  classical  particU  .  n  mixed  systems  at  finite  temperatures  through 
the  numerical  solutions  of  iii^  ;ime-dependent  Schrodinger  equation  for 
quantum  particles  and  Newton’s  equations  of  motion  for  classical  particles.  To 
understand  this  technique,  consider  a  quantum  particle  of  mass  m  described  by 

;he  wave  function,  \}/(r,t),  interacting  with  a  classical  system  of  N  particles  with 
masses  Mi  and  positions  {R-,(t)}.  For  the  quantum  particle  we  have 


(1) 


.<fy(r,t) 
1  dt 


+  V(r)  v|/(r,t); 


N  __ 

V(r)  =  X  v(f-Ri) 

i=i 


where  v  is  the  interaction  potential  between  the  quantum  and  classical  particles 
and  we  use  atomic  units, H  =  e  =  1.  In  the  Born-Oppenheimer  approximation,  the 
positions  of  the  classical  particles  evolve  from  Newton’s  equations  of  motion: 

MiR,  =  -ViU«Ri))  -  V,  f  dr  |y(r,t)|2  v(r-Ri)  (2) 

where  U  is  the  potential  energy  for  the  classical  particles.  Choosing  a  small 
time  step  At,  the  solution  of  Eq.  (1)  can  be  written  as6 

\[i(r,t+At)  =  eiAt?2/4m  e-i4tYeiAtv2/4m¥(f.o  +  o[(&t)3]  (3) 

Equation  (3)  is  solved  with  fast  Fourier  transform  (FFT)  techniques6.  First,  note 

e i A t V2/4m¥(?>t }  =  ¥(kit )  e‘iAt k2/4me 1  kf  (4) 

k 


which  means  multiplying  the  Fourier  transform  of  i.  e.,  \y(k,t)  with  exp(-iAt 
k2/4m),  followed  by  an  inverse  FFT.  Next,  the  outcome  of  Eq.  (4)  is  multiplied  by 
exp(  -  iAt  V).  Finally,  the  FFT  of  the  resultant  of  the  last  step  is  multiplied  with 

exp(-iAt  k2/4m)  and  then  the  inverse  FFT  is  taken.  These  three  steps  are 
repeated  to  obtain  the  time  evolution  of  the  wave  function.  The  classical 
equations  of  motion  can  be  integrated  numerically  with  one  of  several  available 
algorithms7. 

For  systems  with  broken  symmetry  due  to  the  presence  of  surfaces  or  an 
external  electric  field,  the  use  of  periodic  boundary  conditions  in  the  broken- 
symmetry  direction  is  inappropriate  and  so  the  first  and  third  steps  in  the  time¬ 
stepping  operation  in  Eq.  (3)  cannot  be  executed  with  the  FFT  method.  For 
these  physical  situations  a  new  QMD  algorithm8  has  been  developed. 

3.  AN  EXCESS  ELECTRON  IN  A  DENSE  HELIUM  GAS 

An  injected  electron  in  a  sufficiently  dense  helium  gas  tends  to  localize  in  the 
form  of  a  bubble  as  a  result  of  the  strongly  repulsive  electron-helium  interaction9 
at  short  and  intermediate  distances.  Recently  we  investigated  the  electron 
bubble  formation  using  the  QMD  scheme10.  For  the  electron  in  helium  problem, 
the  helium  particles  interact  with  each  other  via  a  two-body  Lennard-Jones 


potential  with  parameters  s=  10.22K  and  c=  2.576  A,  and  with  the  electron  via  a 
pseudopotential9. 

The  QMD  simulations10  were  performed  at  77K  for  reduced  helium  densities 
n  =  po3  =  0.1 , 0.1 7,  and  0.25  which  correspond  to  p  =  0.61 , 1 .0,  and  1 .46  x  1 022 
cm'3.  Systems  with  n  =  0.1  and  0.17  contained  512  helium  atoms  while  the 
simulations  at  the  highest  helium  density  were  carried  out  with  64  and  140 
particles.  These  gave  similar  results  though  t‘:a  sm  life*'  svstem  showed 
evidence  of  finite-size  effects.  The  electron  wave  packet  was  propagated  with 

fast  Fourier  transforms  on  323  grid  points  and  with  time  step,  At,  in  the  range  of 
0.2  -  0.5  a.u.  The  total  energy  was  conserved  to  better  than  0.1%  over  106  time 
steps.  Some  simulations  with  643  grid  points  gave  identical  results.  Classical 
molecular  dynamics  for  helium  atoms  was  performed  in  the  canonical 
ensemble11. 

4.  RESULTS 

Figure  1  shows  the  electron-helium  radial  distribution  function,  G(r), 
measured  relative  to  the  center-of-mass  of  the  electron  wave  packet  for  n  =  0.25. 
At  this  high  helium  density,  there  are  no  helium  atoms  up  to  r  ~  12  a.u.  and  this 
excluded-volume  effect  is  an  indication  of  an  electron  bubble.  The  excluded 

volume  is  also  present  at  n  =  0.17.  The  shape  parameter10  for  the  electron 
bubble  indicates  that  the  bubble  is  spherical.  From  the  excluded-volume  region 
in  G(r)  we  estimate  the  size  of  the  bubble  to  be  ~  12  a.u.  At  the  lowest  density,  n 
=  0.1,  G(r)  does  not  display  the  excluded-volume  behavior  as  the  wave  packet 
de-localizes. 


Figure  1 :  Electron-helium  radial  distribution  function,  G(r),  for  n  =  0.25. 


Next,  let  us  examine  the  time  variation  of  the  extent  of  the  self-trapped 
electron.  At  n  =  0.17,  the  participation  ratio  remains  small  and  almost  constant 
over  the  entire  simulation,  indicating  that  the  electron  is  localized.  The  behavior 
is  quite  similar  at  n  =  0.25,  but  very  different  at  n  =  0.10.  Figure  2  shows  the 

participation  ratio,  p(t)  =  (njdr  |  \j/(r,t)  I4)'1,  normalized  to  the  volume  of  the  MD 


cell  for  n  =  0.1 .  This  simulation  was  started  by  expanding  the  length  scale  of  the 
final  configuration  of  the  n  =  0.17  run  to  a  value  corresponding  to  n  =  0.1  and 
projecting  out  the  new  electron  ground  state,  which  corresponds  to  an  almost 
spherical  bubble.  Therefore  the  participatin  ratio  starts 


t(ps) 


Figure  2:  Participation  ratio  for  n  =  0.1  as  a  function  of  time. 


out  at  a  value  determined  primarily  by  the  bubble  size  at  n  =  0.17  scaled  by  the 
increase  in  the  cell  size  in  changing  the  system  density  from  n  =  0.17  to  0.1 .  It  is 
observed  that  the  bubble  quickly  expands  to  a  volume  that  is  a  sizeable  fraction 
of  the  cell  volume.  After  4  ps  the  electron  wave  packet  undergoes  large 
expansions  and  contractions  and  finally  after  8  ps  the  wave  packet  occupies 
almost  the  entire  volume  of  the  cell.  This  expansion  indicates  de-localization  of 
the  electron  that  is  very  likely  limited  by  the  size  of  the  system.  The  wave  packet 
attempts  to  find  a  region  of  the  cell  where  the  helium  density  is  less.  Thus,  as 
the  simulation  proceeds,  the  wave  packet  localizes  again  temporarily  because 
of  thermal  fluctuation  in  the  system.  This  partial  localization  and  de-localization 
is  expected  to  continue  indefinitely. 

Figure  3  shows  the  current-current  correlation  function  of  the  electron,  x’’(co), 
for  n  =  0.25  and  0.1  obtained  by  Fourier  transforming  the  time  correlation 

function  <V\j/fVv0  >  and  taking  the  configurational  average  within  the  Franck- 

Condon6  approximation.  For  n  =  0.25,  x”(a),  which  is  simply  related  to  the 
optical  absorption  spectrum,  displays  a  significant  structure  which 
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Figure  3:  Imaginary  part  of  the  current-current  correlation  function  as  a  function 
of  energy  in  electron  volts  for  helium  densities  n  =  0.25  and  0.1.  The  peak  at  n  = 
0.25  is  due  to  intra-bubble  transitions. 


shifts  to  a  lower  energy  as  the  helium  density  is  reduced  to  0.17.  These  peaks 
are  present  throughout  the  simulations.  At  n  =  0.1 ,  when  the  wave  packet  is  de¬ 
localized,  the  structure  disappears  and  only  a  background  is  observed. 

At  T  =  77K  we  have  also  calculated8  the  electron  mobility  at  a  helium  density 
of  1.25  x  1 022  cm'3.  Applying  an  electric  field  of  2.6  x  105  volt/cm  in  the  x 
direction,  it  is  observed  that  the  electron  drifts  with  a  velocity  of  (2.2  ±  0.5)  x  10  4 
cm/s.  The  applied  field  is  small  since  it  does  not  produce  any  noticeable 
changes  in  G(r).  From  this  simulation  we  obtain  the  electron  mobility  to  be  (0.08 
±  0.02)  cm2/volt-s  which  is  in  good  agreement  with  the  extrapolated 
experimental  value12  (0.1  cm2/volt-s). 

5.  CONCLUSION 

In  this  paper  we  have  presented  the  QMD  simulation  technique  for  studying 
the  dynamics  of  quantum  particles  in  mixed  systems  at  finite  temperatures.  This 
technique  is  applied  to  simulate  the  behavior  of  an  excess  electron  in  helium 
gas  at  77K.  It  is  found  that  the  electron  localizes  in  the  form  of  a  nearly  spherical 
bubble  of  approximately  7  A  in  radius  above  a  critical  density  of  0.6  x  1022  cm*3. 
The  bubble  possesses  quasi  bound  excited  states  and  intra-bubble  dipole 
transitions  between  these  states  give  rise  to  a  pronounced  structure  in  the 
optical  absorption  spectrum.  Below  the  critical  density  the  electron  percolates 
through  i’ne  helium  gas  and  displays  a  featureless  excitation  spectrum. 
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This  contribution  reports  on  recent  progress  in  the  investigation  of  the  structure  and  excitations 
of  quantum  liquids  adsorbed  to  a  solid  surface,  the  states  of  impurity  atoms  in  such  films,  and  the 
interaction  between  them.  The  subtitle  of  my  talk  is  UA  Generic  Many-Body  Problem  >\  it  indicates  that 
I  will,  to  some  extent,  emphasize  the  methodological  aspect  of  the  problem. 

We  start  the  microscopic  description  of  a  manv-body  system  of  N  particles  in  a  given  volume  fl  in 
an  external  potential  with  an  empirical  Hamiltonian 


H  = 


rjl)- 


(1) 


To  be  specific,  we  consider  a  system  of  4IIe  atoms  at  zero  temperature.  Ideally,  one  would  like  to 
solve  the  Schrodinger  equation  for  this  Hamiltonian.  But,  even  with  substantial  computational  resources, 
this  is  only  possible  only  for  a  few  systems  with  simple  geometries  and  simple  interactions.  Therefore, 
we  must  resort  to  approximations.  However,  obvious  physical  consequences  of  the  exact  problem  are  not 
necessarily  satisfied  by  an  approximate  theory. 

A  microscopic  theory  must  be  able  to  deal  with  the  short-ranged  repulsion  between  particles  and 
include  the  long-ranged  correlations  correctly.  Considering  N  particles  in  a  given  volume  the  theory 
must  also  be  able  to  decide  whether  particles  fill  the  given  volume  uniform!},  or  only  a  part  of  it.  There 
may  be  a  regime  of  average  densities  N/Q  where  both  states  are  possible.  This  is  where  the  pressure 
of  the  uniform  phase  is  negative,  but  the  compressibility  is  positive.  However,  if  the  uniform  phase  is 
diluted  to  the  density  point  where  the  compressibility  becomes  negative,  the  theory  for  the  uniform  phase 
should  cease  to  have  solutions.  In  other  words,  the  theory  must  describe  the  response  of  the  system  to 
long  wavelength  excitations  correctly. 

The  correct  treatment  of  the  bulk  phase  is  prerequisite  to  determining  whether  a  system  of  particles 
will  or  will  no!  bo  adsorbed  to  a  substrate.  The  requirement  of  the  correct  inclusion  of  short-  and 


v  long-ranged  correlations  is  needed  when  wc  consider  the  structure  of  the  adsorbed  atoms.  Hcuristicaily, 
one  would  think  of  the  4  He  atoms  as  hard  spheres.  These  spheres  would  first  form  a  layer  closest  to  the 
surface,  then,  as  their  number  is  increased,  a  second  layer,  and  so  on.  Eventually,  the  zero-point  motion 
of  the  atoms  will  win  over  the  attraction  to  the  substrate,  and  far  from  the  substrate  the  4  lie  particles 
will  behave  as  if  they  are  in  the  bulk.  Next,  we  add  some  attraction  to  the  ‘‘hard  spheres’',  which  is 
strong  enough  to  form  a  many-body  bound  state.  The  particles  will  wish  to  stick  together  even  when  the 
attractive  substrate  is  not  there.  In  the  presence  of  the  attractive  substrate,  it  is  not  clear  whether  the 
plane  geometry  (i.e.  a  geometry  that  is  translationally  invariant  in  the  directions  parallel  to  the  surface) 
is  energetically  favourable.  The  question  is  to  what  extent  the  liquid  will  “wet”  the  surface.  Formally, 
we  again  encounter  the  problem  of  stability  against  long-wavelength  excitations. 

Summarizing  the  above  qualitative  arguments,  we  find  that,  in  order  to  adequately  describe  the 
physics  discussed  above,  the  main -body  theory  to  be  used  must  at  least  qualitatively  implement  the 
following  features: 

•  The  theory  must  describe  short-ranged  correlations  in  order  to  account  for  a  layer-structure  of  the 
adsorbed  liquid. 

•  The  theory  must  include  long-ranged  correlations  in  order  to  give  us  the  right  geometry. 

•  Finally,  our  theoretical  description  of  the  system  should  be  flexible  enough  to  take  any  symmetry 
breaking  fully  into  account. 

These  physical  requirements  translate  into  the  simple  dictum  that  a  theory  should  be  internally 
consistent  such  that  it  has  only  stable  solutions.  Phrased  in  the  formal  language  of  many-body  theory, 
our  theory  must  contain  a  self-consistent  description  of  short-  and  long  ranged  correlations,  which  is 
accomplished  by  a  self-consistent  summation  of  ring-  and  ladder  diagrams.  For  simpler  systems  like  the 
bulk  liquid  one  may  get  away  with  simpler  theories  for  a  while,  but  u  '!/  the  “crimes”  will  catch  up 
and  lead  to  unphvsical  predictions  like  a  negative  compressibility.  Even  if  numerical  comparison  with 
experiments  look  good,  it  is  clear  that  something  is  fundamentally  wrong  with  such  a  theory. 

Of  course,  everything  w*e  do  should  eventually  be  compared  with  experiments.  Unfortunately, 
quantities  that  can  be  calculated  easily  are  often  hard  to  measure  and  vice  versa.  From  the  wav  the 
experiments  are  done,  ground-state  energy  measurements  are  difficult.  Neutron  scattering  experiments 
probing  the  liquid  surface  are  difficult  since  ‘he  penetration  depth  of  the  neuton  is  about  100  A,  whereas 
the  surface  width  of  4IIe  is  of  the  order  of  10  A.  ir  c  ue  is  interested  in  exploring  the  layer  structure, 
one  has  to  look  for  experiments  that  are  sensitive  to  the  surface  structure.  Therefore,  it  may  also  be 
necessary  to  extend  the  theory  in  order  to  address  the  quantities  that  can  be  measured  more  easily. 

Going  back  to  the  formal  problem  of  a  self-consistent  summation  of  ring-  and  ladder  diagrams,  it 
has  been  known  for  some  time  that  this  goal  is  accomplished  by  both  the  optimized  hypernctted-chain1 
(HNC/EL)  and  the  parquet-diagram2  theory,  which  ultimately  lead  to  the  same  equations  to  be  solved3. 
Here  we  adopt  the  IINC/EL  version  of  the  theory  since  it  is  more  widely  developed  and  can  more  easily 
tested  by  variational  Monte-Carle  calculations.  One  starts  with  a  variational  ansatz  for  the  ground-state 
wave  function  of  the  N  particle  system  of  4lle  atoms  with  coordinates  rj, 


(2) 


*o(ri,...lrA')  =  exp-  <  ui(r,)  +  ti2(ri,ry) 


l<i<N 


!<•<;<  A 


The  one-body  function  tii(r)  describes  the  spatial  structure  of  the  system,  and  the  two-body  function 
W2(r»,r;)  describes  the  short-  and  long-ranged  correlations  between  pairs  of  particles.  These  functions 
are  determined  by  minimization  of  the  ground-state  energy-expectation  value  1 Eq. 
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=  0. 


(3) 


In  many  cases  the  calculation  of  the  variational  energy  expectation  value  Eq  cannot  be  performed 
exactly.  Using  approzimaie  energy  functionals  may  of  course,  render  the  Euler  equations  (3)  meaningless 
It  is  therefore  important  to  use  an  energy  functional  that  has  meaningful  stationary  points  under  the 
variational  problems  (3).  It  has  been  known  for  some  time  that  the  IINC  hierarchy  of  approximations  has 
meaningful  variational  minima.  To  show  that  the  HNC/EL  equations  also  provide  a  correct  description 
of  the  short-  and  the  long  ranged  correlations,  we  give  the  HNC/EL  equations  for  the  homogeneous  Bose 
liquid4.  These  equations  are  most  conveniently  formulated  in  terms  of  the  pair  correlation  function  <?(r), 
and  the  static  structure  function  S(k)  =  1  4-  p /  d3r[g(r)  -  l]e*k  r: 

=  *  (4) 

with 

tWr)  =  +  L|VV^)|=  -f  (<?(r)  -  1  }wj{r).  (5) 

The  “induced  interaction5”  tu/(r)  is  most  conveniently  represented  in  Fourier-space: 

*(*)“ («) 


(The  three-dimensional  Fourier-transform  is  denoted  by  a  tilde).  It  is  worth  noting  here  that 
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Vpfc(|r-r'|)  = 


6p\(r)6pi(r,y 


(?) 


where  the  variational  derivative  is  taken  for  constant  U2(r,r/).  In  Eq.  (4)  we  recover  the  RPA  expression 
for  the  static  form  factor.  The  HNC/EL  theory  supplements  the  RPA  with  a  microscopic  theory  of  the 
particle-  hole  interaction 

An  alternative  way  to  formulate  the  Euler-Lagrange  equation  is5 


m 


v{r)  -f  wj{r) 


s/g(r)  =  0. 


(8) 


Eq.  (8)  is  the  Boson- version  of  the  Bethe-Goldstonc  equation,  in  which  the  bare  interaction  has  been 
supplemented  by  the  “induced  interaction”  wj(r). 

Eqs.  (4)  and  (8)  are  algebraically  equivalent,  they  merely  suggest  different  iteration  paths  for 
the  determination  of  the  optimal  g(r).  The  self-consistent  solution  of  Eqs.  (4-6)  takes  three  Fourier 
transforms  per  iteration,  compared  with  two  rouiier  transforms  per  iteration  for  the  solution  of  the 


UNO  equations  with  a  parametrized  JasUou  function  U2(r).  Thus,  the  efficiency  of  tlic  optimization 
algorithm  makes  the  use  of  a  parametrized  trial  function  completely  pointless.  The  same  is  true  for 
three-body  correlations  and  for  the  Fermion  version  of  the  IINC/EL  approach. 

Returning  to  the  physics,  wc  see  that  the  HNC/EL  theory  includes  both  the  ring  and  the  ladder 
diagrams  exactly,  and  the  mixed  diagrams  approximately2,  in  particular,  Eq.  (4)  shows  that  the 
theory  has,  as  required,  no  uniform  solution  with  negative  compressibility,  c.f.  Eq.  (7).  The  appropriate 
symmetry  breaking  (a  droplet  or  a  plane  surface  geometry,  for  example)  is  put  into  the  wave  function  (1) 
by  including  a  one-body  function  ui(r)  and  breaking  the  translational  invariance  of  «2(rnr;).  If  system 
is,  within  the  assumed  geometry,  stable  against  small  perturbations,  then  the  theory  has  solutions.  We 
conclude  that  the  HNC/EL  theory  satisfies  all  the  minimum  requirements  that  are  needed  to  deal  with 
the  problem  of  a  quantum  liquid  surface. 

The  simplest  symmetry  breaking  that  can  be  treated  with  reasonable  computational  effort  is  the 
plan;  surface  or  the  spherical  geometry.  The  IINC/EL  equations  for  the  inhomogeneous  system  are 
straightforward  generalizations  of  Eqs.  (4-G),  the  equations  and  the  algorithm  for  their  iterative  numer¬ 
ical  solution  has  been  given  in  Ref.  6. 

The  physical  system  considered  here  is  a  number  of  helium  atoms  interacting  via  the  Aziz  potential'. 
The  atoms  are  adsorbed  to  a  substrate  which  is  described  by  an  external  field  U9Ub(z)‘  A  simple  form  for 
Usu i(-)  is  the  potential  obtained  by  averaging  Lennard-Jones  interactions  between  helium  and  substrate 
atoms  over  a  half  space8.  One  obtains 


Given  the  substrate  potential,  the  two-bodv  interaction,  and  the  surface  coverage 

n  =  J  dzpi(z),  (10) 

the  physical  proolem  is  completely  defined.  In  our  geometry,  all  two-body  quantities  depend  on  the 
distances  zt  of  both  particles  from  the  substrate,  and  their  separation  r( j  parallel  to  the  surface.  The 
Euler  equations  for  the  4 He  background  were  solved  numerically  in  the  HNC  approximation,  specifying 
only  the  surface  coverage  n.  Some  typical  examples  of  density  profiles  are  shown  in  Fig.  1.  Most 
remarkable  are  the  stong  density-oscillations  of  the  4IIe-background,  which  are  due  to  the  geometric 
core-exclusion  between  the  4  He  particles. 

It  was  already  pointed  out  that  a  direct  measurement  of  the  surface  profile  by  neutron  scattering 
experiments  is  very  difficult.  It  also  turns  out  that  the  energy  per  particle  does  not  depend  strongly  on 
the  layer  structure  of  the  system9,  and  the  chemical  potential  has  only  a  weak  modulation,  c.f.  Fig.  2. 

One  might  expect  that  the  sound  velocity  exhibits  a  stronger  dependence  on  the  layer  structure  of 
the  surlace,  but  the  analysis  of  our  theoretical  calculations  and  the  experiments  shows  that  the  third 
sound  velocity  c$  (or  the  derived  quantity  mc§)  depends  strongly  on  the  underlaying  substrate  potential 
(Fig.  3)  which  makes  conclusions  on  the  structure  of  the  film  very  difficult. 

So  far,  we  have  seen  that  it  is  difficult  to  probe  the  layer  structure  of  an  adsorbed  4 He  film  directly. 
A  moderately  easy  probe  of  the  surface  structure  is  a  3He  impurity.  Experimental  interest  focusses 
presently  on  the  following  areas: 


Fig.  1:  The  density  pi  (2)  of 
the  background  film  of  Tie  atoms 
(dotted  lines)  and  of  the  3He 
impurity,  p[(z)  (solid  lines)  are 
shown  for  surface  coverages  of  n 
=  0.15,  0*20,  0.25,  0.30,  and  0.35 
^  He  aloms/A3.  Tlie 4 He  densities 
of  all  films  shown  are  virtually 
identical  within  the  first  layer, 
whose  density  maximum  is  about 
0.08  atoms/ A3.  The  3IIe  impu¬ 
rity  density  is  normalized  such 
that  f  dzp[(z)  =  1. 


Fig.  2:  The  chemical  potential 
//  of  A He  atoms  is  shown,  as  a 
function  of  the  surface  coverage 
n,  for  three  different  substrate 

*  i 

potentials  of  different  strength. 
The  upper  curve  corresponds  to  a 
weakly  attractive  glass  substrate, 
the  lowest  one  to  graphite,  and 
the  middle  curve  to  a  potential  of 
average  strength. 


(a)  Measurements  of  the  binding  energy  and  the  specific  heat  of  3He  impurities11,  and 

(b)  Measurement  of  magnetic  properties  of  3IIe  impurities  adsorbed  to  Tie  surfaces12. 

Qualitatively,  one  should  expect  the  following  effects: 

♦  For  a  Him  film  of  Tie  atoms  (one  or  two  layers  thick),  the  substrate  potential  dominates,  and  the  3He 
impurity  atom  will  be  in  the  outermost  layer. 

♦  For  a  thick  film  of  Tie  atoms  (five  or  more  layers),  the  larger  zero-point  motion  of  the  3He  impurity 
will  dominate,  and  the  impurity  atom  will  appear  to  “swim”  on  the  Tie  background. 

♦  If  the  outermost  layer  of  the  4IIe  film  is  full,  it  will  be  hard  for  the  3Ile  impurity  to  move.  The 
“backflow”  effect  and  hence  the  effective  mass  of  the  impurity  atom  will  be  large. 

♦  If  the  outermost  layer  of  the  Tie  film  is  incomplete,  the  backflow,  and  hence  the  effective  mass,  will 


be  smaller. 


Figure  3:  11103  for  tin  medium 

(weakly  oscillating  curve)  and  the 
strong  (strongly  oscillating  curve) 
substrate  potential  are  compared 
with  the  experimental  data  of 
Maynard  and  Chan  (Ref.  10, 
crosses).  See  Ref.  9  for  details 
on  the  analysis  of  the  experiments 
and  the  calculational  procedures. 


The  last  two  points  also  depend  on  how  deep  the  impurity  atom  penetrates  into  the  4He  surface. 

The  above  experimental  situations  require  the  development  of  three  levels  of  theoretical  tools. 

(a)  Single-impurity  theories  aim  at  the  calculation  of  impurity  binding-energy  and  effective  mass.  The 
static  correlations  between  particles  determine  if  and  how  many  bound  states  the  impurity  particle 
can  have  at  the  surface,  and  whether  it  penetrates  to  the  substrate.  The  calculation  of  the  (complex) 
self-energy  of  the  impurity  particle  gives  access  to  the  specific  heat  and  the  mobility. 

(b)  Mixture  films  in  the  dilute  limit  correspond  to  a  system  of  two  static  impurities13.  To  the  extent 
that  the  concentration  of  the  impurities  is  small,  one  can  ignore  all  higher-order  correlations  between 
the  impurity  particles  and  possible  dynamic  effects  due  to  momentum-dependent  correlations. 

(c)  A  theory  for  two-component  systems  is  needed  for  large  “impurity”  concentrations.  This  enables 
us  to  study  the  structure  of  quantum-liquid  interfaces. 

Formally,  impurities  may  be  added  to  the  system  by  extending  the  wave  function  to  include  impurity- 
background  and  impurity-impurit}  correlation  functions  Uj(r;),  U2B(r7,r;),  and  U27(rJ  ,rj),  For  Fermion 
imparities,  the  wave  function  is  multiplied  with  a  Slater  determinant  to  ensure  the  antisymmetry  with 
respect  to  exchange.  To  include  a  momentum  dependence,  one  may  either  add  a  “backflow”  function  to 
the  one-impuritv  correlations,  or  calculate  the  self-energy  in  CBF  perturbation  theory.  The  variational 
wave  function  of  a  system  of  N  4Ile-background  atoms  and  one  3He  impurity  of  momentum  q  is 


/v’ 

w{(r/)  +  jq-r/  +  ^[ 

aq(r7,r,)  is  the  backflow  correlation  function14  describing  the  current  of  4Ile  atoms  flowing  around  the 
moving  impurity.  The  functions  u7(r7),  ti27?(r7,rl),  and  aq(r7,rt)  are  determined  by  minimizing  the 
energy  £q  of  the  system  consisting  of  the  4IIe-background  and  an  impurity  particle  with  momentum  q. 


-u^(r7,r,-)  -f  jaq(r7, r,-)j  }  vf/0(ril . . . ,  r/;).  (11) 


'Mr/,r1>...rAr)  =  exP  j £ 


furl  ty!  \ 


=  0, 


=  0, 


6Eq 

firt.Jr1  1\1 


=  0. 


(12) 


The  total  energy  of  the  whole  system  can  then  be  written  as 


£q  =  £o  +  <o  +  ^,  (13) 

where  Eq  is,  as  above,  the  ground  state  energy  of  the  4 He  background  system,  Co  is  the  binding  en¬ 
ergy  of  an  impurity  with  zero  momentum,  and  mu  is  the  so-called  hydrodynamic  effective  mass.  It 
is  the  contribute  to  the  effective  mass  arising  from  the  interaction  of  the  imj  *Hy  particle  with  the 
background. 

Similarly,  the  variational  wave  function  for  two  impurity  Fermions  with  coordinates  r{,  rj  is 


*#+2(rn  r?i  ri>  *  •  -rAr)  =  «p  ui(r0  +  "1(^2)  +  «2f(r{  >r2)+ 


+  23  lu2B(ri*rO+«2B(r2»r.O]>*o(ri,--MrAO^(l»2)i 


where  4>(1,2)  is  a  2-particle  Slater  determinant.  We  leave  out  the  backflow  correlations  since  we  will 
consider  only  small  momenta  of  the  3He  impurities.  The  only  new  unknown  function  is  the  impurity- 
impurity  correlation  function  Uo;(r{»r2)>  which  is  again  determined  by  minimization  of  the  total  energy. 

Instead  of  describing  the  further  manipulations  in  detail  (see  Refs.  15-17),  let  us  turn  to  the 
results  of  our  calculations.  The  solid  lines  in  Fig.  1  show  the  impurity  density  in  comparison  with  the 
background  density.  We  see  that  the  3Ile  particle  is,  as  predicted,  inside  the  outermost  4He  layer  for 
our  calculation  of  a  double-layer  system.  As  the  thickness  of  the  background  increases,  the  3He  particle 
is  pushed  outward  into  the  low-density  regime  of  the  film. 

A  more  instructive  pic*  ^re  is  given  by  considering  the  distribution  of  4lle  particles  in  the  vicinity  of 
the  impurity.  Figs.  4  and  5  show  the  A He  density  assuming  that  the  3He  atom  is  located  at  (xy  y)  =:  (0, 0), 
i.e.  the  quantity 

PS(^r\\)~  J (15) 

where  pib{y!  ,  r)  is  the  impuntv-background  two-booy  density.  For  the  very  thin  film  we  see  that  the 
3 He  atom  is  located  within  the  outermost  layer,  "here^s  for  the  thick  film,  the  4He  background  is  only 
modestly  deformed. 

We  have  earned  out  extensive  calculations  for  the  static  ground  state  of  one  3He  impurity  on  a 
family  of  4He  backgrounds  ranging  from  a  double  layer  (n  =  0.12A  ~)  to  a  system  of  about  five  helium 
layers  (n  =  0.35A  2).  The  hydrodynamic  effective  mass  has  been  calculated  using  the  uniform  limit1 
for  the  Euler-equation  for  the  backflow  function  a^r7,!*,).  In  this  limit,  one  obtains  a  closed-form 
expression  for  the  effective  mass  which  is  identical  to  a  second-order  approximation  for  the  self-energy 
in  terms  of  the  emission  and  re-adsorption  of  ripplons  and  phonons.  Our  theoretical  description  of  the 
first  layer  should  not  be  considered  very  realistic5,  but  the  uncertainty  in  the  description  of  the  first 
layer  does  not  seriously  alTect  our  conclusions.  This  part  of  the  system  has  a  very  high  density  and 
impurity  hardly  penetrates  (c.f.  Fig.  1).  Effectively,  the  impurity  atom  “sees”  the  first  lav' ,  i  1  wcf . 


Our  results  are  compared  in  Fig.  6  with  the  experimental  results  of  Refs.  11.  We  find  that 
the  hydrodynamic  mass  is  somewhat  lower  than  the  experimental  ma~s.  We  notice,  however.  that  the 


Fig.  4-  The  density  of  4He  atoms  in  the  vicinity  of  the 

3He  impurity,  pB(zy:  |j),  is  shown  for  the  dou¬ 
ble  layer  background  film  with  a  surface  cov¬ 
erage  of  u  =  0.1 5A  The  density  is  normal¬ 
ized  to  the  calculated  bulk  equilibrium  den¬ 
sity,  p™c. 


Ftg .  5.  Same  as  Fig.  4  for  the  background  films  with 
the  largest  surface  coverage  n  =  0.35 A  \ 


expenment  has  not  been  done  in  the  dilute  limit,  but  for  a  3IIe  density  of  0.3  layers.  Therefore  we  must 
include  the  effects  of  the  quasiparticle  interaction  between  the  3He  atoms. 

The  interaction  between  3He  impurities  is  intelligently  ’^cussed  in  terms  of  a  local  effective  inter¬ 
action  Vtjj( r,r').  This  effective  interaction  includes  both  the  direct  Van  der  Waals  force  between  the 
impurities  and  the  exchange  of  ripplons  and  phonons  through  the  4IIe  background.  It  can  be  obtained  b\ 
generalizing  the  theory  of  the  impurit) -impurity  interaction11  to  the  inhomogeneous  case.  In  the  case  of 
a  dilute  system  of  impurities,  we  may  take  the  impurity-impurin  interaction  as  an  effective  interaction 
between  any  two  3He  atoms  for  a  finite  impurity  population  ("quasiparticle  interaction”).  The  relation 
to  Landau  s  Fermi-Liquid  theory  in  two  dimensions18  is  drawn  in  momentum  space  by  identifying 

/cO?t|)  +  fc(<l\\)<r\  '  c-,  =  V{0+)  -  - V($,|)(l  4-  ci  tr2) 


J  dz}dz2d:ir^p[(zi)pI](z2)Vcjj(zi.zuii^)  1  -  , 


with  the  quasiparticle  interaction.  Given  the  interaction  (16),  wc  can  calculate  the  Fermi-Liquid  param¬ 
eters  of  the  two-dimensional  Fermi  liquid,18 


A,(0)/e(|qi  -  <12 1)  =  Fm  COs(»?ip) 

m=°  (1) 

OO  '  ' 

W(0)/<,(|qi  -  qal)  =  F>”  cos(7V<?)> 

m= 0 

where  ^  is  the  angle  between  <u  and  q2,  and  Ar(0)  =  the  density  of  states  at  the  Fermi  surface 


*  1'ioiTi  these,  we  obtain  Ur  magnetic  MistepUbiJUN  in  units  oi  tn<  susceputulHA  ol  llu  l r^t  two- 

\ 

dimensional  Fermi  gas,  ;*30, 

X'(0)/X’3o  =  (»>i//»»3)(l  4  F,72)/(l  4  Fq)  (18) 

and  the  total  effective  mass 

m*  =  m//(l  +  F{ /2).  (19) 

Figuie  6  (solid  line)  shows  the  total  effective  mass  r?r  calculated  from  Eq.  (19).  For  one  active  )a\er 
of  4He  the  hydrodynamic  mass  is  about  1.8  m3}  in  reasonable  agreement  with  both  the  data  of  Refs. 
11  and  the  conclusion  of  Ref.  12  with  their  thinnest  4  He  film.  In  this  regime  the  correction  due  to  the 
quasiparticle  interaction  becomes  quite  sizable  wLh  increasing  3Ile  concentration,  but  it  is  difficult  to 
make  a  quantitative  statement  due  to  both  the  high  density  of  the  background  and  the  rapid  variation 
of  the  effective  mass. 

Fig.  6:  The  calculated  effective  mass 
m*  is  shown  in  units  of  the  bare  3Ile  2 
mass  7713  as  a  function  of  the  inverse 
film  Ihickness  D~l  (dashed  line).  The 
circles  with  error  bars  are  the  experi¬ 
mental  data  of  Refs.  11.  The  upper 

S 

scale  shows  the  surface  coverage  corre-  \  1 

c 

sponding  to  the  film  thickness  D.  The 
dashed  line  shows  the  hydrodynamic 
effective  mass  mu  for  0.3  layers  of  3He. 

d"  (n 

For  all  cases  with  thicker  4IIe  films,  the  accuracy  of  our  theoretical  prediction  is  quite  satisfactory. 
The  microscopic  theory  predicts  a  hydrodvnamic  mass  that  is  consistent  with  the  estimate  m^/m^  = 
1.26  ±  0.15  given  by  Valles  et  al.  (Ref.  12),  but  below  the  results  of  Refs.  11.  The  agreement  with 
the  latter  data  is  improved  when  the  corrections  due  to  the  quasiparticle  interaction  are  included.  The 
quasiparticle  interaction  between  3Ile  atoms  in  the  surface  gives  a  density-dependent  correction  to  the 
effective  mass.  The  contribution  is  about  10%  for  the  case  of  0.3  layers  of  3Ile. 

There  are  slight  oscillations  of  the  3He  effective  mass  as  a  function  of  the  thickness  of  the  underlying 
film.  A  weak  plateau  is  seen  around  D~}  %  O.lA  This  coincides  with  the  regime  where  the  third 
liquid  layer  is  formed.  Thus,  we  conclude  that  the  impurity  effective  mass  has  weak  la\er  structure,  but 
that  these  oscillations  dar  k  i  out  rather  fast.  The  effect  of  the  migration  of  the  impurity  atom  into  the 
surface  is  much  stronger  than  the  modification  of  the  backflow  for  different  degrees  of  filling  of  the  last 
layer  so  that  it  would  be  much  harder  to  observe  fluctuations  of  the  3Hc  effective  mass  for  very  thick 
films. 

A  comparison  of  the  theoretical  and  the  experimental  magnetic  susceptibility  ratio  ,\(0)/A3o  is  shown 
in  Fig.  7.  He*1  we  have  taken  a  hydro*4  , namic  mass  of  1.2G77)3,  which  gives  the  best  fit  to  experimental 


refer  to  a  much  smaller I *  3IIc  coverage  \Yc  find  a  quite  satisfactory  agreement  berween  theory  and 
experiment. 

The  agreement  is  not  quite  so  good  for  smaller  4 * *Ile  coverages.  This  is  partly  due  to  the  fact  that 
no  attempt  was  made  to  re-adjust  the  effective  mass.  One  must  also  be  concerend  about  the  accuracy 
of  the  theoretical  description  when  the  3He  impurity  penetrates  into  a  regime  of  higher  4He  density, 
where  HNC  methods  are  intrinsically  less  accurate.  However,  HNC-typc  microscopic  theories  are,  while 
using  modest  computational  resources,  quite  capable  of  giving  at  least  a  semi-quantitative  picture  of  the 
structure  of  complicated  systems  like  the  ones  discussed  here. 

Finally,  we  stress  that  the  good  agreement  between  the  experimental  results  and  our  calculation 
has  been  obtained  with  a  siaiic,  momentum-dependent  quasiparticle  interaction.  The  calculation  of 
the  magnetic  susceptibility  involves  phase-space  integrals  of  the  quasiparticle  interaction  over  momenta 
between  <7  =  0  and  q  =  2 kp.  Therefore,  with  increasing  impurity  density  one  probes  the  momentum- 
dependence  of  the  quasiparticle  interaction.  We  believe  that  measurements  of  the  type  of  Ref.  12  are 
extemely  useful  to  enhance  our  theoretical  understanding  of  the  quasiparticle  interaction  in  liquid  3He. 

Fig .  7:  Magnetic  suscepti¬ 

bility  x/X30  of  the  3IIe  film  as 
a  function  of  areal  density  n  in 
atoms/A2,  for  4lle  coverages  between 
0.15  atoms/A2  (uppermost  curve), 
and  0.30  atoms/A2  (lowest  curve). 

1  Solid  squares:  experimental  data12 
for  a  surface  coverage  of  9.5  4 He  lay* 
ers.  Circles:  experimental  data12  for 
a  3IIe  density  of  0.088  layers,  for  4He 
coverages  of  2.8,  3.4,  4.75,  and  5.24 
atoms/A2. 


Layers 


I  hope  that  I  have  succeeded  in  this  contribution  to  present  the  study  of  quantum  liquid  surfaces  as 
an  exciting  field  foi  both  experimental  theoretical  research.  Let  me  conclude  by  highlighting  the  most 
interesting  physical  and  conceptual  aspects  of  this  field  of  research: 

♦  Physically,  one  is  able  to  construct  very  clean  two-dimensional  Fermi  systems  and  investigate  these 
systems  over  a  wide  regime  of  densities.  This  is  in  contrast  to  bulk  3Hc,  where  experiments  can  be 
performed  only  in  a  narrow  regime  around  the  equilibrium  density.  B\  carrying  out  measurements  on 

3He  impurity  films  over  a  large  density  regime,  one  should  be  able  to  measure  the  momentum  dependence 

of  the  quasiparticic  interaction. 

•  Conceptually,  quantum  liquid  films  are  many-bodv  systems  where  state-of-the-art  manv-body  theory 

is  necessary  for  a  thorough  microscopic  understanding.  At  this  time  I  do  not  intend  to  suggest  that 


quantum  liquid  films  be  studied  with  other  many-body  methods.  But  looking  back  at  the  obvious 
requirements  that  must  be  satisfied  by  a  theory  in  order  to  successfully  deal  with  such  systems,  it  may 
be  worth  asking 


Can  your  favorite  many-body  theory  deal  (in  prin¬ 
ciple  and  in  practice)  with  these  systems  ? 
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ABSTRACT 

The  effective  Hamiltonian  for  a  pair  of  Fermion  derived  from  the 
basic  many  body  Hamiltonian  allows  for  the  interaction  of  a  pair  with 
its  medium  in  a  number  of  ways  e.g.,  once  via  the  mean  field  and  another 
via  the  pairing  part  of  the  interaction.  The  two  state  model  developed 
earlier  has  been  solved  in  a  better  approximation.  The  model  has  then 
been  applied  to  calculate  pairing  energy,  critical  temperature, 
coherence  length  and  critical  field  for  with  X  =  0,  S,  Cl, 

F,  Br,  I.  The  model  can  account  for  the  general  trend  and  magnitude  of 
these  observables,  particularly  the  very  large  value  of  the  critical 
field.  Using  a  model  Hamiltonian  of  Frdlich  type,  the 
electron-phonon  coupling  strength  is  estimated.  The  electron- phonon 
coupling  strength  for  our  particular  system  is  more  than  an  order  of 
magnitude  larger  than  the  one  in  normal  superconductors.  The  Eliashberg 
dimensionless  coupling  parameter  X  is  found  to  be  about  2.0  to  3*0. 


1 .  INTRODUCTION 


Transition  from  a  normal  state,  be  it  a  conductor  or  a 
semiconductor  or  a  perovskite,  to  a  superconducting  state  involves  a 
phase  transition.  Whereas,  electrons  in  a  normal  state  behave  like 
incoherent  Fermions,  the  formation  of  3oson-like  coherent  pairs  is  a 
necessary  criterion  for  a  superconducting  state.  Near  critical 
temperature,  Tc,  one  has,  therefore,  two  types  of  states  of  a  system, 
one  representing  incoherent  Fermion  and  the  other  coherent  quasi  boson 
like  pairs.  A  set  of  coupled  equation  commensurate  with  this  situation 
has  been  derived”^  and  the  zeroeth  order  solution  has  been  successful 
in  predicting  the  trend  of  TP  in  Y-Ba-Cu-O-X  with  X  =  0,  S,  Cl  and  F. 


*A  travel  grant  from  the  U.  S.  Army  Research  Office  is  thankfully 
acknowledged. 

•'/ Dedicated  to  Professor  Don  Lichtenberg  in  appreciation  of  his  human 
qualities  and  contribution  to  physics. 


The  key  to  the  pairing  has  been  assumed  to  be  the  formation  of  doubly 
ionized  negative  ion  X^-*  Qualitatively,  one,  therefore,  expects  that 
the  substitution  by  Br  and  I  both  of  which  form  negative  ions  should 
give  rise  to  superconducting  behavior  and  this  has  recently  been 
confirmed3).  In  sections  2,  we  present  a  better  approximation  of  the 
two  state  model  and  calculate  Tc,  energy  gap,  coherence  length  and  the 
critical  field.  In  section  3  we  discuss  the  relation  of  the  calculated 
energy  gap  to  electron-phonon  coupling  constant  using  an  interaction 
similar  to  the  one  proposed  by  Frdlich.  Hamiltionian  containing  this 
interaction  has  been  diagonalized  using  the  method  of  Huang  and  Rys^. 

We  estimate  the  electron- phonon  coupling  strength  in  a  simple  model  and 
present  there  a  discussion  of  the  Eliasberg  coupling  constant. 


2.  THE  THEORY 

As  noted  in  ref.  1,  the  equation  for  a  pair  of  Fermion  in  state  a 
in  a  medium  is  given  by 


[e(i)  +  0(j)  +  v(i,j)  -  e^]  <i>a(i,j) 


+(1/2)1  J 

8=1 


dkdl<j>  (k ,1 )  W<j>  (i,j)$Q(k,l)  =  0 

P  U  D 


(D 


where  V(i,j)  is  the  effective  two-electron  interaction  obtained  by 
integrating  over  all  lattice  coordinates  R,  of  the  bare 
electron-electron  and  electron-lattice  interaction  hee  and  he_x, 
respectively. 


V(i,j)  =  <  hee(ij  )  +  he_x(R ,  ij  )>R  (2) 

0(i)  is  the  effective  one  body  operator  which  could  be  bare  kinetic 
energy  or  an  effective  one-body  interaction.  W  is  given  by 

W  =  [ V( i  ,k)  (l-X(i.k))  +V(i,l)  (1-X(i,l)) 

+  v(j,k)  (1-X(j,k))  +  V(j,l)  ( l-X(j  ,1)]  (3) 

X(i,j)  is  the  exchange  operator  transposing  all  coordinates  i  and  j.  <p 
is  an  antisymmetric  function.  W  is  the  average  potential  generated  by 
pairs  other  than  4>a.  In  case  all  the  pair  functions  are  replaced  by  2x2 
determinants  of  single  particle  orbitals,  (1)  reduces  to  the  set  of 
Hartree-Fock  equation.  Above  critical  temperature  the  total  electronic 
wave  function  <jjp  is  given  by  properly  ordered  permutation  P  of  of 
Fermion  type  of  pair  function 


(J;  =  (1//2n!)  I  eP$  (12)  . b  (2n-1,2n)  (*0 

P  a1  an 

Below  critical  temperature  <J»S  is  a  linear  combination  of  Schafroth  pair 
functions*^)  and  leads  to  condensation  necessary  for  superconductivity 

*  -  c  [(-1)P)[<J.(12)4.(3'0 

-  p 


4*(2n-1,2n)] 


(5) 


Because  of  the  difference  between  these  two  types  of  wavef unctions, 
the  average  potential  seen  by  a  pair  i.e.,  the  terms  containing  W  in 
(1),  differs  significantly  in  two  cases.  That  might  result  in  a  change 
in  total  energy  signaling  a  phase  transition.  Near  critical  temperature 
we  expect  significant  components  of  both  types  of  states  to  be  present 
and  the  total  wavefunction  may  be  written  as 


X  =  (6) 

Thus,  near  critical  temperature,  one  is  to  diagonalize  the  matrix 

•/  «n  h12^ 

^  h21  h22  j  (7) 

where  H ^ i ( i=1 ,2)  contains  diagonal  matrix  elements  in  states  1, 
identified  as  Fermion  like  states  and  in  state  2,  identified  as 
quasi-Boson  like  states.  2  =  H21  are  the  coupling  terms.  H-^  is  a 
sum  of  two  terms,  e ^  containing  effective  average  field  seen  by  a  pair 
and  ejxi  the  interaction  between  a  pair.  T.he  coupling  term  H^j  contains 
at  least  the  effective  two  body  interaction  V ( i , j )  but  might  also  take 
into  account  the  difference  in  mean  field  in  two  cases. 


The  diagonal ization  of  (7)  leads  to 


(8) 


E_,  being  lower  in  energy,  is  to  be  identified  with  the  superconducting 
state  and  E+  with  the  normal  state.  Since  (Hu  -  H22)2  and  |Hi2|2 
are, respectively ,  of  the  order  of  a  few  eV  and  a  few  tens  of  meV,  it  is 
reasonable  to  assume  jH-j  2|  2<<  | H -|  i-H22i  2.  Then 


E_  =  H22 


H11-H22 


ere2+e11-G22 


G11"G22  usually  much  less  than  (ei~e2).  Therefore, 


E  s  e2+e22  - 

<-.cc  epe2 


(10) 


.♦.  Pairing  energy  in  the  superconducting  state  AE 

AE  =  (E  -  e2)  =  e22  -  (11) 

Both  e22  jH12j2  involve  integration  over  all  lattice  coordinates  and 
represent  a  coupling  between  electrons  and  the  medium.  These  depend  on 
various  parameters  of  media  or  lattice  and  are  collectively  denoted  by 
vc. 


Since  this  is  a  zero-temperature  theory  one  can  relate  AE  to  Tc 
using  the  relation7* ^ 

(ir/Y)  kTc  =  AE  (12) 

(Here  k  and  Y  are,  respectively,  Boltzmann  and  the  exponential  of  Euler 
constants  and  n/Y  1.761!).  Assuming  |Hi2|2  to  depend  on  mass  of  ion  X 
as  1//M,  we  may  write 


a(<)  -  b(<)//M(e2-ei ) 
*  a(<)  -  bU)  AW 


(13a) 


T 


c  ” 


(13b) 


Using  Tc  °  93°K  and  90°K  for  X  =  0  and  S,  respectively,  we  find 
a(>0“310K  and  b(«)  =  2521  °K.  Since  AW  for  0  and  S  are  very  close,  we 
expect  them  to  have  similar  Tc  which  is,  indeed,  the  case.  As  noted  in 
refs.  1  and  2,  the  superconducting  state  in  this  system  is  identified 
with  the  localization  of  electron  pairs  in  forming  doubly  ionized 
negative  ion,  X2“  and  the  normal' state  with  singly  ionized  negative  ion, 
X^-.  In  table  1  we  have  compared  calculated  Tc  with  the  observed  ones. 
For  X  «  Br  and  I,  we  provide  theoretical  estimates  using  a  range  ol 
electron  affinities  of  Br-  and  I-.  The  theory  can  reproduce  the  trend. 

There  is  some  uncertainty  about  the  validity  of  (12)  in  this  type 
of  superconductors.  Nevertheless,  we  have  used  it  to  calculate  the  gap 
parameter  and  tabulated  them  in  table  2.  The  measurements  for  the 
oxygen  case  is  in  accord  with  the  calculated  number  but  indicate  that 
the  constant  of  proportional ity  between  A  and  Tc  is  about  twice  the  one 
in  (12). 

Coherence  length  t,Q  is  related  gap  parameter  by 

C0  =  nVp/TtA  (I'D 

where  Vp  is  the  velocity  near  Fermi  surface..  Assuming  the  Vp  in  this 
type  of  superconductors  is  close  to  those  in  type  II  superconductors,  we 
may  determine  this  from  Nb^SnCTc  =  18°K,  c0  =  50  A).  Calculated  values 
are  tabulated  and  compared  with  data  in  table  2. 

Since  A  and  c0  for  this  type  of  superconductors  are,  respectively, 
substantially  larger  and  smaller  than  those  in  type  I  and  II 
superconductors ,  we  may  call  these  type  III  superconductors.  The  large 
A  and  small  coherence  length  are  indicative  of  very  strong  coupling 
between  pairs  and  the  environment  and  the  coherence  is  more  localized. 
Thus,  the  BCS  approximation  of  a  constant  weak  coupling  between  a  pair 
of  electron  via  phonons  is  invalid  for  these  materials. 

The  difference  between  the  total  energy  of  the  normal  and 
superconducting  states  is  related  to  critical  field  Hc  by  Hc2/87:  = 
energy  difference  per  unit  volume 

H^/8ti  =  x  /(H^_H22)^  +  (15) 

where  p  is  the  number  of  electrons  per  unit  volume.  We  have  used  (8)  to 
get  the  energy  difference  per  pair  between  the  normal  and  the 
superconducting  state. 

Assuming  that  (en~e22)  <<  we  8et 

(Hc2/8tt)  (2/p)  s  |(ere2)  +  Hj^-  | 

G 1  G2 

=  |  (ei-e2)"2(7i/Y)k  b(<)  AW J 


(16) 


The  last  approximation  is  valid  since  the  second  term  is  of  the  order  of 
a  few  tens  of  meV  but  |  ( e ■(— e 2 )  |  is  about  5  to  lOeV.  Eq.  (16)  is  an 
important  result  and  is  very  different  from  the  expression  obtained  in 
the  Schafroth  or  BCS  theory.  In  latter  two  cases  instead  of  (15),  one 
has 9  (Hc2/8u)  (2V/pg)  =  A^.  Thus,  the  theory  predicts  that  critical 
field  in  this  types  of  superconductors  is  much  larger  than  that  in  type 
I  superconductors.  We  may  estimate  p  from  normal  superconductors  e.g., 
tin  which  has  a  Tc  -  3.73°K  and  Hc  =  306  Gauss.  Noting  the  p  per  meV  = 
P[r/V,  we  get  p  -  lO^/cmS  from  the  tin  data.  Calculated  critical  field 
using  this,  is  noted  and  compared  to  the  observed  value  in  table  2.  The 
theory  can  account  for  the  large  critical  field  seen  in  high  Tc 
superconductors. 


3.  EFFECTIVE  COUPLING  STRENGTH  IN  A  FROHLICH  MODEL 

Describing  the  medium  and  the  electrons,  respectively,  by  Boson 
operators,  bk  and  Fermion  operators ak  we  may  write  the  Hamiltonian  for  a 
linear  coupling  between  the  two  as  follows 

H  =  ^AeAaIaA  +  ^  ha,(^  Cbkbk  + 

k  ^ 

+  lx  £  a*a^  [v(k,A)bk  +  v(k,A)bk]  (17) 

This  is  basically  the  Hamiltonian  proposed  by  Frdhlich1*  and  used  in 
the  Eliashberg11  formalism.  Taking  v(k,A)  to  be  the  matrix  element  of 
(c.VV)  where  e  and  V  are,  respectively,  unit  polarization  vector  and 
potential  at  lattice  sites,  one  may  define  a  dimensionless  coupling 
parameter  A  within  the  framework  of  Eliasberg  formalism.  The  averaged 
electron-phonon  interaction  a^(w) ,  averaged  electron- phonon  matrix 
elements  <g2>  and  the  dimensionless  coupling  constant  A  are  related  as 
follows: 


2(  w)F(w)d<jj/w 


(18) 


N(0) 


»< g2>  _  f 

2M  0J 


ua^(  w)F(  (o)dm 


(19) 


In  the  above  F(w)  is  phonon  density  of  states,  and  N(o)  is 
electronic  density  of  state  at  the  Fermi  energy  and  M  is  the  atomic  mass. 
In  case  =  constant,  as  suggested  by  McMillan11* 


-  N(0)<82>  _  2a?Lv1 

M  <<^>  "  2a  Jl  U/(u 


d0J 


(20) 


with  <u)2>  s  |F((ij)wdu/|F(oi))a)_1  dm 


(21) 


Evaluation  of  and  <g2>  require  a  knowledge  of  A  and  either  <cj2>  or 


<w"1>. 


A  can  be  determined  from  observed  T~.  For  small  and  moderate 


values  of  A,  McMillan's  relation  holds 


1H 


t  v-  avnr  l-OMl  +  X) 

Tc  =  K  exp[  x_mi((1+0<62A)-) 


(22) 


Here  p*  is  McMillan’  3  Coulomb  parameter  and  K  is  either  logarithmic 
average  phone  frequency  wx0g/1.2  or  Debye  temperatue  Td/1.^5. 

For  very  large  X  i.e.X>  five  times  (phonon  frequency)2,  Allen  and 
Dyne’s^ 3  relation  holds 

T c  -  0.180(1+2.60  u*)_1/2  X1/2  (23) 

Implication  of  (23)  is  that  there  is  no  upper  limit  of  Tc.  x,  is  not 
expected  to  be  very  large  for  the  1-2-3  superconductors.  Since  wiog  is 
not  very  well  known,  we  can  estimate  X  using  K  =  Tq/  1 .  -M5.  For  X  =  oxygen 
Tc  =  93.  and  Tp  =  812  (ref.  1 6) .  Calculated  values  for  X  are  1.38, 

1 .89,  2.09,  2.H5  and  2.81  for  p*  =  0,  0.10,  0.13,  0. 1 8  and  0.22, 
respectively.  Thus,  X  is  estimated  to  be  between  2  and  3.  Both  a2  and 
<g2>  can  be  evaluated  once  phonon  spectra  become  available. 

One  may,  however,  follow  a  different  approach  noting  that  the 
Hamiltonian  (17)  can  be  diagonalized  exactly  following  the  method  of 
Huang  and  Rys5*!2.  xne  expression  for  diagonalized  energy  E  is  given 
by 


E^(n,k)  =  -fiuKk)  [nk  +  (2^) 

a  k 

with 

AX  =  IRlv('<,X)12ffiu)(k)  '  (25) 

v(i?,X)  could  take  various  forms,  including  the  matrix  element  of  (e.VV). 
We  can,  however,  make  a  simple  estimation  of  the  electron- phonon 
coupling  strength  by  setting  v(it,X)  =  g-ftwtk).  In  this  case2^ 


Aa  =  8a  Ikfcw(k)  8 \  k  tD 


(26) 


Using  the  observed  values  of  A  =  32  mev!5  and  50  meV1?  and  a  Debye 
temperature^ 6  Tp  =  812°  K,  we  get  g2  s  0.^6  and  0.71,  respectively. 


Both  (23)  and  (25)  in  conjunction  with  (12)  imply  that  there  is,  in 
principle,  no  limit  of  Tc.  The  key  to  high  Tc  is  strong  electron- phonon 
coupling  as  suggested  by  Frdhlich^. 


Since  the  model  should  be  valid  for  type  I  and  II  superconductors, 
we  can  calculate  g2  for  those  cases.  In  Table  3,  we  have  compared 
effective  coupling  strength  g2  for  type  I,  II  and  III  superconductors. 
The  electron- phonon  coupling  in  high  Tc  superconductors  is  two  orders  of 
magnitude  larger  than  that  in  type  I  superconductors.  It  is  about  an 
order  of  magnitude  larger  than  that  in  type  II  superconductors. 


1J.  CONCLUDING  REMARKS 

Although  in  this  analysis,  particularly  in  section  3,  one  has 
mentioned  electron- phonon  coupling,  the  theory  outlined  in  section  2 
doe3  not  restrict  superconductivity  to  effective  pairing  caused  by 
electron- phonon  coupling  case  only.  IN  FACT,  ANY  MECHANISM  THAT  MAKES 
EFFECTIVE  INTERACTION  (2)  ATTRACTIVE  AND  ALLOWS  TO  FORM  SCHAFROTH'S  PAIR 


COULD  GIVE  RISE  TO  SUPERCONDUCTIVITY.  A  coupled  electron  pair  in  a 
singlet  state  and  with  opposite  momenta  is  usually  the  lowest  energy 
state  in  many  potential  wells.  Hence,  it  is  natural  to  use  this 
coupling  scheme  along  with  the  coherence  state  condition  for 
superconducting  states.  One  can  form  such  a  state  in  many  ways  e.g.,  by 
surrounding  a  pair  of  electrons  by  positively  charge  ions.  Given  the 
right  condition,  that  system  could  also  exhibit  superconducting 
behavior,  even  though  no  electron- phonon  coupling  is  directly  involved. 
There  is  no  upper  limit  of  Tc  an.d  a  strong  electron-phonon  or 
electron-environment  interaction  is  a  necessary  requirement  for  high  Tc, 
as  noted  by  Frdhlich^). 
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Table  1:  Calculated  Tc  (col.  5),  are  compared  to  observed  one  (col.  6). 
Electron  affinities  in  col.  2  and  for  0  and  S  in  col.  3  are  taken  from 
ref.  18.  Electron  affinities  for  F"  and  Cl“  are  theoretical  calculation 
and  Br  and  I  are  estimation,  (a),  (b),  (c),  (d),  (e)  and  (f)  are  refs. 

1 9— 23  and  3,  respectively. 


X 

e.a(x)(eV) 

ea(x~)eV 

AW 

Tr,(cal.)0K 

To(expt)0K 

0 

1.47 

-8.73 

-0.02451 

93 

93a 

F 

3-45 

-4.81 

-0.02778 

102 

s 

2.07 

-0.02332 

90 

90e 

Cl 

3.61 

mammm 

-0.01739 

75 

72e 

Br 

3-36 

-0.01202 

-0.01346 

-0.01528 

61 

65 

70 

-  75- 

T 

3.06 

S£k3 

i 

-0.00979 

-0.01100 

-0.01259 

56 

59 

63 

-  50 1 

Table  2:  Calcul  ted  gap  A  (col. 2),  coherence  length  c,0  (col.  4), 
critical  field  (col.  6)  are  compared  with  respective  observed  values,  in 
columns  3,  5  and  6. (a),  (b),  (c),  and  (d)  are,  respectively,  ref.  24-27. 


X 

A(cal.) 

MeV 

A(expt) 

MeV 

CQ(calc. )A 

C0(exp)A 

_ j: _ 

H^(cal.)G 

Hc(exp)0e 

0 

14 

10 

m 

II 

F 

16 

- 

9 

- 

- 

S 

14 

- 

10 

- 

- 

- 

Cjj 

11 

- 

12 

- 

- 

- 

B>  ■'-11 

- 

13-15 

- 

- 

- 

I  10 

- 

15-17 

- 

- 

Table  3:  Coupling  strength  g2  in  some  type  I,  II  and  III 
superconductors . 


Substance  -*■ 

TV  in  °K 

k  x  debye  Temperature  (meV) 

B2 

■ 

3 


Re 

1.69 

35.8 

0.004 


Sn 

3-72 

17.20 

0.01 


Nb 

9.22 

23.90 


9 


0.03n 


y-Ba-Cu-0 

93.00 _ 

70.06 _ 

0.46  to  0.71 


A1 

1.16 

36.9 
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INTRODUCTION 

In  the  last  workshop  in  this  series,  I  presented  some  results  obtained  by  ex¬ 
pressing  the  Baym-Kadanoff  algorithm[l-2]  in  the  language  of  the  parquet  theory[3- 
14).  That  paper[15]  showed  (i)  that  the  parquet  theory  is  not  a  conserving  theory 
in  the  sense  of  the  Baym-Kadanoff  theory,  (ii)  applied  the  Baym-Kadanoff  theory 
in  this  form  to  some  simple  self-energy  diagrams,  and  <  Ill  •  observed  that  none  of 
the  resulting  vertex  functions  was  antisymmetric.  This  contribution  presents  some 
new  results  obtained  using  the  same  formalism.  The  main  result  is  a  proof  that  any 
set  of  diagrams  which  is  both  <  unserving  and  antisymmetric  under  exchange  of  the 
outgoing  legs  must  include  all  parquet  diagrams  obtained  using  the  bare  interaction 
as  well  as  an  infinite  number  of  irreducible  diagrams. 

NOTATION 

The  notation  has  been  changed  slightly  since  the  last  workshop,  so  a  definition 
of  terms  is  in  order.  The  diagrammatic  constructs  of  interest  are  the  two-body 
Green’s  function  and  the  single-particle  self-energy.  Diagrams  for  the  two-body 
Green’s  function  are  to  be  written  in  terms  of  one-body  Green’s  functions  which 
use  the  self-energy;  hence  explicit  inclusion  of  self-energy  insertions  in  the  diagrams 
for  the  two-body  Green's  function  should  be  avoided.  The  lowest-order  contribution 
to  the  two-body  vertex  is  taken  to  be  the  bare  interaction  and  its  exchange;  the 
terms  in  the  two-body  Green’s  function  which  are  just  the  direct  and  exchange 
products  of  two  one-body  Green’s  functions  are  not  included  in  this  category.  A 
direct  diagram  is  one  in  which  the  line  of  propagators  entering  one  side  at  the 
bottom  ultimately  can  be  traced  through  to  the  same  side  as  it  goes  out  the  top. 


equationsjS 

r  =  V  +  S  4  V  -  C  +  L  +  R 
5  =  r©(V  +  V  -hC  +  L-rR) 

U  =  r®(V  +  S  +  C  +  L  +  R) 

(1) 

c  =  r©(v  +  5+i/) 

L  =  T(D{v  +  S  +  U) 

R  =  T<I)(V  +  S  +  C  +  L  +  R). 

The  apparent  asymmetry  in  these  equations  arises  from  the  fact  that  some  diagrams 
generated  by  the  ©.  ©,  and  (D  operations  could  be  generated  in  several  different 
ways.  The  equations,  as  written,  guarantee  that  each  such  diagram  will  be  generated 
but  once.  A  more  symmetric  set  of  equations,  which  is  equivalent,  is  written  in  the 
form  [12,14,16] 

t=i+s+u+t 

S  =  (T-S)®T 
U  =  (T-U)®T 

T  =  B  +  T(D(I  +S  +  TJ  +  B)  +  {I  +  S  +  U  +  J9)©r  +  T©(7  +  S  +  U  +  B)®T 
B  =  (J  +  S  +  U)(B>(I  +  5  +  V  +  B), 

(2) 

where 

a©/?  =  a©0  +  a©r©/3.  (3) 

Each  form  has  its  utility;  eq.  2  is  more  symmetric,  but  it  also  has  two  operations 
explicitly  on  the  right-hand  side. 

A  convenient  way  to  compute  the  single-particle  self-energy  E  is  to  take  some 
approximation  for  the  exact  vertex  I\  add  an  interaction  V*  using  the  operation  ©, 
and  then  close  off  the  two  right-hand  lines  with  a  one-body  propagator.  If  T  is 
antisymmetric,  the  self-energy  can  be  computed  from  the  direct  part  of  T  by 

£  =  d(r©V')+e(r©V),  (4) 

where  the  direct  and  exchange  lines  contributions  are  formed  by  applying  the  two 
operations  illustrated  in  fig.  2. 

If  T  is  i\ot  antisymmetric,  then 

E  =  d(t®V).  (5) 

The  exact  theory  and  the  parquet  theory  both  have  an  antisymmetric  T,  so 
that  r  =  T  and  eq.  4  is  appropriate.  However,  the  Baym-Kadanoff  algorithm 
generates  an  approximate  vertex  which  is  not  necessarily  antisymmetric. 


BAYM-KADANOFF  THEORY 

The  Baym-Kadanoff  (BK)  algorithm  takes  an  approximation  to  the  self-energy 
S  which  is  taken  to  be  constructed  diagrammatically  from  one-body  Green’s  func¬ 
tions  which  use  the  same  E.  The  source  of  the  approximate  E  need  not  be  speci¬ 
fied;  it  could  be  computed  from  some  approximate  vertex  using  eqs.  4  or  5,  but  it 
need  not  be.  The  conserving  vertex  is  obtained  as  the  functional  derivative  of  the 
self-energy  with  respect  to  a  weak  external  two-point  non-local  perturbation.  The 
one-body  Green’s  function  is  itself  a  functional  of  this  interaction,  and  the  func¬ 
tional  derivative  of  the  one-body  Green’s  function  can  be  expressed  back  in  terms 
of  the  functional  derivative  of  the  self-energy  with  respect  to  this  small  non-local 
perturbation.  We  note  that  any  self-energy  contribution  can  always  be  expressed  in 
the  form  of  fig.  2,  where  the  T  may  be  irreducible.  The  algorithm,  in  the  parquet 
notation,  proceeds  in  two  stages. 

First,  one  generates  a  two-body  kernel 


(6) 


where  the  prime  denotes  the  functional  derivative  with  respect  to  the  Green’s  func¬ 
tion.  Careful  inspection  shows  that  5  can  have  both  direct  and  exchange  parts,  but 
they  cannot,  be  of  all  reducibility  types.  Specifically,  only  reducibility  types  J,  S, 
U ,  /,  S\  6’,  L ,  R  and  A  can  come  out  of  this  operation.  To  show  this,  we  first  write 

d'(a®0)  =  d{a'®0)  +  d(a®0')  +  a®0  +  a©0  +  ^©0 

e'(a®0)  =  e(a'®0)  -f  e(a®0')  +  a®0  +  a®0  +  a©0 
The  immediate  functional  derivatives  are  all  of  the  types  indicated.  For  computing 
E,  0  —  V,  so  that  0'  =  0.  The  only  remaining  terms  are  the  d(a'®0 )  or  e(a'®0). 
If  a  is  irreducible,  then  these  terms  are  as  well.  If  a  is  reducible,  then  the  functioned 
derivative  can  be  transformed  using 

d((a®0)'®  7)  =  d((a'®0)®7)  +  d.((a®'0)®7)  +  d((a®0')®7) 
e((a®0)'®7)  =  e((a'®0)®7)  +  e((a®'0)®7 )  +  e((a®0')®7). 


(8) 


The  last  term  on  the  right-hand  sid**$  of  oq.  8  is  always  irreducible.  The  middle 
term  can  be  expressed  in  reducible  form  using  table  1.  All  of  these  table  entries  are 
in  one  of  the  stated  categories. 


©  d{{a®'/3)®  7)  d{(  a®'0)<§h) 


®  a©(/3® 7)  a©((3®~))  a®(3®~))  o®(d®y) 


©  <*©(7®/?)  a®(ft©’i)  a©{ 7©/?)  a©(d®y) 


©  a©(y©ft)  a®(y©{3)  ot®h<D0)  o.©(y®3) 


©  a©(7©/?)  a©( 7©/?)  a®( 7©/?)  a®(y©0) 


a©(7©/?)  a©( 7©/?)  q©(/?©7) 


Table  1:  Diagrams  from  the  middle  term  of  the  right-hand  sides  of  eq.  8. 


The  final  case  is  the  first  term  on  the  right-hand  sides  of  eq.  8.  These  can  be 
rearranged  using  the  diagrams  in  fig.  3  which  correspond  to  the  rigln-hand  sides  of 
the  identities 

d((a®(3)®  7)  =  d(a®((3® 7)) 
d((a@/?)©  7)  =  d(a©(/?®  7)) 
d((  a©/3)®y)  =  <£(a©(  7©/?)) 
d((a©/3)®y)  =  d(a®(  7©/?)) 
d((a©/?)©  7)  =  e(o'©(/?©7)) 
e((a©/3)©7)  =  e(a®{fi®y)) 
e((a©/3)®  7)  =  e(a©(/3©7)) 
e((a©/S)®y )  =  d(  a®(y®0)) 

e((or©/3)©7 )  =  e(a©(7®/3)) 
e((a(Z)/?)©7)  =  e(a©(/?©7)) 

These  diagrammatic  identies  are  equally  valid  if  a  is  replaced  by  a',  and  hence  the 
first,  terms  on  the  right-hand  sides  of  eq.  8  can  be  transformed  into  the  form  of  the 
left-hand  sides  of  eq.  8..  This  process  cannot  go  on  indefinitely,  since  each  time  eq.  9 
is  used  to  make  a  transformation,  the  diagram  a  is  of  lower  order  in  perturbation 
theory. 

The  full  BK  conserving  vertex  is  generated  from  E  by  solving  the  equation 

X  =  £  +  Z©X +E®X +±©X.  (10) 

Sorting  out  the  direct  and  exchange  parts  and  doing  some  diagrammatic  rearrange¬ 
ment.,  we  cast,  the  resulting  components  of  the  BK  vertex  in  the  form  of  eq.  2, 


Fig.  3  Some  diagrammatic,  rearrangements 
obtaining  exchange  components 

r* 

ix 
Sx 
Ux 

Tx 

and  direct  components 

r.v  =  Ix  4  Sx  4  Ux  +  Tx 
I  x  =  h 
Sx  =  SE 

Ux  =  Us  (12) 

Tx  =  Bx  4  t(D(Ix  4  S*  4  Ux  4  Bx)  4  (Ix  4  Sx  +  Ux  +  Bx)®t 
4  t(D(Ix  +  Sx  4-  Ux  +Bx)®t 
Bx  =  (Ix  +  Sx  +  Ux)®(Ix  4  Sx  4  Ux  4  Bx) 


=  Ix  4  Sx  4  Ux  4  T 
=  h 

=  Ss  (11) 

=  (tx  -  Vx)®tx 
=  TW 


The  ©  operation  in  eq.  5  requires  the  Lindhard  bubble  to  be  dressed  by  f\ 

These  expressions  for  the  conserving  vertex  basically  can  be  easily  understood: 
the  self-consistency  implicit  in  the  BK  algorithm  is  a  t-channel  operation  as  seen 
in  eq.  10.  Thus,  the  integral  equations  build  the  t-channel  chains  and  vertex  cor¬ 
rections,  while  the  exchange  channel  is  the  u-cliannel.  The  s-channel  diagrams  can 
only  come  from  the  explicit  functioned  differentiation  with  respect  to  the  one-body 


Green  V  function,  since  they  can’t  be  generated  by  a  direct  or  exchange  t-channel 
operation. 

As  we  noted  before,  the  vertex  generated  by  the  BK  algorithm  need  not  be 
antisymmetric;  indeed,  eq.  11  and  eq.  12  show  what  the  diagrammatic  content, 
actually  is. 

CONSERVING  AND  ANTISYMMETRIC  IS  MORE  THAN  PARQUET 

In  this  section,  we  prove  that  any  approximation  to  the  self-energy  which  is 
used  to  generate  a  vertex  which  is  both  antisymmetric  and  conserving  and  which 
includes  the  Hartree  term  must  include  a  class  of  diagrams  which  includes  all  par¬ 
quet  diagrams  generated  from  the  bare  interaction  and  a  large  class  of  irreducible 
diagrams.  The  key  to  this  argument  is  that  antisymmetry  of  the  conserving  ver¬ 
tex  implies  no  need  to  distinguish  between  T  and  f\  This  provides  a  significant 
constraint,  since  eq.  11  and  eq.  12  then  can  be  used  to  identify  certain  functional 
derivative  results  with  the  results  of  integral  equations  in  the  ©  and  ©  channels. 

At  the  Hartree-Fock  level  (Fock  is  included  if  Hartree  is),  I  will  include  V. 
Now  assume  that  all  of  the  parquet  terms  are  generated  up  to  order  n  in  perturbation 
theory  by  some  set  of  self-energy  diagrams.  (A  minimum  set  of  self-energy  diagrams 
can  always  be  found  by  using  the  algorithm  of  the  following  section).  The  exchange 
channels  (eq.  11)  then  generate  U  correctly  to  order  n  +  1  using  only  the  V  in  /, 
since  all  of  the  inputs  are  correct  to  order  n  and  the  first-order  term  is  correct.  The 
direct  channels  (eq.  12)  generate  T  correctly  to  order  n  +  1  for  the  same  reason. 
Finally,  we  note  that  the  functional  differentiation  of  any  diagram  which  produces 
or ®/3  also  produces  a©/?,  so  that  since  U  is  correct  to  order  n  -j- 1,  so  must  be  S. 
But  if  all  of  the  diagrams  are  correct  to  order  n  +  1,  then  they  are  correct  to  all 
orders.  Hence  all  parquet,  diagrams  must  be  generated. 

As  we  argued  previously,  parquet  is  not  conserving;  a  set  of  self-energy  di¬ 
agrams  which  generates  parquet-  must  also  generate  an  infinite  set  of  irreducible 
diagrams.  As  a  consequence,  parquet  theory  is  not  a  sufficiently  large  class  of  dia¬ 
grams  to  be  both  conserving  and  antisymmetric.  We  speculate  that  the  only  set  of 
diagrams  which  satisfies  both  criteria  is  the  complete  set  of  diagrams. 

PARQUET  DIAGRAMS  GENERATE  SELF- ENERGY  DIAGRAMS 

Here  we  present  the  algorithm  which  determines,  for  any  reducible  diagram  p., 
the  parquet  diagram  whose  direct  or  exchange  contribution  to  the  self-energy  gives 
fj  on  functional  differentiation  with  respect  to  the  Green’s  function.  It  is  convenient 
to  use  the  notation  of  eq.  1  in  this  section.  The  first,  step  is  to  determine  whether 
the  diagram  is  a  member  of  the  acceptable  types  {5,  £/,  S,  C,  L,  R}. 

If  so,  then  fi  =  axA,  and  the  starting  point  is  obtained  by  looking  in  tables  2- 
3. 

In  tables  2-3,  the  — >  entries  indicate  that  the  direct  diagram  can  be  rearranged, 
and  one  should  look  up  the  transformed  form. 


4  I 
» 


Table  2:  The  starting  point  for  direct  diagrams. 


A 

(g©A) 

(a©A) 

(a©  A) 

(a©A) 

V 

e(a©V*) 

d(a©F) 

e(a©U) 

e(a©F) 

(3®7 

-»  ((a©/?)©7) 

d((a®j3)®  7) 

e((a©/?)©7) 

e((a®/3)®j) 

(3©y 

e((a©7  )©/?) 

d{{a©y)®j3) 

e((a©7  )©/?) 

-4  ((a®/3)®7) 

8©  7 

d((a(Df3)®  j) 

-*  ((a©/?)© 7) 

-*  ((a©/?)© 7) 

d((a(D  7)®fl) 

/?©7 

e((a(D/?)©  7) 

->  ((a©/?)© 7) 

— >  ((a©/?)©7) 

e((a©7)®/3) 

/?©7 

e((a©8)®y) 

->  ((a® (3)© 7) 

— >  ((a®/?)© 7) 

e((a®7)®P) 

Table  3:  The  starting  point  for  exchange  diagrams. 


The  second  kind  of  entry  is  indicated  by  the  first  line:  if  A  is  V’,  then  we  are 
really  done,  since  then  a  is  the  parquet  diagram  whose  functional  derivative  gives 
the  desired  diagram. 

The  final  case  is  a  little  more  complicated:  the  table  entry  is  used  as  a  starting 
point  to  carry  out  a  reduction  using  table  4.  That  is,  a  direct  or  exchange  diagram 
from  table  2  or  3  is  used  to  look  up  a  new  diagram  from  table  4.  If  this  diagram 
has  the  form  d(u®V)  or  e(n©V),  the  procedure  terminates  with  v  as  the  parquet 
diagram  whose  functional  derivative  gives  the  desired  diagram.  If  not,  the  table 
entry  is  used  to  generate  a  new  table  entry.  For  parquet  diagrams,  this  procedure 
eventually  must  terminate. 

DISCUSSION 

We  are  at  present  at  work  determining  a  good  form  for  the  three-body  integral 
equations  which  sum  all  irreducible  diagrams  generated  by  the  BK  algorithm  using 
parquet  as  input.  The  present  form  of  the  equations  is  not  particularly  enlightening, 
and  will  not  be  presented  here. 

The  main  conclusion  to  be  drawn  is  that  two  properties  of  the  exact  vertex 
which  one  would  desparately  like  to  maintain,  antisymmetry  and  the  conserving 
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Table  4.  Iterate  using  this  table  until  the  right-hand  term  is  V. 

property,  are  incompatible  at  any  level  less  than  parquet.  This  has  been  pursued 
for  the  purpose  of  finding  out  a  satisfactory  approximation  scheme  for  fermion 
parquet,  theory.  It  appears  that  to  satisfy  these  properties  where  they  are  needed, 
it  will  be  necessary  to  use  different  approximations  in  different  contexts. 
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